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1. INTRODUCTION 


1.1 Problem 


Improving thermal efficiency and the specific performance by 
increasing a turbine’s inlet temperature and by increasing the 
compression ratio is an important element in the further development 
of gas turbines. In current engine designs it is possible to control 
compression ratios of approximately 30:1 and turbine inlet 
temperatures of almost 1,700 K Cl, 23. However, the goal of future 
developments are turbine inlet temperatures which are as close as 
passible to the stoichiometric combustion temperature (s 2,400 K) of 
the combustion chamber (see [1,2]). Standard values for an increase 
of the compression ratio are contained in latest thermo-dynamic 
circulation procedure studies performed by NASA C13. In accordance 
with these studies and according to the turbine’s inlet temperature, 
compression ratios of a 64:1 magnitude will be required in order to 
achieve minimum specific fuel consumption. 

Despite further development of available materials, cooling of 
components which are subjected to high thermal and mechanical stress 
will remain a necessity in the forseeable future. Therefore, the 
main task in achieving these goals is an optimization of the cooling 
requirement. By removing that portion of the exhaust which is 
necessary for cooling, the mass flow is reduced from the combustion 
chamber. This, in connection with further losses, results in a 
decrease of performance and thermal effeciency. Therefore, an 
appropriate design of the cooled components will avoid a situation 
where, with the intended additional temperature and pressure 
increases, an expected gain is compensated by an additional cooling 
air requirement. 
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Components with the highest cooling requirement are the stationary 
and rotating blades at the initial stages of a turbine. They are 
exposed to high thermal and mechanical stress. With these components, 
the cooling requirement and the necessary design o-f the cooling 
con-figuration on the inside o-f the blade depend on the development o-f 
heat transfer along the blade surface. To predict this heat transfer 
in a reliable manner is presently the main difficulty in the design 
and the further development of cooled gas turbine blades. Available 
calculation procedures are still very limited with respect to their 
efficiency. According to latest studies C 1 3 , available calculation 

methods can only be used for advance computation of blade surface 
temperatures with an exactness of i 100 K. This uncertainty in the 
theoretical determination of the outer heat transfer and the 
resulting blade temperatures would require an unacceptable 
overcooling of the blade to guarantee the required strength and life 
expectancy. In order to avoid this inefficiency, the calculation is, 
in a comprehensive design process, always supplemented by an 
experimental test on a model blade. This, however, requires 
considerable additional effort with respect to cost and time Eli. As 
a result, the strategy of modern gas turbine construction strives to 
reduce additional efforts for experiments by constant improvement of 
calculation methods. 

The causes for the above problems with blade design can be found in 
the complex and partially inexplicable heat transfer and flow 
procedures along the blade surface Cl, 23. In addition, a number of 
parameters and flow phenomena have a decisive influence on the 
boundary layer behavior and the heat transfer to the blade. The 
following parameters are of great significance: 

- the laminar-turbulent reversal behavior of the boundary layer 

- the blade’s longitudinal surface curvature 

- the highly negative pressure gradients in the flow 

- superimposed free stream turbulence 

- material value influences generated by the cooling of the 
blade, and 

- separation in the boundary layer. 
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Another significant -factor in blade flow of gas turbines is that the 
magnitude of the influence of these effects are always superimposed 
in a non-linear manner. Even when each parameter is studied on an 
individual basis, we presently have no full appreciation of respect 
its effect Cl, 23. Therefore, ignorance of basic knowledge and 

available data are the reasons for incomplete modelling of numeric 
boundary layer calculation procedures and the resulting uncertainties 
for blade design when the outer heat transfer is calculated. This 
was recently demonstrated in parallel studies for the development of 
differential equations to describe heat transfer on cooled gas 

turbine blades C43. 

For the development of calculation procedures, there is a basic lack 

of experimental studies, which provide information on the 

significance of the most important quantities: 

- free stream turbulence 

- free stream acceleration 

- and wall cooling 

and which take into consideration the mutual interaction of these 
parameters. 

This dissertation will meet the requirements for fundamental 
experimental studies of the influence of free stream turbulence, 
pressure gradients and wall cooling. Under turbine-like free stream 
conditions, comprehensive tests of transitional boundary layers with 
laminar, reversing and turbulent flow increments were performed to 
decouple the effects of the parameters and to determine the effects 
during mutual interaction. 


3 



In light of -further development of calculation procedures, a high 
degree of general expressiveness was required of the tests. Although 
blade configuration tests are the real purpose of this dissertation, 
their results are always influenced by individual blade geometry. 
Therefore, a genuine reproduction of basic effects cannot be 
achieved. Because of this, I did not select this type of flow but 
the more general, flat plate boundary layer which is not influenced 
by geometry. In addition, numerous test procedures which were 
developed with the this type of boundary layer provided the advantage 
that, in addition to heat transfer measurements, more complex 
boundary layer measurements could be accomplished. 

The above requirement for more exact ways to calculate heat transfer 
at the surface of gas turbine blades has resulted more and more in 
the use of so-called field methods whereby the differential equations 
are solved numerically in connection with a turbulence model [4-91. 
Only this method was recognized as having the flexibility necessary 
to describe the complex blade boundary layer flow conditions with 
sufficient accuracy. 

In contrast, integral procedures - despite a certain amount of 
success (CIO-133, see C31, C141) - did not seem suitable for the 
calculation of heat transfer on cooled gas turbine blades. This is 
supported by the fact that, to date, no method has been developed to 
describe the long transfer zones, which are characteristic for blade 
conditions, with integral procedures. In addition, the calculation 
methods for heat transfer, which are decisive in blade cooling, and 
their main analogy equations were only partially developed and 
encouraging proposals for improvement did not evolve. The obvious 
disadvantages of integral procedures are mainly attributable to the 
fact that, with the introduction of field methods, the development of 
integral calculation methods was neglected, although their use is 
still relatively common in industry. 
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With integral procedures, the basic equations, which describe the 
boundary layer, are integrated in a normal manner to the main -flow 
direction by adopting known profile equations. This results in a 
system of common differential equations which can be solved in a 
quick and easy manner. This method is widely used because of its 
economy and simple application. In order to utilize these advantages 
in the future it seemed practical to use integral procedures and 
examine them with respect to their actual efficiency. In order to 
accomplish this, the development of a universal integral procedure 
was introduced which, as an alternative to field methods, 
demonstrates the usefulness of these methods for calculating heat 
transfer coefficients for gas turbine blades and which, at the same 
time, maintains the advantages of economy and simple application. My 
immediate goal was to provide a calculation model based on the 
integral method which calculates the influences of free stream 
turbulence, continuous reversal processes and heat transfer 
coefficients in cooled boundary layers with high pressure gradients 
without using analogy equations or purely empiric additive functions. 

1.2 Focal Points 


This dissertation contains theoretical and experimental aspects. 
Experimental work consists of heat transfer and boundary layer tests 
on the flat transitional plate boundary layer which separates the 
influences of free stream turbulence, wall cooling and negative 
pressure gradients as well as describing them during mutual 
superimposition. The experimental portion is the main component of 
this dissertation. It is the basis for the chapter "Present State of 
Research" and the reference list. 

Chapter 6 discusses the above integral procedure and makes a 
comparison with tests of other authors and those made by me. 
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PRESENT STATE OF RESEARCH 


2.1 Appearance of Transitional Boundary Lavers 

Before discussing the analysis of the influenceing parameters, I 
consider it practical to describe the basic components and properties 
of transitional boundary layers. This will explain the basic 
conditions and the surroundings which causes the interaction of 
boundary layer and influencing parameters. 

2.1.1 Laminar Boundary Lavers 

A transitional boundary layer, which is formed along a fixed surface, 
begins in a laminar manner and becomes turbulent in its further 
development. The laminar-turbulent transition occurs within a 
defined intermediate area, which is referred to as the reversal area, 
artd in which the boundary layer undergoes a continous change of 
significant consequences. 

The boundary layer, whose beginning is laminar, describes a shear 
flow in which turbulence flows are steamed out. This flow is 
controlled exclusively by viscuous impulse and energy exchange 
phenomena. This causes typical boundary layer profiles which, in the 
vicinity of the wall, distinguish themselves by large areas of linear 
velocity and temperature distribution. 


Because there is no exchange process, the calculation of laminar 
boundary layers is relatively simple. With certain conditions for 
the distribution of the free stream velocity u (x), analytically 
exact solutions of the boundary layer differential equations are 
possible. This applies especially to the flat un-accelerated plate 
boundary layer and, generally, to the wedge flow with u (x) = const x 

W 

x . These cases have "similar solutions" with the characteristic 
property that boundary layers are not dependent on the running length 
in the similarity description: 



resp. jj- « f(n = 

OD 



( 2 . 1 ) 
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In boundary layers with arbitrary velocity distribution, which 
includes the -flow around gas turbine blades in the vicinity of the 
wall, similar approximations cannot be used. Solutions to describe 
laminar heat treansfer and tangential stress conditions can only be 
obtained via the numeric integration of the boundary layer 
differential equations. 

2.1.2, The Laminar-Turbulent Boundary Laver Reversal 

After a long running time, laminar flows with total pressure 
gradients become unstable even with very small disturbances. This 
generally occurs after a characteristic value of the Reynolds number 
Re 0 , which is formed with the impulse loss thickness, is exceeded. 
The laminar boundary layer remains stable below this "neutral 
Reynolds number" (Re 0 ~ 160). If the boundary layer has higher 
impulse loss thickness Reynolds numbers, flows are generated which 
cause turbulance further along in the flow, thus generating the 
laminar-turbulent reversal process. Arnal C133 demonstrates this 
process with comprehensive hot wire tests. He shows how, with 
increasing running length, low long-wave flows are transferred into 
high-frequency turbulence which results in pronounced, turbulent 
shear stresses and profile deformations. The key phenomenon of the 
reversal is a process of generation, excitation and propagation of 
eddy systems which occurs in several phases (see C33): 

• 

1. In the first phase - after the neutral Reynold number has been 
exceeded - unstable interference waves are generated and excited in 
the boundary layer. The interference waves are two-dimensional and 
bear the names of their discoverers, i.e. Tollmien-Schlichting 
waves. They can be generated by irregularities in the external flow 
(turbulence, accoustic vibrations, etc.) and by irregularities of the 
surface. Tollmien-Schlichting waves move with approximately 307. of 
the outside velocity in the direction of the main stream; their wave 
lengths can have a thickness of several boundary layers. 
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2. Further along the path, slight irregularities in the boundary 
layer or in the -free stream may cause three-dimensional waves and the 
•formation of a three-dimensional wave system which alternates with 
the main stream ('"peak-valley development"). This results in initial 
eddy systems whose axes are in the direction of the flow and which 
rotate in pairs in an opposed manner. In a diagonal direction, this 
secondary flow generates zones with varying interference wave 
intensity thereby causing a slight three-dimensionality of the basic 
flow. The planes of the maximum interference wave amplitudes are 
called "peak profiles", the planes of minimum amplitudes are referred 
to as "valley profiles". 

3. Further down-stream, the formation of a high shear layer with a 
high eddy concentration occurs in the peak/valley profiles. In the 
breakdown process, the longitudinal eddies disintegrate at this 
location, thus generating higher-frequency eddies with the center of 
rotation in the diagonal direction (hairpin eddies). 

4. These hairpin eddies form 'wedge-shaped turbulence spots which, 
during their up-stream movement, move in and diagonal to the flow 
direction. In experiments, the occurance of turbulence spots 
resulted in boundary layer intermittance along with a chronological 
alternation of laminar and turbulent form of flow. 

5. When the turbulence spots grow together, the boundary layer 
become increasingly turbulent. 

The above processesare, for the most part, extracted from reports by 
Schubauer and Skramstadt C163, Schubauer and Klebanoff C173 and 
Klebanoff et al 1183. Overviews are provided by Tani C193, Morkovin 
C203, Reshotko C213, Schlichting C223, Eppler and Fasel C233 and in 
the dissertation by Ghon C243. Reshotko C213 and Eppler & Fasel C233 
describe current research, while Scheurer C43 and Ghon C243 include 
parameter influences and the theoretical accessibility (calculation 
procedures) into their studies. 
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The process of individual reversal phases can be influenced by 
different parameters (e.g. free stream turbulence, surface roughness, 
surface curvature, pressure gradients, etc.). Under certain 
conditions there may even be situations where one or more reversal 
phases are bypassed. An important example is the presence of high 
external turbulence (Tu > 17.). When such conditions exist, the 
presence of Tollmien-Schlichting waves was, at times, not observed 
(Morkovin C20D). 

In the above, the process of boundary layer turbulence generation was 
emphasized without an evaluation of the resulting deformation of the 
boundary profiles, or the alteration (increase) of the heat transfer, 
and wall tangential stress. Available tests (e.g. C15, 16, 173) show 
that there is a distinct difference between the turbulence process 
and a noticeable profile deformation. The boundary layer profiles, 
which are determined in a chronological manner, and their friction 
and heat transfer coefficients begin shedding laminar behavior and 
their laminar values only after the turbulent spots are formed. 
According to Schubauer and Skramstad (see C25, 263), a minimum 
intermittance of approximately 167 is necessary during turbulent 
external flow, before a deformation or the velocity profiles is 
noted. The external occurence of the reversal is thus concentrated 
on the last phase (phase 5) of the reversal process. This is 
important for the evaluation and assignment of quantitative 
statements for the beginning of the reversal, if reversal Reynolds 
coefficients were not determined by the turbulence but with the aid 
of measured velocity profile deformations or heat transfer 
distributions. 
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2.1.3 Turbulent Boundary Lavers 


In turbulent boundary layers, -flow types occur which are mainly 
generated by a regularly recurring ’burst’ process. The burst 
process occurs in a given pattern. Slower moving -flow sections are 
torn out of the zone near the wall and, during down-stream movement, 
are pushed toward the edge of the boundary layer while creating 
turbulence (see C27-29]). The turbulence of this rotational flow is 
similar to a conglomerate of eddy structures with varying sizes, the 
largest of which has the same thickness as the boundary layer. The 
smallest eddies are dissipated by the flow and are transferred into 
heat. In an energy cascade, the large eddies withdraw energy from 
the main stream; this energy is transferred in stages to smaller 
eddies and to the dissipating eddies. This turbulence mechanism 
depends very little on the viscosity. Higher viscosity will not 
result in higher dissipation. Its only effect is that the dissipation 
starts with large eddy systems (see C301). 


Turbulent boundary layers have distinct velocity, turbulence and 
tangential stress distribution structures which provide this type of 
flow with a similar character in most areas. Based on dimensional 
analyses and numerous experiments, the velocity profiles and 
turbulent shear stress profiles can be described by a universal 
interrel ati onshi p. 


T 

(y+ = -^> » 


( 2 . 2 ) 


The wall distance y is also a linear measure for the largest eddies 
with the most energy and u^, cann be assigned to them as the 
applicable measure of speed. Because of this, the standardized wall 
distance coordinate y + is compared to a turbulence Reynolds number 
which represents the energy carrying eddy systems of the turbulent 
boundary layer. 
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Illustration 2.1 is a schematic of the velocity distribution 
described with equation (2.2). The behavior o-f this velocity pro-file 
extends through several phases which are marked by characteristic 
di-f-ferences in the relative in-fluence of viscous and turbulent 
tangential forces. It is practical to break down the entire profile 
into an internal and external area. The internal area is comprised 
of the areas near the wall (y/6 < 0. 1-0.2), consisting of the laminar 
bottom layer, the intermediate layer and the logarithmic area. This 
boundary layer section distinguishes itself by almost constant 
tangential stress r = r . In the external section of the boundary 
layer, the transition from the highly turbulent flow to the almost 
eddy-free external flow occurs near the wall. 



1 - wake area; 2 - external boundary layer section; 3 - area of the 

log. wall theory; 4 - log. wall theory; 5 - wake portion; 6 - 
internal boundary layer section; 7 - viscous area; 8 - transition 
area; 9 - laminar bottom layer 

« 

Illustration 2.1; Approximation of a Turbulent Velocity Profile 
(according to Bradshaw C301) 


ORIGmr 

OF POO ft 



li 



Depending on the Reynolds number, the laminar bottom layer y + < 5 is 
extraordinai ly thin. It is dominated by viscous -forces, although in 
the chronological disintegration of the flow, a structure occurs 
which is already disturbed by continuous small scale velocity 
fluctuations and fluid elements, which penetrate this layer from the 

* 4 * 

outside. In the buffer layer (5 < y < 30), a continous transition 
occurs from the laminar bottom layer to the fully turbulent flow form 
which, in the subsequent layer, is characterized by the validity of 
the logarithmic wall theory. In the buffer area, the turbulent shear 
stress portion increases by the same amount as the decrease of the 
viscous tangential stress portion. 

Only the turbulent impulse exchange events and tangential stress 
portions are important in the subsequent area (y + > 30). The 
logarithmic wall theory applies: 


The condition t = - P u V * t w * const, is the theoretical basis of this 
relationship. However, according to test results, turbulent boundary 
layer flows are possible (e.g. when free stream turbulence is 
present) where, because of special turbulence structure, the 
logarithmic wall theory satisf actor i 1 y describes the velocity path, 
even for large wall distances with a clearly reduced tangential 
stress (Bradshaw [30]). 

In the external section of the boundary layer (Illustration 2.1), the 
velocity profile is dominated by the statutes of free or wake 
turbulence. Because of the intermittent ejection of turbulence to 
the edge of the boundary layer, turbulent boundary layer flows 
alternate with rotation-f ree shear flows with viscous properties. 
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The velocity in the external portion of the boundary layer is 
described by the wake funcion w (y/6). This function is assigned an 
additional weighting factor p(x)/* and j S then added to the 
logarithmic wall theory (see Illustration 2.1). There are a number 
of empirical approximations for this wake function (see [22, 353). 
Coles, for examples, proposed: 

w(y/«) = 2.sin2(£ . f) j (2 .4) 

In the variable wake force P(x), consideration is given to the fact 
that turbulent boundary layers often depend on their history. The 
longer life of large eddies is the reason for this memory effect. It 
requires that stored information be carried down-stream; changes of 
external conditions like pressure gradients, roughness, etc. are only 
realized at the boundary layer with a certain amount of delay 
(relaxation). Because relaxation effects decrease from the outer 
edge to the wall, the wake area is mainly affected by this. The 
logarithmic wall theory and the flow along the wall remain unchanged 
during semi -logarithmic application (as in Illustration 2.1). After 
the relaxation effects have subsided, the boundary layer is in a 
total equilibrium condition (’equilibrium boundary layer’). This 
condition, like other similar solutions in laminar cases, 
distinguishes itself by totally similar velocity profiles with 
unchanging wake areas and similar turbulence structures. Therefore, 
the convective terms (see C30D) become negligable in the maintenance 
equation which describes turbulent energy exchange events. 

The turbulent boundary layer with low impulse loss thickness reynolds 
numbers (Re,, < 5,000), which is character istic for blade flows, also 
belongs to relaxation flows. In transitional boundary layers it can 
always be found immediately after the laminar — turbulent boundary 
layer reversal. 
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Turbulent boundary layers^it^smal9^impu?se loss thickness Reynolds 
numbers are highly affected by viscosity influences. They are 
present near the wall and are caused by thickened laminar bottom 
layers; they can also be found at the outer edge of the boundary 
layer where they are caused by decreased wake portions. According to 
White C31], additional changes can be observed in the constants of 
the logarithmic wall theory (C and K) when impulse loss thickness 
Reynolds numbers Re„ are < 600. Because of this, events in turbulent 
boundary layers with low Reynolds numbers have a higher Reynolds 
number dependency on the characteristic boundary layer quantities 
like friction coefficient c, and profile shape paramater H.^. In 
addition there is a Reynolds number dependency on the production, 
dissipation, convection and diffusion of turbulent kinetic energy. 
Murlis C33] and Murlis et al [32] quantify these basic phenomena in 
comprehensive analyses with the aid of measured velocity profiles and 
turbulence distribution at the boundary layer. 


Photographs taken by Head et al [34] (visualisation with smoke) 
( illustration 2.2 ) show the differences of the global turbulence 
structure for impulse loss thickness Reynolds numbers of 500 to 
10,000. With lower Reynolds numbers (Illustration 2.2a>, vortex 
loops with highly viscous portions occur at the edge of the boundary 
layer. When the Reynolds number is increased, they are replaced by 
hairpin eddies (Illustration 2.2b), thus making the transition to the 
equilibrium boundary layer. 

■4 ■ StrSmungsrlchtung 



(a) Re 2 ■ 5000 (*>) Re2 = 10000 

1 - flow direction 

Illustrations 2.2a and b; Turbulence Structures in the Boundary 
Layer with Re„ = 5,000 (a) and Re^ = 10,000 (b) (according to Head 

[32] 
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2.2 Overview of Available Literature 


2.2.1 Boundary Lavers under the Influence of Free Flow Turbulence 

The influence of free stream turbulence extends to all areas of a 
transitional boundary layer. In addition to the influence of the 
reversal process, effects occur mainly in turbulent boundary layers. 
In laminar boundary layer zones, events caused by the degree of 
turbulence occur predominantly when pressure gradients are present. 

2.2. 1.1 Laminar Boundary Lavers and Reversal with Varying Free 
Stream Turbulence 

The generation of heat transfer on cooled gas turbine blades is 
determined by the laminar-turbulent boundary layer reversal. 
Beginning and length of the reversal process are subjected to a 
multitude of influencing parameters, with the most important being 
free stream turbulence. Several publications demonstrate their 
significance under flat plate conditions with a mostly adiabatic flow 
[35-463; this was proven with comprehensive quantitative information 
C37-39, 42, 43, 443. Because pressure gradients also change the 
reversal process, many studies determined both parameters together 
[35, 37, 38, 41-453. All authors recognized that, with increasing 
turbulence, an upstream shift of the reversal and a shortening of the 
reversal length occurred. In flows with pressure gradients it was 
determined that with a super imposition of both parameters the prefix 
of the pressure gradient has a decisive significance: positive 
pressure gradients excite turbulence movements in the boundary layer 
and cause, similar to free stream turbulence, an early reversal. In 
accelerated flows, a stabilization of the laminar boundary layer and 
a delay of the reversal start taking place. Negative pressure 
gradients thus counter-act the external turbulence, and the relative 
strength of both components becomes a decisive factor for the final 
formation of the reversal . 
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Experimental test results are based on observations of the heat 
transfer coefficient path or on scanning of the boundary layer with 
hot wire and Pitot tubes. Hall and Gibbings C433 (1972) provide a 
comprehensive summary and evaluation of all experiments made up to 
that time and of models - developed on a semi -empirical basis - for 
the prediction of the impulse loss thickness Reynolds number at the 
initiation of the reversal. In an analysis of their tests and 
previous studies they express the opinion that, with increasing 
turbulence, the reversal Reynolds number approaches a lower threshold 
which is considered a stability limit and which must not be exceeded 
even by the highest degree of turbulence. The Reynolds number at the 
neutral point (Re^ ~ 160), which is used in the stability theory, is 
named as the limit. This theory is supported in recent studies by 
Abu-Ghannam and Shaw C453. As of this date, however, a clear 
explanation of this fact which is based on experimental data, has not 
been possible because of instabilities and the spread of measured 
reversal Reynolds numbers (see C453). In addition, little data is 
available which deal with the high turbulence area in a satisfactory 
manner . 

The significance of pressure gradients in flows with external 
turbulence was expressed by correlations in several publications. 
Abu-Ghannam and Shaw [453 provide an overview of these 
relationships. The impulse loss thickness Reynolds number is 
predominant for the global characterization of the reversal 
initiation. 

In accordance with experimental findings, the boundary layer 
stabilization, which occurs with accelerated flows, results in an 
additional down-stream relocation of the reversal initiation and in 
an increase of the reversal length. With lower degrees of 
turbulence, this event is coupled with an increase of the impulse 
loss thickness Reynolds number at the beginning of the reversal . In 
accordance with recent studies by Abu-Ghannam and Shaw [453 and Blair 
[393, the latter event no longer applies to accelerated boundary 
layers with more than 17. of free stream turbulence: 
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In contrast to earlier explanations by Hall and Gibbings [433, with 
turbulence in excess of IV. the same impulse loss thickness Reynolds 
number is present at the beginning of the reversal as with identical 
external turbulence in the case without acceleration. Therefore, the 
smaller boundary layer or impulse loss thickness increase is the sole 
reason for the down-stream relocation of the reversal. 

Abu-Ghannam and Shaw [45] deal with the significance of the boundary 
layer’s history. They describe the difficulties of a correlation, 
which is based on reversal prediction, if the reversal occurs in 
flows with non-constant acceleration. 

An estimation of the reversal length has only been made in few 
studies (see C 45 3 ) and only for boundary layers without pressure 
gradients. So far, measurements in flows with applied pressure 
gradients have not shown any noticeable correlation possibilities for 
reversal length. 

Buyuktur et al [413 and Hi si op (see E433) pointed out the 
significance of free stream turbulence structure as an additional 
parameter which influences the reversal process. In experiments, 
various turbulence structures with the same degree of turbulence were 
generated by various types of turbulence nozzles. While Buyuktur et 
al C4 1 3 did not realize any effects caused by changed linear measures 
of the free stream turbulence, Hi si op (see 1433) recognized these 
tendencies. Hall and Gibbings [433 believe that the great 
differences and spreads of the reversal Reynolds numbers, which were 
measured by various authors, are caused by linear measurement 
influences. 
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The laminar boundary behavior under the influence of free stream 
turbulence was mostly evaluated with the aid of measured heat 
transfer distributions C 35-41, 441. It was mutually determined that 
laminar boundary layers without pressure gradients are not subject to 
any influences, but that laminar boundary layers in an .acceleratedl 
flow experience a slight increase of the heat transfer. This effect 
in accelerated flows was verified by Junkhan and Serovy C351 in 
boundary layer measurements. In these experiments, fuller velocity 
profiles were observed under the influence of external turbulence; 
however, they continue to show the typical laminar boundary layer 
shape. 

2.2. 1.2 Turbulent Boundary Lavers with Free Stream Turbulence 

The influence of free stream turbulence on turbulent boundary layers 
has been the subject of many experimental, analytical and numeric 
studies Csee 35-41, 46-631. With adiabatic boundary layer flows, all 
authors agreed that free stream turbulence results in an increase of 
the friction coefficient c^ and in a higher boundary layer growth. 
Higher tangential stress coefficients and boundary layer growth are 
in agreement with fuller boundary layer profiles which can generally 

be identified by greater displacement thicknesses 6^, impulse loss 

thicknesses and, especially, by the decreased form parameters 

However, despite a restructuring of the boundary layer 
profiles, the validity of the logarithmic wall theory has proven to 
be unchanged so that the deformation of the profiles with semi- 

logarithmic application (see Illustration 2.1) extends exclusively to 
the external area of the boundary layer (wake portion). The wake 

portion decreases with increasing turbulence. 
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Kline et al C463 and Junkhan et al [353 were the -first to demonstrate 
such effects. Their theory was mostly based on the analysis of 
velocity profiles. With additional registration of turbulence 
structures, Charney et al [473, Huffman et al [483 and Evans C501 
were able to show an increased turbulence level, stronger turbulent 
exchange events and increased turbulent shear stresses within the 
boundary layers as a result of free stream turbulence. Huffman C481 
balanced the turbulent energy exchange within the boundary layer. By 
doing so, he was able to demonstrate that production and dissipation 
of turbulent kinetic energy are increased and, compared to the 
convective portions, dominate in the same manner as in boundary 
layers without external turbulence. In addition, he provides a 
comprehensive analysis about the change of the turbulent mixing path 
and the turbulent dissipation length as a result of free stream 
turbulence. Evans C503 describes the increase of individual 
turbulant fluctuation components in and vertical to the main flow 
direction caused by external turbulence. In his tests he further 
determines that an increase of the impulse loss thickness does not 
necessarily have to occur with increased free stream turbulence. The 
deformation of the profiles near the wall and at the edge of the 
boundary layer can adopt shapes which compensate each other in the 
impulse loss thickness determination. These effects were also 
discovered by Raghunathan et al [603 in more recent studies. 

Green C513 and Robertson S< Holt [523 also perform analytical and 
theoretical studies of boundary layer phenomena during free stream 
turbulence. They were taken up by Bradshaw [533 who summarises them 
in an overview which includes additional experimental results. It 
provides a theoretical reason for the linear increase of the 
tangential stress coefficient with increasing external turbulence 
which was demonstrated by numerous experimental studies (e.g. C47, 
49, 52, 543). Experiments have shown that, compared to bounde'-y 
layers with identical impulse loss thickness Reynolds number, a 2-4 
percent increase of the tangential stress must be expected for each 
percentile of external turbulence. The reasons for this increase of 
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are, 


Cj. are, in addition to technical diff iculties with measuring 
equipment, instabilities which are a result of additional Reynolds 
number effects and varying free stream turbulence structures [49, 
533. 


The significance of turbulence structure was realized in recent 
studies (Simonich and Bradshaw [543, Bradshaw [543, Hancock [563, 
Hancock et al [593, Meier and Kreplin [5B3, Blair [36, 393. In order 
to characterize the structure, a turbulent linear measure was created 
with the attenuation, behaviour or with the energy spectrum of the 
local free stream turbulence. This measure is also referred to as 
dissipation length. The studies show that the strongest influence of 
external influence occurs in those cases where these linear measures 
have almost the same magnitude as the boundary layer thickness. 
According to Hancock et al [593, in turbulent eddy structures with 
large linear measures turbulence motions at the edge of the boundary 
layer are impaired by the vicinity of the wall. This causes their 
interaction with the boundary layer to attenuate . With these 
considerations, Bradshaw and Hancock [55, 563 also studied the 
measurements by Meier and Kreplin [5B3 where, for low degrees of 
turbulence below 17., only very minor increases of the tangential 
stress coefficient c^ were noted. There was, however, a square-law 
dependency on the turbulence intensity. At the same time, linear 
measure influences were cited as possible reasons for findings in 
recent studies [52, 573 (see {53, 55, 593) whereby, for high 
turbulence intensities above 87., the increase of the tangential 
stress coefficient c^. asymptotically strives toward a limit of 
approximately 20-307.. 

In another study, which provides a detailed introduction of boundary 
layer measurements and heat transfer measurements, Blair ([36, 383) 
discusses and correlates the mutual importance of linear measure and 
Reynolds number influences. He states that Reynolds number 
influences have the same effect as linear measures, which are too 
large, and thus contribute to a decrease of the turbulence intensity 
influence. 
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Raghunathan C583 expands the studies about the influence of -free 
stream turbulence in turbulent boundary layers to flows with high 
velocities (ma = 0.8). While circumventing the above Reynolds number 
and linear measure influences he quantitatively demonstrates that, 
even under these conditions, the same turbulence intensity effects 
exist as in previous studies with low free stream velocity. 

Numerous experiments with the adiabatic boundary layer provide 
consistent information about the behavior of the flow boundary 
layer. However, distinct differences of opinion are noted in the /17 

evaluation of the turbulent heat transfer under the influence of free 
stream turbulence. In the publications of Edwards and Furber [613, 

Reynolds et al [633, Kestin and Maeder [413, Biiyiiktiir and Kestin 
[403, Junkhan and Serovy [353 and Cosigny et al [443 it is stated 
that the influence of the turbulence on the turbulent heat transfer 
is negligible. Tests by Sugawara and Sato [623 and recent studies by 
Blair [36, 393, Simonich and Bradshaw [543 and F'edisius et al [573 
show that increased rates of heat transfer were observed which are 
comparable with increases of the tangential stress. The discrepancy 
in this evaluation can be explained by the fact that, in many cases, 
the heat transfer coefficients were compared at locations with 
identical running length or identical running length Reynolds 
numbers. This procedure does not consider that boundary layers with 
increased external turbulence increase at a higher rate. Therefore, 
thicker boundary layers with increased heat transfer resistance occur 
at the same location x or Re w . This process counteracts the actual 
increase of the heat transfer; consequently, no net effect may be 
visible when the heat transfer coefficients are locally applied. The 
heat transfer changes caused by free stream turbulence will only 
occur when boundary layers with identical local characteristics (e.g. 
impulse loss thickness) are used for the evaluation. In 
comprehensi ve comparisons of boundary layers with identical impulse 
loss thickness Reynolds numbers, Simonich et al [543, Blair [36, 393 
and F'edisius [573 demonstrate that there is a close interrelation 
between the behavior of the heat 
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transfer and the wall tangential stress. The increase of the heat 
transfer with increasing turbulence as well as the significance of 
linear measures and Reynolds number effects, which we know from the 
tangential stress coefficients, were thus discovered. 

Detailed boundary layer measurements which show the reaction of 
velocity and temperature profiles to increased free stream turbulence 
are described by Simonich et al [54], Blair C36-393 and Pedisius et 
al C573. All authors demonstrate that, even in the temperature 
profile, the influence of the turbulence is mainly in the wake area. 

2.2.2 Boundary Lavers under the Influence of Varying 
Material Values 

Boundary layers with heat transfer and constant material properties, 
which were used as a basis for the above analyses of the influence of 
free stream turbulence on the heat transfer, are an exception in 
practical engineering. They are often generated in tests to 
demonstrate basic heat transfer phenomena. In contrast, boundary 
layer flows, which are generated around gas turbine blades with wall 
cooling, contain high temperature gradients which result in 
pronounced material value variations within the boundary layer. This 
has a major effect on air and comparable combustible gases with 
respect to density p, viscosity P, and heat conductivity^. The 
specific heat c^ and the Prandtl number Pr experience only minor 
alterations. 

2. 2. 2.1 Laminar Boundary Lavers and Reversal during Wall Cooling 

Statements about the influence of variable material values for 
laminar boundary layers are almost exclusively based on theoretical 
or numerical boundary layer analyses C64-673. An examination of 
available literature revealed no representative or conclusive test 
results which provide a detailed description of the laminar heat 
transfer and the laminar wall tangential stress or which demonstrate 
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conclusive boundary layer pro-file measurements- According to 
theoretical calculations, only minor increases of heat transfer and 
wall tangential stress occur with wall cooling, despite a distinct 
deformation of the velocity and temperature profiles. The reason for 
these insignificant effects is that the velocity and temperature 
gradients near the wall, which have been increased by the cooling, 
are almost totally offset by a simultaneous heat conductivity and 
viscosity decrease. 

Eckert (see [64, 681) developed the concept of considering the 
influence of material values in such a manner that they are 
established for a suitable, empirically determined reference 
temperature; they are then used for the relations obtained in the 
case of constant fluid properties. Eckert proposes the arithmetic 
average between free stream temperature and wall temperature as the 
reference temperature for incompressible laminar boundary layers. 

A second possibility of registering material value influences is to 
expand correlations determined for constant fluid properties by 
temperature correction functions. Proposals on this subject are made 
by Kays [643 and Brand & Wehle [653. 

Schlichting C223 discusses the influence of material value changes in 
the boundary layer on the laminar-turbulent boundary layer reversal 
for flow cases without external turbulence. Based on the linear 
stability theory, a profile theory referred to as reversing point 
criteria, calculations by Cebeci and Smith (see C223), and 
measurements by Liepmann (see [223 ) the opinion is expressed that 
wall cooling increases the stability of the laminar boundary layer, 
thus delaying the reversal of the boundary layer. Shon [223 
summarizes this theory in a general study on the subject of laminar- 
turbulent boundary layer reversal and provides additional background 
information. 
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More recent experiments with supersonic flows made by Watson E693 and 

Fischer E703, however, contradict the idea of boundary layer 

stabilization by wall cooling. According to their measurement 

results, no altered reversal occurs in flows with wall cooling if the 

same free stream Reynolds number is maintained in a varying 

temperature ratio T /T . 

w « 

2. 2. 2.2 Turbulent Boundary Lavers with Varying Wall and 
Free Stream Temperatures 

The influence of changeable material values originally became a 
dominant part of boundary layer research with the study of 
compressible boundary layers. The friction heating, which develops 
within these boundary layers, causes temperature increases which - 
even without heat transfer at the wall - are related to fluid 
property changes of such magnitudes that they can no longer be 
neglected. In addition to experimental gathering of information 
about specific events, it was the goal of studies performed on this 
subject to find those conditions with which the boundary layer 
differential equations could be transformed to a simpler 
incompressible form. Another focal point was to convert these 
transformations to turbulent boundary layer profiles, friction 
coefficients and, at a later time, to heat transfer coefficients. 
The relationships and correlations obtained for constant fluid 
properties would then be used for further research (see E713). In 
studies about the influence of compressibility, Spalding and Chi C72] 
included the aspects of additional heat transfer; Gran et al E733, 
Watson E693 and Lee et al E743 included the specific aspects of wall 
cooling. In contrast to the influence of compressibility, wall 
cooling affects a slight increase of heat transfer and wall 
tangential stress. Together with other research about compressible 
turbulent boundary layers, White C 70 3 summarized these effects in an 
overview. 
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Parallel to research with compressible boundary layers, interest 
increased on the significance of variable material values in 
incompressible boundary layers in subsonic flows. This has a special 
significance for the flow conditions in gas turbines. Because the 
exhaust Mach numbers of current turbine nozzles are only in the area 
of 

0.6 < Ma < 1.2 

and local values in excess of Ma = 1.8 rarely occur, no pronounced 
compressibility influences are expected (see [43). 

Tests were made for boundary layers with wall cooling [75-803 as well 
as for flows with high wall heating [81-843. Comparisons with 

compressible boundary layers were frequently made and the 

relationships and transformation equations obtained for them were 
assumed and further developed. For example, Spalding and Chi C753 
used equations, which were obtained with the description of 

compressible boundary layers, for the analytical computation of heat 
transfer and wall tangential stress and expanded them for the 
calculation of incompressible boundary layers. I performed the test 
for incompressible flows in the preparation of this dissertation with 
highly cooled flat plates. With increased wall cooling, these 
measurements show an increase of the heat transfer which, with a 
running length Reynolds number of Re^ = 10^* and a decrease of the 

temperature ratio to T T„ = 0.4, is between 15 and 20 7.. Similar 
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local increase rates were observed by Back and colleagues C673. 

The effects are more distinct if, instead of wall cooling, the 
surface is heated. Instead of a heat transfer increase, there is a 
decrease in the local heat transfer coefficients which, according to 
the measurement results achieved by Rotta C833 and the analyses by 
Kays [643 can be approximated with the factor 


In further studies of the influence of wall cooling, Back, Cuffel and 
llassier C763 and Sill C783 made detailled analyses of turbulent 
boundary layer profiles and friction coefficients. The illustration 
of velocity and temperature profiles in a semi-logarithmic plot 
demonstrates differences to the isothermal flow; these differences 
are especially pronounced with respect to the logarithmic wall 
theory. To describe these events, Rotta C863 developed expanded 
analytical equations for the logarithmic wall theories of flow and 
temperature boundary layers. 

The development of optical measuring methods (Laser-Bi -Focus and 
Laser-Doppler measuring procedure) provided the opportunity for 
studying the transport events within the boundary layer, which are 
responsible for profile deformations and tangential stress 
alterations, even under high cooling or wall heat conditions. With 
optical measurement procedures Eriksen, Wittig et al [79, 801 
performed tests in highly cooled boundary layers, while Cheng and Ng 
C841 made their tests in highly heated subsonic boundary layers. 
Independent of the type of heat transfer (heat, cooling), it was 
demonstrated that material value variations have no influence on the 
turbulence structure. The turbulent shear stress distribution - p uV 
followed the influence of the distinct temperature field only because 
of its immediate relationship with the density p. A great number of 
details are contained in the test results by Eriksen et al C79, 803. 
In addition to turbulence and shear stress distribution, the authors 
determined turbulent energy balances, triple products of the 
turbulent fluctuation quantity, and linear dissipation measures in 
the boundary layer. 

2.2.3 Boundary Lavers under the Influence of Pressure Gradients 

The pressure gradient is the parameter with the most extensive 
consequences for the boundary layer of the gas turbine blade under 
flow conditions. It can result in flow separation or relaminization 
and plays a major part in the formation of the laminar-turbulent 
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boundary layer reversal. The reasons -for this are the destabilizing 
effects of positive pressure gradients and the stabilizing effects of 
negative pressure gradients (acceleration) . Back C663 demonstrates 
this with the aid of numeric calculations for accelerated laminar 
boundary layers. Stabilization occurs because the profile gradients 
near the wall increase thereby generating fuller boundary layer 
profiles which distinguish themselves by increased heat transfer and 
increased wall tangential stress, but also by decreased boundary 
layer thickness growth. 

The above effects are also present in turbulent boundary layers. In 
addition to the stabilization of the boundary layer, an attentuation 
of the boundary layer turbulence occurs which, in extreme cases, may 
result in relaminization. According to Kays and Moffat C863, 
relaminization occurs in flows without free stream turbulence when 
the local acceleration parameter k (x) exceeds the value 

v du 

k = — • -r-- h, 3-10" b 
u w 2 

For comparison, the real flow in the intial stages of gas turbine 
blades frequently has much higher acceleration intensities. There 
are, at times, k-values with a magnitude of 10 ° . 

For the most part, the analysis of pressure gradient influences is 

determined by the observation of turbulent boundary layers and 
laminar-turbulent reversal events. The interest in purely laminar 
boundary layers without free stream turbulence has, for the most 
part, been exhausted. Due to a lack of turbulence events, this type 
of flow has theoretically been revealed and can be mastered in a 
relatively simple manner with available computers. 
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Turbulent boundary layers were experimentally studied -for delayed as 
well as for accelerated flows [86-921. Clauser [931, Rotta [941 and 
Cebeci [911 published overviews on this subject. Their experiments 
were mainly concentrated on the representation of the boundary layer 
profile [86-921 and the analysis of heat transfer distributions [86, 
89, 901. In these tests it was discovered that the logarithmic wall 

theory which applies to unaccelerated boundary layers can, without 
any changes, be applied to boundary layers with pressure gradients 
even for compressible flows, if there is no danger of flow separation 
or relaminization tendencies. Therefore, the effects of pressure 
gradients occur only in the wake area and near the wall. The 
deformations near the wall can be explained by the wall bond of the 
boundary layer differential equations (see chapter 6). Events in the 
external area can be identified by the fact that the wake portion 
decreases with acceleration. The wake portion increases when a delay 
occurs. Analytical possibilities of describing these transf ormations 
were developed by Townsend [95], Szablewsky [963, Coles [973 and 
McDonald [983. 

Tests performed by Kays and Moffat [863 and Back et al (see [913) are 
particularly interesting, because they deal with basic properties of 
turbulent boundary layer development in flows with pressure gradients 
and because they demonstrate the differences between flow and 
temperature boundary layers. Measurements by Kays and Moffat [863 
show that, for constant accelerated flows with k(x) = const., the 
boundary layer strives toward a condition of equilibrium which 
distinguishes itself by constant impulse loss thickness Reynolds 
numbers, shape parameters, tangential stress coefficients and 
’similar’ velocity profile shapes. This, however, does not apply to 
the temperature boundary layer; consequently, the Stanton 
coefficients continue to fall with the running length. Additionally, 
this report [863 notes that the turbulent Stanton coefficients 
decrease with acceleration; because of this, special consideration 
must be given to the starting conditions of the boundary layer before 
entering the acceleration phase. The enthalpy loss thickness/impulse 
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loss thickness relationship at the beginning of the pressure gradient 
is named as the characteristic parameter of the boundary layer’s 
history. 

Back et al (see Cebeci C91]) performed tests in a nozzle with the 
greatest acceleration (k < 8 x 10 ^) . The report makes a direct 

comparison between the different reactions of temperature and 
velocity profiles when a staggered acceleration increase, which 
culminates in relaminization, occurs. In the velocity profile, 

relaminization occurs earlier than indicated in the distribution of 
the heat transfer coefficients or in the temperature profile. A 
similar divergence of flow and temperature boundary layers was 
discovered by Blair C37, 393 in the laminar-turbulent reversal 
behavior. In spite of a relatively low acceleration (k < 1 >; 10 ^), 
the reversal in the flow boundary layer was terminated noteably 
faster than in the heat transfer distribution along the plate. The 
reason for the sluggish behavior of the temperature boundary layer is 
that the temperature profile is affected only indirectly by the 
pressure gradient via convective events in the boundary layer, which 
depend on the deformation of the velocity profile (see boundary layer 
differential equations in chapter 6.2). 

The importance of pressure gradients for the laminar-turbulent 
reversal process was discussed in the previous chapter on the 
influence of free stream turbulence. Because of this, I will not 
deal with it a second time. Additionally, I will not provide 
information about positive pressure gradients. As explained before, 
I am mainly concerned with accelerated flows. 
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EXPERIMENTAL GOALS 


A reading of the af orementi oned literature revealed a number of reports 
which deal with the influences o-f free stream turbulence, variable 
material values and pressure gradients. These reports show the 
significant effects of these parameters and are supported by 
quantitative statements. 

Well-known experimental tests were preferred for tests of. turbulent 
boundary layers. Therefore, the number of those experiments which 
expand the knowledge of the reversal process to laminar boundary 
layer section is comparatively small. Laminar boundary layers are 
almost exclusively discussed on a theoretical basis; this applies 
especially to the influence of valiable material values on wall 
cooling. 

The measurements themselves are based on flow parameters which are 
rarely found in a turbine. For example, in the examination and 
further development of calculation procedures there are no boundary 
layer studies which determine the range of low Reynolds numbers, 
which is characteristic for blade flow. Therefore, an accurate 
description of the additional effect caused by low Reynolds numbers 
does not exist. Additionally, there is no basic data on boundary 
layers with free stream turbulence which are characterized by 
acceleration intensities similar to those in a turbine. In a 
turbine, the acceleration parameter k (see equation 3.1) achieved 
values of a 10 J magnitude. In current experiments with transitional 
boundary layers with free stream turbulence, only acceleration forces 
below k = 1 x 10 ^ (see [391) were achieved. From the high increase 
of laminar heat transfer (e.g. Krisnamoorthy C993), which was 
observed in blade tests, it was concluded that special effects occur 
when highly accelerated flows with a strong external turbulence are 
appl ied. 
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Another decisive -factor of this study was that, to date, no 
experiments have been published which describe the influence of 
external turbulence and pressure gradients under turbine-like cooling 
conditions. The reversal behavior during wall cooling was as unclear 
as the effect of the cooling on the local heat transfer rates. 
Present knowledge on the influence of free stream turbulence was 
obtained with adiabatic or slightly heated surfaces, i.e. only for 
flows with almost constant material values. Because of this, it is 
also unclear whether turbulence intensity effects in turbulent 
boundary layers experience a change due to strong wall cooling. 

Because of the above conditions, further studies were required which, 
for a direct relationship to the turbine flow, describe transitional 
boundary layers in the area of low Reynolds numbers and, when under 
the influence of high free stream turbulence, intensive wall cooling 
and high free stream accelerations occur. 

To develop experimental data, a flat boundary layer tunnel was 
constructed which allows fundamental studies of transitional boundary 
layers with external turbulence, free stream acceleration and wall 
cooling. Chracteristic flow conditions realised are: 


- free stream temperatures of 

- temperature conditions of 

- free stream turbulence 
intensities of 

- free stream acceleration 
values of 


T = 300 - 600 K. 

W 

T /T = 1.0 - 0.53 
w <# 

T = 1.6 - 9.37. and 

u -6 
k = 0 - 6.3 x 10 


The acceleration parameter k is defined as follows: 

v du 

k= ^7r <3.i) 

Of> 

The running length Reynolds numbers at the end of the plate were 
between 0.25 x 10^ and 1.0 x 10^. This suf f f iciently met the demands 
of turbine-relevant Reynolds number conditions. 
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The experiments were made in individual stages. For reasons of 
comparison, I began by studying the adiabatic boundary layer with 
varying free stream turbulence. I then superimposed wall cooling and 

4 

pressure gradients. In addition to local heat transfer coefficients 
I also measured boundary layer profiles and wall tangential stress 
values for unaccelerated flow. This was accomplished with mechanical 
flow and temperature probes. Due to technical constraints, tests 
with accelerated flow were limited to the determination of local heat 
transfer numbers. In this case, practical boundary layer 
measurements were precluded by testing limitations and boundary layer 
thicknesses which were estimated to be extremely thin. 

Free stream turbulence was determined with the aid of Laser Doppler 
anemometry. The use of this optical test method enabled a reaction- 
free determination of individual turbulence intensity components and 
the registration of the continuous turbulence intensity path over the 
measuring plate. 

4. TEST DESIGN AND EXPERIMENTftL METHODS 


4. 1 Test Unit 


For the experimental determination of heat tranfer coefficients under 
the ideal conditions of a flat plate, we developed a measuring 
section which could be incorporated in the available hot gas unit of 
the Institute for Thermal Dynamics. Illustration 4.1 shows a 
schematic of the entire test unit which consists of a compressor, an 
electric air heater, a steadying zone, the measuring section and an 
exhaust duct. 
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The air -from the compressor is -fed into the measuring section via an 
electric air heater and a steadying zone. The air mass stream, which 
can be adjusted with a bypass valve, and the adjustable -filament 
power (120 kW max.) enable hot gas temperatures o-f up to 600 K. 

A stabilizing zone is located between the air heater and the test 
area. This zone contains a mixing unit, a -flow straightener and 
several wire sieves to steady and equalize the -flow. The hot gases 
are fed into the adjacent measuring section via an inlet nozzle 
(surface contraction ratio 1:25). 

The height and, consequently, the profile of the adjoining tunnel 
outlet are adjustable which enables variation in the pressure level 
within the measuring section. An additional hot air blower will 
drain bottom and side wall boundary layers from the test plate. 



1 - suction filter; 2 - compressor; 3 - air heater; 4 - mixing unit; 

5 - sieves; 6 - straightener; 7 - steadying zone; 8 - particle probe; 

9 - measuring section; 10 - suction fan. 

Illustration 4.1: Hot Air Test Unit with Measuring Section to 

Determine Heat Transfer 
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4.2 


Measuring Section 


My primary goal was the determination of local heat transfer 
coefficients with varying free stream turbulence, free stream 
temperatures and free stream acceleration. Additional boundary layer 
measurements were made for a detailed demonstration of boundary layer 
events assigned to measure heat transfer. In addition, they 
provided an opportunity to compare initial studies of the isothermal 
boundary layer with results available in literature. 


4.2.1 Design Description 

The design of the measuring section is presented in Illustration 
4.2. The measuring section has a rectangular profile of 100 x 125 mm 
and a total length of 800 mm. The principal part of the tunnel is a 
copper test plate with a length of 440 mm and a width of 125 mm. The 
test plate is sub-divided into 28 individual parts which contain one 
cooling hole each. This allows independent cooling of each 
individual segment. At the plate surface, the segments are inter — 
connected with a thin bar (1.5 mm), which enables mutual heat 
exchange to be as low as possible. The length of the bar is 3 mm. 
In the flow directions, the segments have a length of 12 mm, a width 
of 115 mm and a height of 68 mm. With the temperature difference in 
individual segments, which is measured by two 0.5 mm sheathed 
thermocouples at a distance of 1.5 and 40 mm from the plate surface, 
we obtain - under the assumption of a one-dimensional heat conduction 
- the heat flow removed from the hot gas and, consequently, the heat 
transfer coefficient. (Required measures for the consideration of 
the influence of the connecting bars are described in following 
chapters). In our tests, the cooling of the segments was adjusted in 
such a manner that the temperature along the plate surface was 
constant. Thermal insulation mats protected the segments from heat 
from the sides (see Illustration 4.5 ). The measuring plate is 
supported by a base plate. An equalization of thermal expansion 
between the cooled plate and warmer tunnel walls is possible. 



1 - inlet nozzle; 2 - suction; 3 - turbulence grid; 4 - flow 

direction; 5 - position pickup; 6 - boundary layer probe; 7 - 

suction; 8 - cooling water supply; 9 - measuring plate. 


Illustration 4.2; Boundary Layer Tunnel with Measuring Plate 

The generation of varying -free stream turbulence intensities is 
achieved by round bar grids with varying grid diameters (see 
Illustration 4.4) . The grids may be inserted in the tunnel 150 - 200 
mm in front of the leading edge of the plate. Because free stream 
turbulence is determined with optical methods (Laser-Doppler- 
Anemometry), the side walls of the measuring tunnel were equipped 
with large windows. 

The tunnel lid has four openings for the insertion of the flow 
probes. This enabled the registration of boundary layer profiles at 
a distance of 120, 220, 320 and 420 mm from the leading edge of the 
plate. 
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In order to achieve an initial laminar boundary layer, it was 
necessary to expand the test plate with an uncooled, 15 mm long brass 
plate with a elliptic leading edge profile. In addition, the tunnel 
inlet flow on the bottom had to be drawn off at the leading edge of 
the plate to obtain a well defined boundary layer beginning. The 
boundary layer draw-off was adjusted in such a manner that the static 
pressure path over the leading edge of the plate remained almost 
constant. This guaranteed an isokinetic suction. Consequently, 
there were no flow disturbances at the leading edge of the plate 
which could have caused a premature laminar-turbulent boundary layer 
reversal on the plate. We examined a possible dependency of the 
laminar-turbulent separation reversal on the suction. An up-stream 
displacement of the reversal was only noticed when the suction was 
weak, i.e. when the pressure increased at the leading edge of the 
plate. A high degree of suction had no effect on the reversal or on 
the measured heat transfer coefficient. 

To avoid an influence on the heat transfer measurement by side wall 
boundary layers, they were drained in the immediate vicinity of the 
plate’s leading edge. This requires a replacement of the round side 
windows by appropriate suction devices. The required suction 
strength was determined by measuring the side wall boundary layer 
thicknesses and with the aid of heat transfer test measurements. 

4.2.2 Expanded Design to Generate Pressure Gradients 

The desired pressure gradients in the flow were achieved by placing 
inserts into the tunnel. The inserts were designed in accordance 
with given acceleration and retardation paths. The contour inserts 
used were designed in such a manner that the quality of the pressure 
gradient paths is analogous to the blade flow on the suction and on 
the pressure side. The two contours are shown in Illustration 4.3 . 
The contour on the "pressure side" accelerates the flow in an even 
manner almost throughout the entire length of the plate. With the 
"suction side" contour, the acceleration becomes, after 170 mm, a 
light, separation-free retardation. 
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1 - contour, 

plate. 



Illustration 4,5; Pressure Gradient Contours 

a) Contour 2 ('suction side’) 

■for high acceleration <k(x) = const.) 

b) Contour 1 (’pressure side’) 

for low acceleration (k(x) = const.) 

The geometric design of the ’suction side’ contour is a compromise 
between the technical solution of the heat transfer in the 

acceleration range (with an established segment width S) and the 
achieved maximum acceleration. When designing the ’pressure side’ 
contour, I intended an even (k(x) = const.) acceleration of the flow 
on those velocities, which were expected at the plate end on the 
’suction side’ contour. With constant acceleration, acceleration 

values k(x) achieved with the contours were between k * 1.2 x 10 ^ 
and k = 2.5 x 10 6 . Depending on temperature and flow velocity, 

maximum values of up to 6.3 x 10 ^ were obtained on the ’suction 
side’ contour. 


For clarif ication of these fundamental velocity and k-parameter paths 
I would like to refer to illustration 5.45 in chapter 5.6. 
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Additions to Generate Turbulence 


The measurement unit was designed in such a manner that, with t H e 
installation of various tubulence grids, a change of the free stream 
turbulence could be made. In order to correctly simulate conditions 
in a gas turbine, top turbulence intensities of 8 - 10% and higher 
were strived for. Known tests by Rose Cl 00 3 and Baines 8< Peterson 
[1013 (see Wittig, Eriksen, et al [1103) were used for the grid 
design. A turbulence intensity range of approximately 2 - 11 % was 
achieved with the aid of 4 grids. In order to generate isotropic 
turbulence, the grids were manufactured from smooth, round wires with 
square meshes and an obstruction ratio of 0.64. The grids for the 
boundary layer tunnel are depicted in Illustration 4.4 . The 
geometric data of these grids and an approximation of achieved 
turbulence intensities are contained in Table 4. 1 . That portion of 
the surface, which is occupied by the grid bars, is referred to as 
the degrees of obstruction. 


Grid 

Mesh Width 

Wire Strength 

Degree of 

Degree of 


Cmm] 

[mm] 

Obstruction 

Turbulence [%3 

1 

5 

1 

0.64 

2-3 

2 

10 

2 

0.64 

3-5 

•-« 

20 

4 

0.64 

5-7 

4 

24 

8 

o 

m 

4 * 

8-11 


*) With given tunnel dimensions a degree of obstruction of 0.64, with 
a coarse grid 4, proved unfavorable. 

Table 4.1: Specifications of the Turbulence Grids Used 
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Illustration 4.4; Round Bar Grids -for Turbulence Generation 

(Grids 1, 2, 3 and 4) 


4. 2. 4 Cooling the Measuring Plate 


The measuring plate was cooled with the thermostatic water-cooling 
cycle, which provided a chronologically constant coolant temperature 
at the inlets of individual segments. 


B] 




The coaling design required an estimation of necessary cooling mass 
flows and coolant temperatures. The ranges of acceptable mass flows 
and coolant temperatures were determined in such a way that - on the 
one hand - excessive mechanical stress due to excessive charges and - 
on the other hand - high coolant heating, which could disturb the one- 
dimensonal heat conductivity process, was avoided. Therefore, the 
influence of the coolant heating was studied in the segment with the 
aid of two-dimensional temperature field calculations (’finite 
elements’) and by simulating realistic cooling and gas parameters. 
The calculations showed that no influence on the measured result is 
to be expected as long as the coolant heating is not in excess of 
half of the temperature difference measured in the segment. These 
parameters for coolant heating, which were constantly monitored 
during the test, enabled a determination of possible coolant 
temperatures and coaling mass flows with the aid of simplified heat 
balance calculations on the segment. 

A schematic illustration of the cooling cycle is depicted in 
Illustration 4,5. A pump, whose speed can be adjusted, transports 
the thermostabilized cooling water into the 28 cooling segments via a 
distributer. Because of the varying expansion of tunnel walls and 
the plate, the cooling water is supplied from the bottom. Each 
segment is assigned to a control valve which is used to adjust the 
mass flow to achieve a constant plate surface temperature. Via an 
accumulator, the heated cooling water is returned to the thermostatic 
pool. Simple cylinder flow-through restrictors are installed in the 
cooling water return lines to monitor and reproduce operating 
conditions. The differential pressures of these restrictors are 
registered on a multi -tube manometer. The chronological stability of 
the coolant inlet temperature was monitored with the aid of an x,t- 
recorder. It was t 0.2 K during the test. Additional thermo- 
elements near the feed and outlet of several segments enabled the 
above mentioned control of the coolant heating. 
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1 - multi-tube manometer; 2 - restrictors; 3 - cooling water 
preparation; 4 - desired; 5 - actual; 6 - recorder; 7 - heating; 8 
cooling; 9 - test tunnel; 10 - plate segment; 11 - heat insulation; 
12 - control valves; 13 - pump (adjustable speed) with bypass; 14 - 
distributor . 

Illustration 4.5; Cooling Cycle 
4.3 Test Methods 


4,3.1 Determining Free Stream Conditions 

Free stream conditions were tested with combined total pressure and 
total temperature probes. The applicable static pressures were 
determined with pressure holes in the tunnel walls. The measured 
values were used to calculate velocity and static temperature with 
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the known relationship for isentropic condition changes and with the 
aid of the ideal gas equation for air. When flow parameters and 
coefficients were determined, the temperature dependency of the 
material data of the flow medium was considered. 

Because, even with acceleration or retardation, total pressure and 
total temperature in the f riction-f ree free stream flow experienced 
no change, the velocity distributions .'along the plate can be 
determined with applicable static pressure paths. Because of this, 
the static pressure on each segment was measured with a pressure 
hole, while total pressure and total temperature were determined in a 
reference plane. The reference plane was located 100 mm up-stream 
from the leading edge of the plate. Position and dimensions of the 
probe (see illustration 4.6 ) were selected in such a manner that 
there was no influence on the boundary layer flow along the test 
pi ate. 
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1 - ’reference probe’; 2 - ’boundary layer probe’ (laminar boundary 

layer); 3 - ’boundary layer probe’ (turbulent boundary layer) 

Illustration 4.6: Flow Probes for Measuring Total Pressure and 

Total Temperature 
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4.5.2 Determining Heat Transfer Coe-f f icients 


The measurement method -for local heat transfer coefficients is based 
on the supposition of one-dimensional heat conductivity in the coaled 
segment. Illustration 4.7 is a schematic representation of the 

temperature path. The linear temperature distribution over the 
segment, which is assigned to the one-dimensional heat conductivity 
process, is present up to the immediate vicinity of the cooling 
hole. The linear temperature path changes near this hole, because 
the heat flow, which is to be transferred to the coolant, is 
distributed throughout the entire hole diameter. Two-dimensional 
temperature field calculations (’finite elements’) have proven that 
this area of influence of the cooling hole does not come in contact 
with the measurement plane of the lower thermo-element in the 
segment . 



1 - temperature 

Illustration 4.7: Schematic Temperature Path in the Segment 




Because of the identical heat -flows 
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gas 


= Q 


seg 


(4.1) 


the calculation of the heat transfer coefficient with the 
temperatures measured in the segment can be made by using the 
following equations (Illustration 4.7): 


Q 3 a • (T • * T ) • S • b 
w gas g v ei w' 


(4.2) 


Cu 


seg 


(T 0 - T u ) - s . b 


(4.3) 


Q is the amount o-f heat removed from the hot gas by cooling. Q 

describes the heat flow in the segment; b is the tunnel width. Both 

heat flows refer to this width. T and T are the two measuring 

temperatures in the segment, is the plate surface temperature, and 

T ^ is the ambient temperature of the wall. The relationship 

(equation 4.1) provides the following solution for the unknown heat 

transfer coefficient « : 

9 

*Cu T o T u s 

“g " T- • T — r j * 5 (4.4) 

ei w 

The wall temperature is derived from the linear extrapolation of the 
temperature path between the temperature measuring locations: 

T w = T 0 + r * (T o • V (4.5) 


Because of the short distance 1’ = 1.5 mm between plate surface and 

the top thermo-element, there is little difference between T and 

o 

T^. The ambient temperature T is determined with the recovery 
factor r: 


T ei 


T + r 


09 


U 2 

°° T 

' T t . 


^5 
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r = JPr (F'r - F'randtl number = 0.71) applies -for laminar -flows and r 
« for turbulent flows. Because, even when pressure gradients are 

present, available local Mach numbers will not exceed 0.3, inherent 
and total temperatures are almost identical. 

The assumption of one-dimensional heat conduction between the upper 
and lower thermo-element along the entire width of the segment had to 
be studied, especially in the area of the upper thermo-element. The 
reason for this was the presence of a crack in the segment width 
caused by the bars which connect the segments. 

The analysis was made with the aid of two-dimensional temperature 

field calculations in the finite element method. In order to take 

into consideration the reciprocal effect of all bars, the calculation 

was made with three inter-connected segments. The given parameters 

for the calculation were: the coolant temperature, the heat transfer 

coefficient on the coolant side, the gas temperature, and the path of 

the heat transfer coefficient along the three segments on the gas 

side. In Illustration 4.8 the result of this analysis is shown for 

two characteristic critical cases for the path of the heat transfer 

coefficient on the gas side. The illustration depicts the calculated 

temperature distribution for the center segment area in question, 

which is marked by the profile crack. T q and T are the temperatures 

calculated at the thermo-element measuring locations. The example 

« (x)~x J applies to the laminar heat transfer, while « (x) = const, 
g g 

is a good approximation for a turbulent case. Illustration 4.8 shows 

that the minor deviation (2-4% max. of the temperature difference) at 

the measuring locations from the temperature distribution during 

ideal one-dimensional heat conductivity confirms the adaption of the 

one-dimensional heat conductivity process. The reason for this is 

the high heat conductivity of copper, which is used as the segment 

material, which causes a rapid even distribution of the temperature 

field in this area. 
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Verteilung der Wdrmeubergangszahlen * 

i 



FEM Gitternetz : — 



Temperaturverteilung in Thermoelementhohe h: 



1 - distribution of the heat transfer coefficients; 2 - FEM grid 

network; 3 - temperature distribution in the height h of the thermo- 

element; 4 running length x. 

Illustration 4.8; Calculated Temperature Losses in the Cooling 

Segment with a Height h of the Top Thermo- 

Element for « (>;) = const. ( — i) and 

, . „ -0.5 9 , 

<x(x)~x ( ) 

g 




An additional examination of the measuring method was made by 

calculating the heat transfer coefficients with the aid of the 

definitive equation (equation 4.4) and the values T and T , which 

o Lf 

were obtained with temperature field calculations at the temperature 

measuring locations. The resulting « -numbers were almost identical 

9 

with the heat transfer coefficients in the segment center for the 
field calculation. 

Because of their significance for an exact determination of the heat 
transfer coefficients, the heat conductivity of the segment material 

was also measured. This was accomplished at the Nuclear Research 

Center Karlsruhe (Insitutue for Material and Solid State Body 

Research) with the Laser-Flash method whose principle is discussed in 
detail in C1021. The value for the temperature range of 290-320 K 
is X £ = 355 W/mK. This result is within the specifications which 

apply to commercial copper (see C1033). 

4.3.3 Boundary Laver Tests and Evaluations 

In order to determine the boundary layer profile, total pressure and 
total temperature profiles in the boundary layer plus the static 
pressure were measured in a wall hole within the respective test 
plane. The boundary layer probes used ( II lustration 4.6) have two 
separate, parallel probe arms. One arm serves as a Pitot tube for 
total pressure measurement, the other contains a sheathed 

thermocouple (ft = 0.25 mm) and is used for total temperature 

measurement. 

For laminar and transitional boundary layers with mostly very thin 
boundary layer thicknesses of 1.3 mm or less we used a Pitot tube 
which was reduced from an outer diameter of 0.6 mm to a height of 

0.24 mm (see Illustration 4.6). A probe with a circular opening 

proved to be more practical for fully turbulent boundary layers, 
because it permits the use of known calibration curves to determine 
tangential stress. 
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The probe is moved with a micrometer gauge. It was inserted into a 
slotted, -fixed probe shaft, which has a thickness of 4 mm and which 
extends to the immediate vicinity of the plate surface (see 
Illustration 4.2). A fixed shaft arrangement ensured that the 
effective flow cross-section in the measurement tunnel and, 
consequently, the flow conditions, would not be altered by probe 
displacement. The distance of the probe to the wall was electrically 
registered by a calibrated Dhmic displacement sensor; it was then 
compared to the reading of a precision dial gauge. The distance of 
the two probe arms from the plate surface was observed under a 
microscope against the light. 

The evaluation of the measured rough data consisted of the following: 

- correcting the total pressure profiles by considering the 
influence of the shear flow and the wall. 

- converting pressures and temperatures to velocities 
(see paragraph 4.3.1). 

- determining tangential stress coefficients and Stanton 
coefficients. 

calculating integral boundary layer quantities. 

Shear flow and wall influences occur because of the velocity 
gradients and the streamline displacement near the wall caused by the 
probe. These influences are corrected by assigning to the measuring 
area a location, which is displaced in the direction of the greater 
pressure, versus the geometric center of the Pitot tube. McMillan’s 
corrective functions C104] were used for the round Pitot opening, the 
corrections of Quarmby et al C1051 were applied to flattened probes. 


Depending on the type of flow, different procedures were used to 
determine the wall tangential stress. For laminar boundary layer 
flows it was calculated with the analytically integrated form of the 
theorem of momentum, whose derivation is described in detail in E221: 

«2 

c f a 3T (4.7) 

Cj. is the tangential stress coefficient and is obtained with the 
measured impulse loss thickness 6^ and the boundary layer runninq 
length x. This calculation method applies to incompressible and 
compressible boundary layers. 

A determination of the tangential stress with the linear velocity 
gradients at the wall was not made in the laminar case. Because of 
the low boundary layer thickness, the wall gradient could not be 
obtained with sufficient accuracy. 

The following two procedures were used for the determination of 
turbulent wall tangential stress: 

the Preston method [1053 (measurement with applied probe) 


- using the general logarithmic wall theory for boundary layers 

without pressure gradient as an implicit definitve equation 

for r i 
w 


u_ 
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y-u 


+ C 
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T 




(4. B) 


with C = 5.1; k = 0.41 


Because of the non-isothermal boundary layers with both test methods, 
the temperature dependency of the material values (density, 
viscosity) had to be considered when determining the tangential wall 
stress with a cooled plate. When using the Preston method, the 
material values were formed with the Eckert reference temperature 
which was proven adequate by Holmes and Luxton C1073. This reference 
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temperature was, in our case, almost identical with the arithmetic 
mean values of free stream and wall temperatures. 

When the wall theory is used to determine tangential stress, the 
material values were based on a temperature which was determined 
as an integral throughout the area near the wall (y < 0.2 x S). 

This is analogous to a transformation of the velocity profile of a 
boundary layer with variable density into a velocity profile of a 
turbulent boundary layer with constant density (see [76, 1083). 

Both tangential stress determination procedures provided results 
with a maximum deviation of 57.. The tangential stress values, which 
were obtained as the average of both procedures, were used as the 
actual result. 

The measured boundary layer profiles were used to define the (42 
integral boundary layer characteristics such as displacement 
thickness, impulse loss thickness, enthalpy loss thickness, etc. 

(see list of symbols). Due to the finite dimensions of the boundary 
layer probes, however, the boundary profiles can only be recorded 
beginning with a wall distance which is equal to one half the probe 
diameter. Along the direction of the wall, the measured profiles 
must be extrapolated with suitable equations to enable a correct 
integration of the profiles for the determination of integral 
boundary layer characteristics. As with the determination of wall 
tangential stress, a differentiation had to be made between laminar 
and turbulent boundary layer flow. The expansion of the laminar 
boundary layer profiles was accomplished by approximating the actual 
boundary layer profiles with the following equation by Polhausen 
[1093s 


- 1 - (1-n) 3 • (1+a-n) ; n - f (4.9) 

oo 
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For pro-files in the reversal area, the exponential equation 


u 


oo 



(4.10) 


was used with the parameters A and B, which were determined by 
adjustment to the measuring points. In the turbulent case, the 
expansion was made with the aid of the logarithmic pro-file 
illustration. In the area of 0 < y + < 5, the laminar bottom layer 

was supplemented by the -formal statement 


/ with u+ 


y-/r w 77r w 


■'hTe’w 


(4.11) 


+ 

Starting with y =30 and ending with the -first profile measuring 
point, an extrapolation was made with the logarithmic wall theory: 

u + , I . In y + ♦ C ; * » 0.41 ( or according to equation 5.7) (U.12) 

K 


The constant "C" was determined by adaptation to the measured 
profile. With wall cooling, the constant K was, in accordance with 
the temperature influence, determined with the derived temperature 
correction of chapter 5.5. 1.2 (equation 5.7a). The area between 
laminar bottom layer and wall theory (5 i y + < 30) was bridged with 
a polynomial of the 3rd order 

+ + H 2 +3 

u = a Q + aj • y + a 2 • y + a 3 • y (4.13) 

whose coefficients were determined from equations (4.11) and 

(4.12) and derivatives of the validity boundaries with y + = 5 and y + 
= 30. 
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The same method was selected -for expanding temperature pro-files and (43 

analogous equations. (T-T )/(T -T ) was used instead o-f u/u . and 
+ w ot w «’ 

u was replaced by the dimensionless temperature 


T + = 


( T - T w>- p w- c pw’^ 


(4.14) 


which de-fines the wall heat flowc** > the wall tangential stress r , 

* w 

and the material values c and p at the wall. 

p, w K w 

Based on the analogy between enthalpy and impulse exchange, the 
analytical equation -for the temperature pro-file in the bottom 
laminar layer is 


_+ + 
T = Pr x y 


(4.15) 


The -following applies to the logarithmic wall theory (see C361): 


T+ 1 1 r 

T = — • In y + C T 


(4.16) 


with 

Pr t = 0,9 
r T = k /Pr t ; 

C-j. was determined with the aid of measured profiles. 
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4.3.4 Test Values: Pressures and Temperatures 


ORIGINAL RAGE IS 

OF POOR dUALFTY 


A summary of all pressure and temperature measurement locations is 
provided in Illustration 4,9. It is a schematic of the measuring 
plate equipped with thermal elements and static pressure holes, the 
positioning of the reference probe to determine free stream data, 
the static pressure holes in front of the leading edge of the plate 
to adjust boundary layer suction, and the measurement connections at 
the boundary layer probe. A PDP 11-34 by ’Digital Equipment 
Corporation’ with a 128 K core memory was used to register pressures 
and temperatures. With a "Scanivalve" alteration switch, the 
computer can monitor 48 pressure locations; with a relay 
multiplexer, it is able to monitor a total of 66 stress and 
temperature measuring locations. 



1 - air supply; 2 - ref. probe; 3 - 

equalizing distanced- measuring section; 5 - 

outlet; 6 - air heater; 7 - boundary 

layer probe; 8 - segments; 9 - pressure 

measurement; 10 - temperature measurement; 11 - 

relay multiplexer; 12 - computer. 

Illustration 4.9; Schematic of the Test Unit with Measuring Value 

Registration 










4.3.5 Determining Free Stream Turbulence with Laser-Doppler 
Anemometrv 


The Institute for Thermal Dynamics at Karlsruhe University has been 
using and developing optical methods for flow measurements for a 
long time. This applies specifically to the Laser-Doppler 
anemometry (LDA) technique for determining turbulence intensity 
paths throughout the measuring plate. The measuring system used is 
based on comprehensive development work by S. Eriksen and S. 
Wittig. They have introduced this systemin a number of publications 
C79, 80, 110, 111]; it is therefore described only in a general way 
in this dissertation. 



1 - beam separator; 2 - photo multiplier; 3 - 

lens ; 4 - oscilloscope 


Illustration 4.10: Laser Doppler Anemometer (LDA) 







4.3.5. 1 Description of the LDft Procedure 


A schematic of the Laser-Doppler anemometer is shown in II lustration 
4. 10. The unit consists of a 4 W argon/ion laser whose beam is 
directed on the transmitter ’s optical system after a deflection of 
180°. The optical system consists of a beam separator module which 

splits the beam into two parallel beams of equal intensity with a 

distance of 30 mm. A lens with a focal length of 500 mm is used to 
cross and focus the two beams. At this location, the measured 
volume has a length of 4 mm and a diameter of 0.22 mm. 

As shown in Illustration 4.10, the stray light signals are 
registered in a forward scattered direction. The stray light is 
focused by a lens, which also has a focal length of 500 mm, on a 
shutter. A photo multiplier is located behind the shutter to record 
the signals. The shutter will reduce the effective measured volume 
diameter to 0.1 mm. 

4. 3. 5. 2 Evaluation of the Measuring Signals 

For the recording and evaluation of measured signals, a system 
consisting of counter, interface and computer was used (see 

illustration 4.10). With aid of the fast counter (500 MHz, TSI 

£1123) the Doppler frequency (see £110, 1113), which forms the (46 

measuring signal, was determined. The statistical evaluation was 
made on-line with a 128 K PDP-11/34 16 Bit computer. The interface 
adjusts the counter to the computer; it also serves as an 
intermediate storage and control unit for the transmission of 
digitized measured values to the computer. The transmission speed 
determined by the computer is 1 kHz. 
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For statistical evaluation, individual measurements were made for 
each position of the measured volume 7,500. These measurements were 
then used to calculate values, which were determined 
chronologically, i.e, average variation and average variation. 


The velocity determined chronologically was calculated in accordance 
with the following formulas 


N 



whereby t gi represents the real measuring time. It could accurately 
be determined by the evaluation system. In this procedure, the 
average value of the turbulent variation quantities is equal to the 
standard deviation of the measured velocities c., which is derived 

i* 

from the following equation: 



l t B 1 < 4 . 18 ) 


The turbulence intensity of a flow was defined as the average value 
of the velocity vector for the three variation components 


Tu 


(4.19) 


2 2 

u’ , v’ were determined in acco rdanc e with equation (4.18)5 based 

2 2 

on a symmetrical tunnel inlet flow, w’ ■ v’ was used. 


The LDA optic is a single component structure. Nevertheless, the 
velocity component in the main flow direction as well as its 
vertical component can be determined in two-dimensional flows. The 
following equations C i 133 apply to average velocity and standard 
deviation with an angle * to the main flow directions 
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c x - u • cos* + v • sin* (4.21) (47 

<P 

c ' 2 = u' 2 »cos 2 * + 2 uV • cos <j) • sin* + v ' 2 • sin 2 * 

.* I (4.20) 

By measuring r J and e ' 2 i in three different directions in the flow 

<P - | *- 

field, these two equations can be used to calculate the five unknown 

values on the rig ht s ide of the equat ion. In order to determine the 

2 2 

turbulence values u’ , -u’v’ and v' , I therefore used angles of 4 = 

-45°, 0° and +45® to the main flow direction in my study. 

4. 3. 5. 3 Particle Preparation and Particle Feed 

The LDA procedure evaluates scattered light signals of particles 

contained in the flow. In order to achieve reliable data, it was 

necessary to add particles (Si0„ solid particles, AEROSIL R972, Si0_ 

Z Z 

portion = 99.8% C1141) to the flow. 

The particles, which are supplied as a powder, are round and have a 
diameter of approximately 20 nm. To prevent these primary particles 
from agglomerating, which would cause them to no longer follow the 
flow motion in a free-flowing manner, Eriksen et al C110, 1113 
developed a comprehensive particle preparation system where, with 
the aid of a fluid bed storage container, the particles are 
constantly aggitated and di sagglomerated. The particle distribution 
volume achieved in this manner showed particle diameters between 0.1 
and 2 Fm and their maximum frequency was between 0.3 and 0.5 Fm. 

This fulfilled the required conditions for a free-flowing 

determination of the turbulence volume. The volume distribution of 
the particles was examined with several procedures (cascade impactor 
(Anderson sampler), Whisker particle collector, multi ratio 
individual particle counter) (see C 1 10, 1113). 
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Detailed information on this particle preparation and the measuring 
system used can be read in the publications by Eriksen and Wittig 
C110. 1113. 

Comprehensive studies were also made to determine the best location 
for particle insertion into the flow. The pipe used for particle 
insertion was placed immediately behind the last stabilizing grid of 
the delay path (Illustration 4.1). With this arrangement, no 
interference of the flow in the test tunnel was caused by the 
particle insertion. 

4.4 Test and Calibration Measurements 

4.4.1 Calibration of the Measuring Plate Thermal Elements 

An important prerequisite for accurate determination of local heat 
transfer coefficients was the exact calibration of the thermal 
elements on the measuring plate. This was especially important in 
my case because, due to the high heat conductivity of the copper 
used for the measuring plate, temperature differences of less than 1 
K had to be anticipated with the low gas temperatures in the 
segment . 

Illustration 4.11 is a schematic of the arrangement used for the 
calibration of the thermal elements. Together with a resistance 
thermometer as a standard, the thermal elements were imbedded in a 
copper body which, in turn, was placed in a closed-off thermostat 
tank. With this arrangement, the thermal elements and the 
resistance thermometer have the same temperature; possible 
temperature variations in the thermostat tank are balanced. The 
calibration temperature was determined with the voltage decrease at 
the resistance thermometer and with an additional reference 
resistance which was connected in series (see Illustration 4.11). 



ORIGINAL PAGS-tl9 
OF POOR QUALITY 



1 - computer; 2 - DVM connection box; 3 - 

ice element; 4 - resistance thermometer; 

5 - rotary transformer ; 6 - Cu-body; 7 - 

recorder; 8 - reference resistance; 9 - 

voltage source; 10 - 

calibration temperature. 

Illustration 4.11; Calibration of Thermal Elements 

Calibration was made with the aid of a computer controlled (49 

measurement value recording unit. Four sets of 16 each Ni-Cr-Ni 
thermal elements were calibrated for the 0-160® range. The 

calibration curve (polynomial of the 3rd order) was then calculated 
for each thermal element. All calibration curves were very similar; 
therefore, the thermal elements, after separating several extreme 
cases, could be summarized into four groups with one calibration 
curve each. In order to examine the calibration accuracy, test 
measurements on 32 of the 64 thermal elements were performed. With 
the calibration curves obtained, temperatures of the thermal 
elements deviated by less than 0.04 K at a temperature of 

approximately 20®C; the deviated factor doubled to 0.08 K with 
approximately 130® C. 
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4.4.2 Test Measurements at the Boundary Laver Tunnel 


Test and calibration measurements serve for the determination of the 
flow properties of the measuring tunnel. The homogeneity of the 
free flow and the achieved two-dimensionality of the plate boundary 
layer is of special interest to us. Additionally, the laminar 
boundary layer development at the leading edge of the plate must be 
confirmed and the measured heat transfer coefficients must be 
examined. 

The distribution of free stream temperature and velocity throughout 
the height and width of the tunnel can be seen in Illustration 4. 12a 
and b. These results, which were achieved with a traversing probe, 
were applied in a dimensionless manner, whereby the temperature 
paths are referred to the temperature difference between free stream 
and wall. The paths shown were measured with a free stream 

temperature of 370 K, a wall temperature of approximately 300 K, and 
an average free stream velocity of 47 m/s. Illustration 4.12a shows 
the profiles, which were measured in the center of the tunnel, down- 
stream of the leading edge of the plate through the vertical 
dimension of the tunnel. Illustration 4.12b shows the temperature 
and velocity distribution through the width of the tunnel at two 
distances from the leading edge of the plate. This was done to 
demonstrate the development of the side wall boundary layers. The 
low variations of less than 17. and the small portions of the side 
wall boundaries indicate a free flow which is very smooth. 

With the shaped leading edge of the plate, a laminar boundary layer 
clearly existed at the beginning of the plate. This was determined 
by boundary layer measurements with low free stream turbulence in 
the first measuring plane <x = 120 mm). With isothermal boundary 
% layer flow, the profile paths measured are in good agreement with 
Blasius’ C1091 analytical solution. Examples, which support this 
solution, are shown in the following evaluation of measured laminar 
profiles, chapter 5.5. 1.1 (e.g. illustration 5.26a). 
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1 ~ tunnel height; 2 - tunnel center; 3 

tunnel width. 

Illustration 4.12: Temperature and Velocity Distribution through 

Tunnel Height (a) and Tunnel Width (b) in 
Different Measuring Planes <x=15, 80, 400 mm) 
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Indications of the two-dimensionality of the laminar and turbulent 
boundary layers were obtained when determining the tangential stress 
distribution diagonal to the flow direction; these determinations 
were made of varying distances from the leading edge of the plate. 
The measurements were made with a Preston tube (# « 0. 6 mm), which 
was introduced to the tunnel from the side and which was moved in a 
traverse direction on the plate. Results, which were obtained with 
two temperature ratios and flow velocities at 100 and 400 mm of 
plate running length, are shown in Illustration 4.15a and 4.15b . 
Illustration 4.13a shows measurements for a plate flow with a 
completely laminar boundary layer, 100 mm running length and a fully 
turbulent boundary layer at the plate end with 400 mm. Because the 
Preston method cannot always be used for laminar and reversing 
boundary layers, the velocities Uy o , which are calculated from the 
dynamic pressure on the probe, are applied to the free stream value 
and are illustrated as such. Illustration 4.13b shows the results 
for the fully turbulent boundary layer, which were achieved at the 
same distances from the leading edge of the plate. The advanced 
reversal, which occurs with the same flow, is a result of increased 
free stream turbulence. Maximum variations in the laminar and in 
the turbulent boundary layer were below 57. (Illustration 5.13a and 
b) , which sufficiently fulfilled the requirements for two- 
dimensionality. 

In order to evaluate measured heat transfer coefficients, enthalpy 
balances were made which compare the heat transferred to the plate 
with increased thickness in the boundary layer due to enthalpy 
loss. As these balances are suitable for the examination, analysis 
and interpretation of the results, they will be discussed in the 
following chapters (chapter 5.4). 

TR.Note: German pages 52-53 missing 
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Nr. 


dp/dx 

T 

00 ,0 

T w /T »,o 

u «,o 

r *L,o 

[K] 

[-] 

[m/s] 

[-] 

1 . 

k 

- 0,0 

330 

1,0 

41 

1,0 

1.1 


0,0 

350 

1,0 

40 

0,8 

2. 


0,0 

378 

0,8 

47 

0,9 

2.1 


0,0 

375 

0,8 

41 

0,76 

3. 


0,0 

478 

0,64 

55 

0,65 

4 . 


0,0 

570 

0,53 

38 

0,36 

- KONTUR 

1 

k(x) - const. 




5. 

k 

= 1,2-KT 6 

381 

0,8 

35 

0,64 

6. 


1,2-KT 6 

466 

0,64 

49 

0,64 

7. 


2,3- 10" 6 

480 

0,64 

27 

0,33 

8. 


2,5- 10" 6 

571 

0,53 

33 

0,32 

.KONTUR 

2 

k ( x ) = const . 




9. 

^max 

= 2.9- 10 -6 

376 

0,8 

39 

0,73 

10, 


3,2- 10" 6 

463 

0,64 

48 

0,65 

11. 


5,65-KT 6 

385 

0,8 

20 

0,35 

12. 


5,7-10-6 

467 

0,64 

27 

0,35 

13. 


6,3- 10" 6 

570 

0,53 

34 

0,33 


1 - contour 

Table 5. li Test Parameters (applicable velocity paths are shown 
in II lustrations 5.45, 6.15 and 6.19) 

The discussion of results for heat transfer and boundary layer 
measurements is divided into individual paragraphs, which are in 
accord with normal experimental procedures. 

The first paragraph (Chapter 5.3) shows the results of boundary 
layer measurements without pressure gradients or wall cooling for 
varying free stream turbulence intensities. These test measurements 
are used to compare results, which are obtained from literature, and 
provide a starting point for an analysis of additional measurements 
with wall cooling and pressure gradients. 
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The -following chapters will introduce the new heat transfer and 
boundary layer measurements with varying wall cooling (T /T = 0.8; 
0.64; 0.53) for unaccelerated flow (chapter 5.4) and for accelerated 
flow (chapter 5.5). Free stream turbulence was a variable in these 
tests. 

Measurement data obtained under hot gas conditions provides new 
information, which expands the present state of knowledge, of the 
significance of the individual influencing parameters of free stream 
turbulence, wall cooling and pressure gradients, and about their 
effects when these parameters are superimposed. Analyses on this 
subject are summarized in chapter 5.5. They describe fundamental 
influences on reversal events, heat transfer rates, tangential 
stress coefficients, boundary layer characteristics, and profile 
deformation in the boundary layer. 

5.2 Free Stream Turbulence 


Free stream turbulence is one of the most important parameters of 
boundary layer flows. Knowing its exact dimensions and path over 
the measuring plate is an important prerequisite for an accurate 
evaluation of the boundary layer and heat transfer data obtained 
from the experiment. In order to meet these requirements, 

comprehensive turbulence intensity measurements were performed with 
the aid of the optical Laser-Doppler measurement technique. Other 
tests made under the same flow conditions are heat transfer and 
boundary layer measurements with unaccelerated and accelerated 
flows. 

In order to determine the turbulence intensity, the u’ and v’ 

components were me asur ed in a nd perpendicular to the main flow 

2 2 

direction, whereby w’ = v’ . Because the tunnel inlet flow was 
symmetrical in the transverse direction, it was assumed that the 
development of both transverse components of free stream turbulence 
occurs in the same manner and, therefore, their values are almost 
identical . 



Illustration 5.1 shows - for the situation without pressure 
gradients - turbulence intensity as a function of the plate’s 
running length, measured without the four installed turbulence 
grids. The zero point is defined as the leading edge of the plate. 
The measurement values were determined for 5 wall distances (Ay = 5 
mm) within the flow area with a width of 25 mm which is in contact 
with the outer edge of the boundary layer. The entire temperature 
range up to 300° is included in the results. The free stream 
velocity varies between 39 and 55 m/s, and the Reynolds number, 
which is a function of the plate length, is between 0.36 x 10 6 and 
1.0 x 10 6 . 

The measurement values describe turbulence intensities at the 
beginning of the plate. Without a grid, they are 1.67. and with 

grids 1 through 4 they are 2.37., 3.87., 6.5% and 8.7%. As expected, 
the turbulence intensities showed no influence of the free stream 
temperature. This is confirmed by comparative tests, which were 
conducted with the grids installed and with an identical velocity 
but with different free stream temperatures. 



1 - without a grid; 2 - grid 


Illustration 5.1: Turbulence Intensity Matrix in an 

Unaccelerated Flow (Temperature Range T w = 
70-300° C, Leading Edge of the F'late: x = 0) 




The attentuation behavior, which is typical for grid turbulence, is 
present in these measurements. A characterization of this event is 
demonstrated in Illustration 5.5 . In this illustration, the data 
are applied as a function of the dimensionless distance to the 
turbulence grid x/b. The parameter b is the wire strength of the 
grid. As a comparison, the curve with its applicable spread 
according to Baines and Peterson C 101 3 was added to the 
illustration. It serves as an indication for typical grid 
turbulence and is confirmed by my own measurements. 

The spread of the measured turbulence intensity shown in 
Illustration 5.1 is based on the higher variations of the transverse 
component v’ . This becomes obvious in Illustration 5.3 by the 
application of individual turbulence intensity components. It 
applies especially to grid 4, which is the coarsest grid. This 
effect is supplemented by imbalances in the turbulence intensity 
distribution over the wall distance. Due to the small distance from 
the grid, these imbalances can be up to 17. with the high turbulence 
intensities of grids 3 and 4 (in the u’ and in the v' component). 
Additionally, there is no distinct main flow component in the 
transverse direction under the given flow situation. This makes a 
test solution of the turbulence component v’ more difficult and 
results in additional inaccuracies of the measurement result. 
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1 - grid; 2 - Baines and Patterson 

spread. 

Illustration 5.2i Logarithmic Application of the Turbulence 

Dependency Matrix Independent o-f the 
Dimensionless Distance x/b to the Turbulence 
Grid 

In the case of grid 4 it must be noted that, compared to the other 
grids, this very coarse grid has a different obstruction factor 
(0.44 versus 0.64). Because of this, it is assumed that it 
generates turbulences, which are increasingly dependent on the inlet 
conditions. The basis for this presumption is that a variation of 
the turbulence intensity component was observed to be dependent on 
the free stream velocity for this grid only. The greatest influence 
was observed at high velocity. The turbulence became almost 
isotropic with decreasing velocity (Illustration 5.3). 
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All turbulence measurements with pressure gradients were made with a 
•free stream temperature o-f T w = 100°C. However, in order to 
characterize the spread o-f the values, measurements were performed 
•for all relevant operational Reynolds numbers and, consequently, -for 
all acceleration intensities. This was sufficient, because, after 
the measurements with pressure gradients, dependency on the free 
stream temperature was not given. 

Illustration 5.4 shows the turbulence intensity matrices for 
accelerated flow (’contour on the ’pressure side’ or contour 1) with 
an almost constant acceleration values. In spite of varying 
acceleration intensities, identical results were achieved for k = 
1.2 x 10 6 and k = 2.5 x 10 ^ (varying acceleration intensities were 
generated by an inlet Reynolds number variation). However, when 
compared to the results, differences exist in unaccelerated free 
flows (see Illustration 5.1): a distinctly more rapid decrease of 

free stream turbulence is caused by accelerations. 

The results for the case of the contour on the ’vacuum side’ 

(contour 2) with higher tunnel profile constriction in the forward 

portion of the plate are shown in Illustration 5.5. It also 

contains a summary of all measurements with k = 2. 9-6. 3 x 10 

max 

Compared to Illustration 5.4, an even more rapid decrease of free 
stream turbulence is noted in the forward portion of the plate. At 
the end of the highly accelerated zone (x = 250 mm) the turbulence 
decrease is again in accordance with the unaccelerated flow (see 
Illustration 5.1). 
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grid. 


Illustration 5.3: Turbulence Component Distribution in 

Unaccelerated Flow 
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Illustration 5.4; Turbulence Intensity Matrices in the 

Boundary Layer Tunnel; Acceleration with 
Contour 1 (Contour on the ’Pressure Side’) 



1 - grid. 


Illustration 5.5; Turbulence Intensity Matrices in the 

Boundary Layer Tunnel; Acceleration with 
Contour 2 (Contour on the ’Vacuum Side’) 




5.3 Boundary Lavers on Adiabatic Walls 


OR!G!N*U 

OF POOR QUALITY 

In accordance with experimental procedures, this paragraph will 

introduce the measured results on a flat, uncooled plate. They will 

be used for a basic characterization of available boundary layer 

types and for the demonstration of the significance of free stream 

turbulence independent of wall cooling. The measurements refer to 

free stream velocities of 40-41 m/s and to Reynolds numbers of Re. . 

(h L ■ U 

= 0.8-1. 0 x 10 , which were achieved along the entire running length 

of the plate (see Table 5.1, flow conditions 1 and 2). With low 

Reynold numbers, the respective boundary layers are already showing 

characteristics which are relevant for a turbine. In addition to 

reversal events, which can be found on turbine blades, there are 

typical relaxation effects (low Reynolds number effects) in 

turbulent boundary layers. The measurements will characterize these 

effects by considering the additional influences of free stream 

turbulence, which is typical in a turbine. 



1 - surface. 

Illustration 5.6; Velocity Profiles and Boundary Layer 

Development along a Flat, Adiabatic Wall 




5.3.1 Transitional Boundary Lavers with Varying Free Stream 
Turbulence 

All boundary layers, which were observed in the measuring section, 

were transitional. Their appearance at the measurement surfaces, 

which were detected with probes or heat transfer measurements, was 

predominantly defined by their Reynolds number and the respective 

free stream turbulence. Illustration 5.6 shows two critical cases 

of boundary layer forms and boundary layer development caused by low 

and high free stream turbulence (Tu = 1.6% and 6.5%). The 

o 

measurements of the lowest possible free stream turbulence in the 

measuring tunnel (Tu = 1.6%) describe events above the laminar- 

o 

turbulent boundary layer reversal . The profile in surface 1 (x = 

0.11 m) is a pronounced laminar profile. The matrices in surfaces 1 
and 3 (x = 0.21 m and 0.22 m) illustrate reversal profiles. A fully- 

turbulent boundary layer profile is present in surface 4. In 

contrast, all profiles measured with grid 3 (Tli q = 6.5%) are 

turbulent. 

The development of both boundary layers with different external 

turbulence is accompanied by a varying boundary layer thickness 

increase and changes in the profile shape. These changes not only (62 

exist in the transitional case (Tu = 1.6%) - because of the 

o 

reversal - but also in the turbulent flow with grid 3 (Tu = 6.5%). 

o 

A comparison of both profile developments in Illustration 5.7a and 
b clearly shows these events. The velocity profiles, which are 
identified with the wall tangential stress velocity u r , are applied 
in a semi-logarithmic illustration. As a comparison, the 
dimensionless linear velocity theory u + = y + (Illustration 5.7a) is 
shown for the laminar profile in surface 1 and the respective 
logarithmic wall theory is shown for turbulent profiles. In a 
diagram, Illustration 5.7 summarizes the turbulent boundary layer 
profiles of the 3 measured surfaces by displacing the u + coordinate. 
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In Illustration 5.7a, a step-by-step description is provided of the 
transition -from the laminar (surface 1) to the f ul ly-turbulent 
boundary layer profile form (surface 4). In surface 1, the 
measurement confirms the linear velocity distribution u + * y + near 
the wall for the laminar case. In surface 4, achievement of the 
turbulent boundary layer condition is documented by the good 

agreement of the measured data with the logarithmic wall theory in 
the area near the wall (y + < 200). 

In comparison to the fundamental profile deformation of the boundary 
layer during reversal (Illustration 5.7a), only an increase of the 
wake portion can be recognized with progressive running length in 
the laminar case (Illustration 5.7b). The wake increase illustrates 
a basic property of turbulent boundary layer profiles with 

relaxation. This condition can also be observed in turbine flows. 
It describes fundamental deformations of the profile path in the 
outer portion of the boundary layer, which are not completed until 
sufficiently high Reynolds numbers have been achieved. In boundary 
layers without external turbulence, the relaxation behavior which 
occurs with the increase of the wake component is a sole consequence 
of the Reynolds number effects. In this case, similar to blade 

boundary layers, turbulence intensity influences are an additional 
factor. This additional effect will be discussed in detail in 
paragraph 5.3.2. In addition, it will be demonstrated (para. 5.5) 
that these relaxation effects, which are superimposed by external 
turbulence, will also have an effect on heat transfer and the wall 
tangential stress. 

The profile deformations of two boundary layers with high and low 
free stream turbulence described in Illustrations 5.6 and 5.7 are 
based on characteristic differences of the boundary layer value 
(impulse loss thickness), (form parameter), and the coefficient 
of friction c^. These parameters are not addressed in this chapter. 
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Illustrations 5.8 and 5.9 show the boundary layer characteristics 

•for available measurements. The tangential stress coe-f-f icients of 

the reversal profile for Yu = 1.6% on measurement surfaces 2 and 3 

o 

were approximated for that velocity gradient which was formed with 
the measurement point closest to the wall. 




Illustration 5.8; Comparison of Transition and Turbulent 

Boundary Layer Characteristics (see 
Illustration 5.7 for Symbols) 
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Compared to the laminar comparison curve by Blasius C1093, the 
reversal process (Tu q = i,s *) is characterized by an impulse loss 
thickness increase, the decrease of the form parameter = 2.6 for 
laminar profiles to values around 1.4 (Illustration 5.8) which are 
typical for turbulent boundary layers, and an increase of the local 
tangential stress coefficients to magnitudes which are above the 
values of the laminar boundary layer by a factor of 3-4 
(Illustration 5.9). The higher turbulence of grid 3 (Tu q = 6.57.) 
will move these reversal events up-stream, to the vicinity of the 
leading edge of the plate. Therefore, the respective measurement 
values in Illustration 5.8 describe only turbulent paths which, 
compared to laminar boundary layers, are characterized by a higher 
increase of the impulse loss thickness and form parameters with 
values around 1.4. 



Illustration 5.9: Coefficients of Friction 




In Illustration 5.9, the coefficients of friction c^ are applied 

depending upon the impulse loss thickness Reynolds number. For a 
general evaluation, they are compared with known, semi-empirical 
correlations by Blasius C1093, Prandtl C22D and Kays C863 for 
laminar and turbulent boundary layers with a constant-density flow. 

My own measurements were in good agreement with these correlations. 
Average results for grid 3 (Tu q = 6.5*/.) are about 6-10’/. above the 

value which, according to Kays, applies specifically for turbulence 
free external flow. These increase of the c^-values can already be 
attributed to the influence of free stream turbulence. 

Reversal events always represent flow processes which were formed 
three-dimensional ly. However, current results relate to 
measurements which only describe boundary layer events in the center 
of the tunnel. In order to have an understanding of the three- 

dimensional formation of the reversal area, transverse measurements 
were performed; they are contained in Illustration 5.10 . In this 

illustration, velocity distributions uy Q (z), which were obtained 
when low free stream turbulences (Tu o = 1.67.) were applied. They 
were measured diagonal to the flow direction with a Pitot tube, (66 
which had a thickness of 0.6 mm and which was in contact with the 
plate surface. Measurements were taken at varying distances from 
the leading edge of the plate (x = 130, 205, 255, 320, and 380 mm). 

A laminar boundary layer, which is predominantly two-dimensional, 
occurs when x * 130 mm. The subsequent measurement surfaces 

indicate a three-dimensional, symmetrical reversal development 
whereby, in the tunnel center and on the sides, zones develop where 
the beginning of the reversal occurs at various running lengths. 

The area of highest velocity (tunnel center) can be compared to 
increased tangential stress values, thus representing a stage of 
advanced reversal process. When x = 380 mm, the end of the reversal 
process is indicated. In accordance with Illustration 5.10, the 
boundary layer regains its two-dimensional character. In the plate 
center, no noticeable increase of the velocities can be noted in 
comparison to the values of the previous measuring surface. 
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1 - tunnel center 

Illustration 5.10; 3-Dimensional Reversal Development 

(Tu = 1.6%) 
o 

Because, in the tunnel center, the reversal with the given velocity 
distribution is visibly separated -from the fringe areas, a guarantee 
is provided that the reversal events in the plate center, which are 
critical to the results, are not determined by the corner boundary 
layers at the contact points between tunnel side wall and measuring (67 
plate, but by the free stream turbulence itself. 

In addition to the diagonal measurements, Illustration 5.10 shows 
the measurements uy Q <x), which were detected by a probe displacement 
mechanism specially designed for this procedure. The measurements 
were taken in the tunnel center by continuously moving the probe in 
a longitudinal direction between x = 80 mm and x = 300 mm. The 

reversal process in the tunnel center is sketched by the 
longitudinal movement of the probe. With the minimum registered 
velocity values (uy /u M = min!) it was possible to approximate the 

O w 

location of the reversal initiation with the respective reversal 
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Reynolds number and to make a comparison with available results -from 
previous measurements. The reversal Reynolds number -for the -free 
stream turbulence of 1.6% is Re^ = 3.5 x 10®, which is in good 

agreement with the compiled values from literature by Abu-Ghannam. 

5.3.2 Turbulent Boundary Lavers with Changeable Free Stream 
Turbulence 

Due to the low acceleration intensities in the rear blade section on 

the vacuum side of gas turbine blades, there frequently are expanded 

areas with turbulent boundary layer flows. It is therefore 

necessary to know the influence of the free stream turbulence on 

turbulent boundary layers in these areas of the blade. 

Illustrations 5.11 and 5,13 show how this turbulence influences the 

form of turbulent velocity profiles, as it was determined under 

ideal conditions of a flat, uncooled and unaccelerated flow boundary 

layer (flow condition 2, Table 5.1). The velocity profiles were 

recorded in the 4 boundary layer measurement surfaces for three 

different turbulence intensities (Tu = 2.3%: 6.5%: 8.7%). In the 

o ’ 

case of low free stream turbulence (Tu = 2.3%), the illustration of 

o ’ 

the measurement results starts with surface 2 (x = 220 mm), because 
surface 1 (x = 120 mm) still contained distinct reversal profiles. 

Illustration 5.11 shows the velocity paths u/u w , which were measured 
with varying free stream turbulence, dependent on the wall distance 
y. In order to clarify this point, the measurements for all 3 
turbulence intensities were only entered for surface 4. 
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Illustration 5.11 show that, with increasing -free stream turbulence, 
the boundary layer thickness increases and the profile form fills 
out. However, these effects are relatively small. The velocity 
gradients at the wall increase at the same time, which results in an 
increase of the wall tangential stress. Under the given conditions, 
the higher free stream turbulence causes an overlap of the profile. 

This does not change the impulse loss thickness, but it decreases 
the displacement thickness and thus the form parameter H,„. 

Illustration 5.12a confirms these events with the aid of the (68 

matrices of the boundary layer values & 2 and H 12 , which are measured 
over the plate, and the friction coefficients. 



Illustration 5.11; Turbulent Velocity Profiles with Varying 

Free Stream Turbulence 

For a general evaluation of existing tangential stress increases, 
Illustration 5.12d again applies the measured tangential stress 
coefficients as a factor of the impulse loss thickness Reynolds 
number Re 2> A comparison with a correlation 

.25 


8i 


c^ = 0.024 x Re 2 ~ 


(5.1) 



is made. Because of its good agreement with the measurement results 
by Murlis, Wieghard (see £31]) and Coles C 1 153 (max. deviation 27.), 
this correletion is especially suitable to describe the conditions 
without external turbulence. 



7.0 
Cf * 10 3 

6.0 

5.0 

4.0 

3.0 

7.0 

Cf-JO 3 

6.0 

5.0 

4.0 

3.0 
20 

0 500 MHW tSOO 2000 2500 3000 

Re 2 

Illustration 5.12s Characteristic Quantities in Turbulent Boundary 

Layers without Heat Transfer with Varying Free 
Stream Turbulence (see 111. 5.11 for Symbols) 
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In this generally similar illustration, increases of the wall (70 
tangential stress of up to 10’/. can be observed with increasing free 
stream turbulence. 

With semi -logarithmic application of the dimensionless velocity 

profiles - Illustration 5.13 - an excellent agreement exists for all 

turbulence intensities near the wall with the known, universal wall 

theory for turbulent boundary layers with constant material values. 

Free stream turbulence thus does not change the validity of the 

logarithmic wall theory (even with low Reynolds numbers; Re^ < 

2,000). However, in the external area of the boundary layer the 

increased free stream turbulence causes a distinct decrease of the 

wake component. With our flow, this process overlaps with the 

Reynolds number effects mentioned in the previous paragraph. As 

shown by Illustration 5.13, free stream turbulence decreases the 

increase of the wake force AU/ut , which is connected with an 

r 

increasing Reynolds number (or running length). With sufficiently 
high turbulence intensities (Tu q = 8.7X), the wake portion is 

completely surpressed. 



Illustration 5.13: Logarithmic Velocity Profiles with Varying 

Free Stream Turbulence 
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Illustration 5.14: Reynolds Number and Turbulence Intensity 

Dependencies of the Wake Portion (AU/u^ = 

Wake Force at the Boundary Layer Edge 

Illustration 5.14 quantifies this Reynold number dependency of the 
wake component with and without free stream turbulence. With the 
aid of my own measurements and the measuring data for turbulence- 
free external flows by Murlis C333 it was determined that Reynolds 
number dependencies below Re„ = 2,000 are especially pronounced. 

They are thus increasingly connected to profile deformations, if the 
wake portion is not completely eliminated by very high external 
turbulence (as in this case, with Tu q = 8.7’/.). Additional profile 
deformations in the wake area in excess of Re 7 = 2,000, which - 
with turbulence-f ree external flow - cease when Reynolds numbers Re,, 
are around 5,000 (see C333), seem comparatively small. Reynolds 
number effects are thus concentrated in the area of the impulse loss 
thickness Reynolds numbers below 2,000. The additional fact in 
Illustration 5.14 is that, with very high free stream turbulence 
(Tu q = 8.77.) there is no increase of the wake component increase as 

a factor of the Reynolds number; this leads to the conclusion that, 
under these conditions, free stream turbulence dominates over 
Reynolds number effects in the external area of the boundary layer. 
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The studies with adiabatic boundary layers created a basis of 
comparison, which confirms known results from literature (see 
Chapter 2) about the influences of free stream turbulence on 
transitional and turbulent boundary layers. These results, together 
with additional knowledge, are applied to three-dimensional reversal 
behavior and Reynolds number effects with free stream turbulence, 
which is also relevant to turbines. Based on this, I would now like 
to introduce and analyze the tests with wall cooling and pressure 
gradients, which are the subject of this dissertation. 

Special significance is assigned to a comprehensive documentation of 
measured data. Without them, a verification of the calculation 
program (and an examination of the developed integral procedure) 
would not have been possible. 

5.4 Boundary Laver Flows with Wall Cooling 

To date, only very few detailed measurements on the influence of 
high wall cooling in sub-sonic flows have been performed. In 
addition, such measurements were only made for turbulent boundary 
layers. This includes the studies by Back et al E753, and the 
measurements made by Sill C773. In contrast, the significance of 
wall cooling in transitional boundary layers has not yet been 
experimentally studied. 

Because of this, I conducted comprehensive basic experiments on the 
subject of heat transfer on cooled surfaces with transitional 
boundary layers. 

In order to achieve realistic cooling conditions as they are found 

in cooled gas turbine blades, the free stream temperature was 

increased step-by-step to 100, 200 and 300®C (flow conditions 2-4, 

Table 5.1) with almost unchanged wall temperature (T = 25 - 30°C). 

w 

Heat transfer distributions with and without all available 
turbulence grids were measured for each free stream temperature. 
Characteristic cases with high and low free stream turbulence were 
supplemented by boundary layer measurements in the 4 measuring 
surfaces. 


c-^ 
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5.4.1 Boundary Laver Development and Heat Transfer with Varying (73 

Mall Coaling and Free Stream Turbulence 

Heat transfer measurements on a flat plate provide the best 
information on the influence of varied free stream turbulence and 
wall cooling on the boundary layer reversal and on the heat transfer 
itself. Illustration 5,15 shows Stanton numbers for various 

turbulence intensities with moderate wall cooling (T TV = 0.8) and 

w ot 

constant free stream velocity as a function of local running length 
Reynolds numbers (Re^ = 0.8 x 10 6 ). Illustrations 5.16 and 5.17 
show comparable measurements with more intense cooling conditions 
with T /T = 0.64 and T /T = 0.53. Because these flows register 

lower Reynolds number ranges (Re = 0.64 x 1G 6 and 0.36 x 10 6 ), 

L, O 

the reference measuring points are increasingly shifted into the 
area of laminar boundary layer sections. 



1 - start of reversal 

Illustration 5.15: Heat Transfer and Boundary Layer Development 

during Moderate Wall Cooling (T /T = 0.8) 

W W 
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1 


start of reversal 


Illustrations 5.16 and 5.17: Heat Transfer and Boundary Layer 

Development with Intensive Wall 
Cooling with T /T^ = 0.64 and 0.53 

3 Ul 








As a comparison, all three illustrations (Illustrations 5.15, 5.16, 

5.17) show the matrices of the Stanton numbers according to two 
relationships by Kays C643, which describe the laminar and turbulent 
heat transfer in plate boundary layers with constant material values 
and negligible free stream turbulence. The relationship of the 
laminar boundary layer takes into consideration the uncooled leading 
edge of the plate, which has a length of 15 mm. The dotted lines (75 

are the measuring points of the laminar boundary layer with lowest 
free stream turbulence (Tu q = 1.6%). Their purpose is a general 

evaluation of the laminar heat transfer and a determination of the 
Reynolds numbers at the beginning of the reversal (see equation 5.2, 
chapter 5.5. 1) . 

Measurements taken during the three temperature conditions basically 
show identical characteristics of the boundary layer behavior. As 
expected, with increasing free stream turbulence Tu^, an advancing 
of the laminar-turbulent reversal initiation and a shorter reversal 
length were observed. In the case of greater turbulence intensity 
(Tu q > 6.5%), the reversal range is already in the vicinity of the 

uncooled leading edge of the plate. 

A special case occurs in the tests with highest wall cooling 

(Illustration 5.17): the reversal matrix is the same for Tu =1.6% 

o 

and Tu q = 2.3%. This phenomenon must not be construed as a 

contradiction to fundamental up-stream movement of the reversal with 
increasing free stream turbulence. It occurs because secondary 
influences of corner and side wall boundary layers generate an early 
reversal process, which - with Tu o = 1.6% - is located far down- 

stream, almost at the end of the plate. In this case, the reversal 
is co-induced by the side walls. Because the applied measuring 

method integrates the heat transfer over larger areas diagonal to 

the flow direction, the Stanton numbers are registered earlier than 

the boundary layer in the center of the plate. This is confirmed by 

the reference boundary layer measurements in measuring surface 3 (x 
* 310 mm, Re = 2.6 x 10 u ) . Contrary to heat transfer, which is 

already reversing at this location (see Illustration 5.17), the 
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boundary layer measurements in the center of the tunnel still show 
distinctly laminar profile forms and boundary layer characteristics 
(see Illustration 5.23 a and b). 

The point of intersection of the laminar fitting curve, which 
represents the case of smallest free stream turbulence, is the 
criteria for the reversal initiation. Its straight line is adjusted 
to the steeply increasing path of the Stanton number in the reversal 
area. The reversal Reynolds numbers, which are summarized in Table 
5.2. were thus determined for the turbulence intensities Tu q = 1.6%, 

2.3% and 3.8%. A determination of the reversal initiation with 
higher turbulence intensities was not possible, because it occurs in 
the area of the uncooled leading edge of the plate which is outside 
the measuring area. The Reynolds numbers at the end of the reversal 
were taken from the matrix of the Stanton numbers. The results are 
contained in Table 5.2. 

The analysis of this reversal data, which was determined under 
cooling conditions similar to those found in gas turbines, is 
discussed in Chapter 5.5.2. This analysis requires comprehensive (76 
comparisons with reversal measurements performed for adiabatic 
boundary layers and the inclusion of boundary layer characterstics, 
which are present in the reversal (impulse loss thickness Reynolds 
numbers) . 
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Table 5.2s Reversal Reynolds Numbers 
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The measured heat transfer distributions (Illustrations 5.15 - 5.17) 
showed that, in agreement with previous measurements (see Chapter 
2.2), -free stream turbulence in unaccelerated flows has no 

significant influence on laminar heat transfer. I would like to 
point out that this applies only in the unaccelerated case. With 
accelerated flow there is a significant influence of free stream 
turbulence in laminar boundary layers (see Chapter 5.5.3). 

In the areas of turbulent flow there is an indication of heat 
transfer coefficient increases as a function of turbulence 
intensity. However, a distinctly even increase with increasing 
turbulence is not noted with the given spread of the measuring 
values (@ 37.). In the case of highest wall cooling (T /T = 0.53; 

Illustration 5.17) the local turbulent heat transfer coefficients of 
varying free stream turbulence are completely identical in their 

path. No influence of free stream turbulence whatsoever is thus 

noted. 

This ambiguous behavior was observed in previous studies (see 
Chapter 2.2) and resulted in different evaluations of the influence 
of free stream turbulence on turbulent boundary layers. It was (77 

frequently not considered that, with the increase of free stream 
turbulence, the laminar-turbulent boundary layer reversal takes 
place further up-stream. This results in additional consequences 
for the subsequent turbulent boundary layer. Because the reversal 
occurs earlier, there are greater boundary layer thicknesses with 
increased heat transfer resistance, although the running length is 
the same. Because of this, no or only a minor net effect becomes 
obvious when the heat transfer coefficients are observed locally. 

Other reasons are the Reynolds number effects, whose importance for 
the formation of turbulent boundary layer profiles has already been 
discussed, and characteristic dependecies of turbulent boundary 
layers on free stream turbulence structure (influences of 
longitudinal measurements) (see Chapter 5.5.2). 
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During local observation, turbulence intensity influences in 
turbulent boundary layers can thus only be recorded in a less than 
desireable manner. Therefore, a complete analysis of turbulent 
boundary layers requires a technique which is directly related to 
boundary layer characteristics. For example, such action could 
compare boundary layers with identical impulse loss thicknesses with 
and without free stream turbulence. Such measurements are made in 
Chapter 5.5.2. 

The comparisons of measured heat transfer coefficients with laminar 

and turbulent correlations according to Kays C643 for flows with 

constant material values (T /T = 1.0) (Illustrations 5.15 - 5.17) 

w o» 

provide a concept for the influence of wall cooling on local heat 

transfer. A distinct increase of the Stanton numbers with 

increasing wall cooling can be observed in laminar flow areas. With 

a minor temperature decrease at the wall (T /T = 0.8; Illustration 

w w 

5.15), the dotted curve through the data points of the laminar 
boundary layer confirms a heat transfer which is approximately 7.57. 
higher than the correlation. A further decrease of the temperature 

conditions to T /T = 0.53 (Illustration 5.17) results in an 

w « 

increase of this amount to almost 207 due to increased material 
value influences. So far, in theoretical studies the significance 
of wall cooling for the laminar heat transfer was considered minor 
or negligible (see C64, 65, 673). However, the increased rates 

observed in these new tests show very distinct efffects. 

In contrast, the temperature dependencies of the local turbulent 

heat transfer are visibly reduced. In the case of highest free 

stream turbulence (Tu = 8.77.) the increase of the Stanton 

o 

coefficient is only about 6% during the transition from T /T = 0.8 

w ** 

to T /T * 0.53. 
w ot 
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These heat transfer tests, which were conducted under ideal (78 
conditions with a smooth, unaccelerated boundary layer, provided 
quantitative information and pointed out significant influences of 
free stream turbulence and wall cooling on boundary layer 

development and heat transfer. However, the study of heat transfer 
coefficients, which are computed with reference to the boundary 
layer’s running length, is not sufficient to demonstrate the total 
importance of these influence parameters. As pointed out at several 
locations in this chapter, more in-depth boundary layer analyses and 
data, which pertain directly to the boundary layers, are necessary. 

Because of this, the previously determined heat transfer 

distributions for high and low free stream turbulence intensities 
were supplemented with boundary layer tests for all temperature 
conditions. 

Illustrations 5.18 - 5.23 introduce the boundary layer 

characteristics obtained in these tests. In order to enable re- 

calculations of the test results, a plotting method was selected 
which refers to the boundary layer’s running length, as was 
accomplished with reference to heat transfer distributions. 
Illustrations 5.18 and 5.19 supplement the heat transfer 

measurements when T /T = 0.8 (flow conditions 2 and 2.1 of Table 

w « 

5.1). It must be considered that, with maximum free stream 

turbulence (Tu = 8.7%), the free stream velocity is somewhat lower 
o ' ' 

than with tests where Tu = 1.6, 3.8 and 6.5% (41 m/s versus 47 

o 

m/s). Illustrations 5.20 through 5.23 characterize the respective 
boundary layer data with intensive wall cooling with T /T = 0.64 
and 0.53. 


92 



The applied tangential stress coefficients c^., the kinematic impulse 

loss thickness matrices , and the distribution of the form 

parameter Hj 2 <x) show identical characteristics to those in uncooled 

boundary layers (see Chapter 5.3) in all temperature conditions. 

For example, with increasing outer turbulence , impulse loss 

thicknesses increase due to the advanced location of the boundary 

layer reversal; at the same time, the kinematic form parameters 

decrease. In the transition from Tu = 6.5% to Tu = 8.7%, the 

o o ' 

differences are only minor because, in both cases, the reversal 
occurs in the immediate vicinity of the plate’s leading edge. In 
comparison to the indication of an increase of local turbulent heat 
transfer (Illustration 5.15 - 5.17), there are, with high turbulence 
intensities, relatively minor turbulence-dependent changes of the 
local turbulent tangential stress coefficients c^ (see specifically 
Illustrations 5.20 and 5.22). This corresponding behavior of slow 
speed and temperature boundary layers can be attributed to the 
Reynolds analogy, which applies to unaccelerated boundary layers. 

The enthalpy loss thickness 6^ (see Illustration 5. 19d, 5.21d, (79 

5.23d) is an additional boundary layer quantity which characterizes 
the temperature boundary layer. Its significance is on the same 
level as the impulse loss thickness 6^. Because of the analogy, 
which applies to impulse exchange and heat exchange (Reynolds 
analogy) in unaccelerated flow, the development of both quantities 
along the plate is identical. 

Before I make a general analysis of the influencing parameters of 
boundary layer and heat transfer data introduced in this chapter, a 
more accurate description of the quality of these values is 
required. An expert application of this data for the examination 
and further development of two-dimensional calculation procedures 
requires proof that maximum two-dimensionality existed in individual 
boundary layer flows. Because of this, data obtained with heat 

transfer and boundary layer tests are subjected to a mutual 

examination in Illustration 5.24 . The results of two-dimensional 
enthalpy balance calculations, which are summarized in Illustration 
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5.24, compare the increase of the enthalpy loss thickness between 

two measuring surfaces to the amount of heat transferred to the 

plate surface by cooling. With a maximum deviation of 8-10%, the 

average values fluctuate around the ideal value of 1. If 

consideration is given to the sensitivity of such integral balances, 

the results can be regarded as an excellent confirmation of the 

results achieved with heat transfer and boundary layer measurements. 

Additionally, these results can be considered as a guarantee of a 

two-dimensional boundary layer. An exception to this is the case 

with the lowest free stream turbulence, which shows average 

deviations of up to 30%. The reason for this is the three- 

dimensional (no longer two-dimensional) formation of the boundary 

layer reversal which, in this case, occurs throughout a wide range 

of plate sections. The non-ideal balance values for Tu = 1.6% thus 

o 

emphasizes again the three-dimensional character of the reversal 
process, as was shown in Illustration 5.13 (see Chapter 5.3) with 
the aid of traverse measurements. 

In Illustration 5.24c (T^/T^ = 0.53), the outstanding average values 

for Tu ■ 1.6% with x = 0.21m can also be used as a confirmation of 
o 

the laminar measurement result (specifically the increase of 
laminar heat transfer caused by wall cooling). With this flow, the 
respective reversal occurs after a greater length of plate running 
length. 
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5.5 Analysis of Influencing Parameters 
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The -following chapters provide a detailed evaluation of the 
measurements. An additional component of this evaluation is the 
measurement results with pressure gradients. 

A fundamental description of the influences of free stream 
turbulence, wall cooling and pressure gradients is provided and a 
quantative analysis is made. This includes a mutual super imposition 
of the parameters. The experimental results will be separated from 
their respective situations and will be generalized with dimensional 
analyses. The main focus of my analysis is the influence of heat 
transfer, wall tangential stress, boundary layer profile and 
reversal process. 


5.5.1 Material Value Influence due to Wall Cooling 

5.5. 1.1 Laminar Boundary Lavers 


The laminar heat transfer: 

The heat transfer tests described in chapter 5.4 clearly showed that 

increased wall cooling results in a distinct increase of the laminar 

heat transfer. This was visually demonstrated by a comparison of 

the resulting Stanton coefficients with a relationship proposed by 

Kays C643 for the laminar heat transfer in flows with constant 

material values (Illustrations 5.15 - 5.17). The following 

paragraphs will make a quantitative evaluation of this observation. 

The laminar curves (see Illustrations 5.15 - 5.17), which are 

assigned to the measurements for all temperature ratios T /T are 

w w 

used as a basis. These laminar curves were determined in 

conjunction with the relationships by Kays C643: 


St = A-Re 


-0,5. pr -2/3 


- (|)°’ 75 )- 1/3 


(5.2) 


99 


£ refers to the uncooled take-off distance of the boundary layer. 
The factor "A" was determined with lowest free stream turbulence 


1 . 6 %). 


(Tu = 
o 

which is dependent on the temperature ratio T /T 

w ™ 


It was then summarized to a general relationship, 

and which 

W w 

describes the influence of wall coolings 


A = 0,342-CyTJ 


-0,25 


(5.3) 


Illustration 5.25 shows the temperature dependence of the laminar 

heat transfer, which is given with equation 5.3, in a standard 

application as a factor of the temperature ratio T /T . The laminar 

w w 

measurement values with their respective fluctuations are 
illustrated as confirmation of these results. 
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1 - calculation 

Illustration 5.25s Influence of Wall Cooling on Local Laminar Heat 

Transfer (St Q describes the isothermal, critical 
case in accordance with equation 5.2) 
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According to the test results, the local heat transfer increases up 

to 187. as a result of the decrease of the temperature ratio T /T = 

w 

0.5. In contrast, available information from literature C64, 65, 

67] shows values which are considerably lower. According to the 
theoretical boundary studies by Back et al C673, cooling influences 
can be expected (Illustration 5.25) which, with identical 
temmperature decrease at the wall, are approximately 107. lower. The 
increases of laminar heat transfer, which were calculated in 
accordance with the reference temperature schematic by Eckert (see 
[67, 643) with incompressible laminar heat transfer correlations 

(Illustrations 5.15 - 5.17), are also lower by the same 

percentages. (Eckert proposes the average value between free stream 
and wall temperatures as a reference temperature) The results 
according to the reference temperature procedure - as the 

calculations by Back et al - can be described by the potency law 

<T /T w ) with sufficient accuracy. According to new tests, 

” -'0 25 

however, there is a higher proportional dependency - (T /T ) 

W w 

on the wall cooling (Illustration 5.25). The biggest contrast to my 
tests are the proposals by Kays C643, which are based on estimated 
material value analyses, where cooling influences in boundary layers 
with constant free stream velocity can be completely neglected. To 
date, the low or negligible effects of the cooling were based on the 
following: the effect of the boundary layer profiles, which became (86 

fuller as a result of cooling, and the steeper wall gradients are 
almost completely compensated by the decreased heat conductivity and 
viscosity near the wall (see Back et al C673). However, recent 
tests have shown that there is an imbalance of the profile 
deformation, which results in distinctly higher increases of the 
laminar heat transfer as a function of cooling than was presumed in 
previous theoretical analyses. 
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In order to clarify the pro-file deformations on which the resulting 
heat transfer increases were based, Illustration 5.26a shows laminar 
velocity profiles with an identical running length Reynolds number 
Re for different temperature conditions. In this application, 
which was made with the similarity coordinates of laminar balancing 
boundary layers, fuller velocity profiles can be observed with 
increasing wall cooling. These profiles are different from the 
ideal solution of isothermal boundary layers according to Blasius 
(see C109]). As a confirmation of this procedures, they experience 
a decrease of the kinematic form parameter from H.^ = 2.55 (T /T = 

Ljl W W 

1.0) to Hj'ij = 2.39 = 0.53). The fuller profiles show steeper 
wall gradients which, despite reduced viscosity and heat 
conductivity, result in a measured increased in the laminar heat 
transfer at the wall. 

Suitable transformation equations provide the opportunity for 
applying incompressible laminar profile equations to flows with heat 
transfer. An example of this is shown in Illustration 5.26b . The 
Levy transformation [116] shown in Illustration 5.26b supplements 
the standard wall distance coordinate f> with an additional material 
value weighting, which is directly dependent on the wall distance y: 

-i y 

7 ‘ f dy 

x J 

o 

Because of this transformation, the profiles which were changed by 
cooling are converted to the standard condition for adiabatic 
boundary layer flows, according to Blasius. With the Eckert 
transf ormation, which was also studied for this profile evaluation, 
the influence of wall cooling was compensated considerably less than 
with the Levy procedure (Illustration 5.26b). This was also 
confirmed by the inaccuracy of the Eckert transformation, as is 
shown in Illustration 5.25, where dependency of the laminar heat 
transfer on the wall cooling was not adequately considered. 
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5.5. 1.2 Turbulent Boundary Lavers 


Heat transfer and wall tangential stress: 

Events observed in turbulent boundary layers with respect to running 
length are always subjected to additional influences, which are 
caused by the formation of the respective boundary layer reversal. 
In order to eliminate resulting side effects, Illustrations 5.27a 
and 5.27b provide an analysis of the turbulent friction coefficients 
measured by me for the temperature ratios T /T = 1.0; 0.8; 0.64 and 

W W 

0.53 as a factor of a similar value, which has a direct relationship 
to the boundary layer: the impulse loss thickness Reynolds number 
Re^. Both illustrations describe cases with constant free stream 
turbulence <Tu = 6.5% and 8.7%). With increased wall cooling, the 
friction coefficient increases by up to 16% (Illustration 5.27). 
The increased rates in both cases of varying and, at the same time, 
high free stream turbulence are almost identical. Also, they are in 
good agreement with the results by Sill C783, which were determined 
in previous tests for turbulent boundary layers with very low outer 
turbulence. The increases of the turbulent heat transfer as a 
result of wall cooling are thus independent of the amount of free 
stream turbulence. 


In order to quantify the influence of wall cooling, comprehensive 
analytical studies were made. Illustration 5.17 shows the result of 
these analyses with the aid of curve systems which, for given 
temperature conditions, describe the behavior of tangential stress 
coefficients in turbulent boundary layers with heat transfer. The 
relationship is as follows: 


= A- Re 


- 0,25 


2 ' 0,58 

wv * 


Re 2 


u -6 2 


(5.4) 


This equation is based on the statement for the friction coefficient 
in accordance with equation (5.1), as it applies to the flow on 
adiabatic walls. By introducing a suitable reference temperature 


(mean value between wall and free stream temperatures) to determine 
the material values contained therein, this statement was expanded 
to flows with variable density. The connection described by 
equation (5.4) with the values Re 0 and c^. which, in turn, relate to 
the free stream material values, was then achieved with a formula 
substitution proposed by Back et al C763 of these newly introduced 
reference material values. 
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1 - grid 


Illustrations 5.27a and b: Turbulent Friction Coefficients 

Linder the Influence of Wall Cooling 
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To confirm the theoretical temperature influences of equation 5.4, 
the empirical constant "A" was determined by approximating this 
relationship to the respective measurement points of the uncooled 
boundary layer flow. With this determination for the constant "A", 
the relationship 5.4 describes the measured influence of wall 
cooling in both cases of varying free stream turbulence with 
sufficient accuracy for all temperature conditions (see Illustration 
5.27 a and b). Compared to the test, the calculation shows slightly 
increased tangential stress coefficients. However, the difference 
is less than ZV. in all cases. Equation 5.4 thus provides a 
practical condition which describes the basic cooling influence in 
turbulent boundry layers in a simple manner. 

Because comparative isothermal tests (T^/T^ = 1.0) were not possible 
for the heat transfer, the above result for tangential stress 
coefficients was linked to measured analogy factors in order to 
determine the coaling influence on the heat transfer. The analogy 
factors determined in my experiments (see Illustration 5.43, Chapter 
5.5. 2. 2) 2 x St/c^ fluctuate around the mean value 1.1. They 
neither show an identifiable dependency on the wall cooling nor a 
direct dependency on the free stream turbulence. This proves that 
the increases, which are a factor of cooling and which exist in the 
heat transfer, are identical to those in the wall tangential 
stress. Therefore, the temperature dependency, which was obtained 
in equation 5.4 for the friction coefficient, is also suitable to 
describe the increase rates in heat transfer. 

Velocity and Temperature Profiles: 

The formation of heat transfer and wall tangential stress in flows 
with varying density is mainly determined by the properties and the 
behavior of the boundary layer profiles. Illustration 5.28 shows 
typical velocity and temperature profiles whose type and development 
were measured along the plate surface in all temperature conditions 
and turbulence intensities. The measurement in Illustration 5.18 


refers to turbulent boundary layers with high free stream turbulence 
(Tu q * 6.5V.) and the highest wall cooling = 0.53). Measuring 

surface 1 is located in the immediate vicinity of the end of the 
laminar-turbulent boundary layer reversal. As can be expected with 
unaccelerated boundary layers, there are almost identical velocity 
and temperature profile forms. This can be attributed to a 
synchronous development of slow speed and temperature boundary 
layers. With an increasing Reynolds number Re£, Reynolds number 
influences become less significant and the profiles visibly approach 
a form which is characteristic for equilibrium boundary layers. Due 
to Prandtl numbers for gases (in our case: air, Pr = 0.71), which 
deviate from 1, the temperature boundary layer thickness is 
always greater than the flow boundary layer thickness 6, and the 
respective value of the temperature profile - with the same wall 
distance - is always lesser than the value of the velocity profile 
(see Illustration 6.14b; Chapter 6.7. 1.3; and Riid et al [1233). 
Drastic difference in the development of velocity and temperature 
profiles will only occur if the flow has considerable pressure 
gradients, or if the heat stream paramater 

«i l ' 

n c,/2 ’ • * dx (5.5) 

% 

which is analogous to the Clauser parameter, has very high values. 
According to studies by Taylor [1213 and Sill C783, pressure 
gradients and distinct heat stream paramaters A have a considerable 
additional influence on the temperature profile. 
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1 - surface 

Illustration 5.28: Velocity and Temperature Profiles in a Highly 

Cooled Turbulent Boundary Layer 

The Logarithmic Wall Theory: 

The boundary layer profiles in Illustration 5.28 have the same basic 
matrices as those measured in adiabatic boundary layers with 

identical free stream turbulence but without wall cooling (see 
Illustration 5.6). Therefore, changes which affect coaled boundary 
layer profiles cannot be recognized with this application. A 
demonstration of the influences caused by wall cooling requires a 
similarity analysis of the profiles in the form of the logarithmic 
wall theory. Because of this, Illustration 5.29 shows dimensionless 
velocity and temperature profiles with identical free stream 

turbulence (Tu q = 8.7'/.) for varying cooling intensities. The 

profiles were measured in surface 4. With falling temperature 
conditions, the applicable Reynolds number Re^ has the values: 



(T - TwMP • Cp)wUx 
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Re 0 - 1,640.1 1,550, 1,530 and 940. As shown in Illustration 5.29, 
the profile properties which exist when wall cooling is present can 
now be clearly seen. Even with high wall cooling, the measured 
profiles show a distinct logarithmic wall section. However, with 
increased wall cooling, the matrices show higher velocities and 
temperatures plus a slight profile twist. 



Illustration 5.29s Dimensionless Velocity and Temperature 

Profiles 
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Only a very -few experimental studies (Sill C78] and Back et al C761) (93 

have been conducted on the effects of highly changeable material 
values in turbulent sub-sonic boundary layers. These tests showed 
an increase of velocity and temperature values with semi -logarithmic 
application. However, the simultaneous rotation of the logarithmic 
wall theory was, in most cases, not demonstrated. For an evaluation 
of our test results I will therefore use the theoretical study by 
Rotta C853 on the influence of Mach number and heat transfer in 
turbulent flows. By integrating the boundary layer differential 
equations, Rotta C853 deduced an asymptopic wall theory which 
addresses the influence of wall cooling. Illustration 5.29 compares 
this wall theory to the profile measurements. Mach number influences 
are neglected, as they are not important in this case. If the 
sensitivity of the semi -logarithmic application method is 
considered, there is a very good agreement of the theory with the 
test results. The velocity and temperature value increases as well 
as the rotation of the wall theory are correctly recorded by Rotta’ s 
equations. This also confirmed Rotta’s assumption that the cooling 
effects will not change the basic turbulence structure but are 
generated exclusively by the material values, which vary throughout 
the boundary layer. 

When revising calculation models, it is often important to have a 
simple and quick test method for computing profile matrices. To 
date, no suitable methods have been developed; therefore, an attempt 
was made to apply, based on the reference temperature principle, 
the simple logarithmic wall theory for adiabatic flows in highly 
cooled boundary layers. At the same time, it was to be determined 
whether reference temperature methods used in integral procedures 
are suitable to expand the profile matrices, to a situation with 
wall cooling (see Chapter 6; development of the integral 

procedure) . The mean temperature between free flow and wall, which 
was used in the tangential stress analysis, was again considered 
valid. 


no 



Illustration 5.30a shows, -for the representation of two cases with 
high wall cooling, the logarithmic velocity and temperature profiles 


with the coordinates 


y Ref ’ 


Ref 


and T_ ,, which were formed with 
Ref’ 


material values at reference temperature. The wall theories for 
adiabatic boundary layers are also illustrated for a direct 
evaluation with of the useability of this reference temperature. As 
expected, the material value transf ormation near the wall again 
shows agreement of the measurement points with the logarithmic wall 
theory for adiabatic boundary layers (see Illustration 5.30a). 
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Based on this determination, the logarithmic wall theory for 
adiabatic flows can also be applied to boundary layers with high 
temperature variations, if the coordinates T + , u + and y + are 
generated from the material values (density, viscosity) at reference 
temperature. If, with the aid of suitable material value 
relationships for density and viscosity (p s 1/T; v = >*/p ; H s 
T^'^), we perform a simple re-transf ormation of these coordinates 
into quantities formed with wall material values, the wall theory 
for adiabatic boundary layers is validated in a general context. 
This context has a common notation, which is independent of wall 
coordinates, and, at the same time, describes the influences of wall 
cooling. The wall theory for velocity and temperature profiles is 
as follows: 


u + = 1_ . m / + c * 

K 


T+ ■ • 1n > + + S* 


with 


* l °* 5 

‘ = K /[± (1 + T,/T w )] 

c * = [i (1 ♦ yrj] 0,5 (M • In (T-rf^) + c > 


(5.6) 


(5.7a) 


(5.7b) 
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The relationships of the logarithmic wall theory with wall cooling 
gained in this manner are considerably simpler and clearer in 
comparison to the Rotta equations <see £853). Their validity is 
emphasized in Illustration 5. 30b by a good agreement with the 
applicable profile test data. The deduced relationships are a good 
reproduction of the profile rotation and the increased temperature 
and velocity values with increasing wall cooling. 

With the newly developed wall theories there is a similarly good 
agreement as in Illustration 5.30b of all other test situations with 
wall cooling. This is confirmed in Illustration 5.31 in a 
comparison of all measured profile constants c land c£ \ which were 
calculated in accordance with equation 5.7. The measurement values, 
which are a factor of the temperature relationship, were determined 
by adjusting the wall theory in accordance with equation 5.6 to the 
measured profiles. and *? were assumed as a parameter in 
accordance with equation 5.7. Additionally, Illustration 5.31 shows 
that, based on the profile rotation due to wall cooling, there are 
only minor changes in the profile constants ° k and c* . In addition 
to the simple method of calculating the wall theory in accordance 
with equation 5.6, this provides another method of simplification. 
In accordance with Illustration 5.31, the values c and c T with 
temperature conditions in excess of T T = 0.5, can also be 

W « 

satisf actori ly determined with their adiabatic values. 
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1 - equation (5.7); 2 - test values 

Illustration 5.31; The Pro-file Constants of the Logarithmic Wall 

Theory in Accordance with Equation 5.7 and 
Test Values 

5.5. 1,3 Integral Boundary Laver Characteristic Values for (97 

Wall Cooling 

A favorable comparison of kinematic boundary layer characteristic 
values to material values in laminar and turbulent boundary layers 
in Illustration 5.32 is the preliminary conclusion of my 
comprehensive and detailed analysis of the influences of wall 
cooling. Material value influence occurs in all characteristic 
values. When wall coolings are present, they effect an increase of 
the impulse loss thickness. At the same time, they cause an even 
greater decrease of the displacement thickness, which results in a 
decrease of the form parameters in laminar and turbulent 

boundary layers. In our temperature range, impulse loss thickness 
ratios increase by up to 207.. Despite very different profile forms, 
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there are no striking differences between laminar and turbulent 
boundary layers. In contrast, the differences in the displacement 
thickness and, consequently, in the form parameters are more 
distinct. 


citgimal PASS -IS 
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1 - turbulent boundary layer 

2 - laminar boundary layer 

Illustration 5.52; Laminar and Turbulent Boundary Layer 

Characteristic Values under the Influence 
of Wall Cooling (Test Values) 
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5.5.2 Influences of Free Stream Turbulence 
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This chapter will provide a quantitative analysis of the influences 
of free stream turbulence. A comparison of test data with and 
without wall cooling will determine whether turbulence intensity 
effects are subjected to additional wall cooling influences. 

5,5. 2.1 The Laminar-Turbulent Boundary Laver Reversal 

In the design of cooled gas turbine blades, the question is always 
raised whether the reversal processes caused by free stream 
turbulence are also subjected to dramatic influences from wall 
cooling. The linear stability theory C223 for boundary layers 
without external turbulence lets us presume that this is the case. 

The same applies to the measured laminar profiles in Illustration 
5.26, because the profile matrices, which have become fuller as a 
function of cooling, point to a stability increase of the boundary 
layer. In a careful comparison of reversal events observed in 
experiments with data of other authors, which have all been 
determined in adiabatic flows, I would like to identify the actual 
circumstances. 

An evaluation of the laminar-turbulent reversal is made with the aid 
of the measured heat transfer distributions along the plate. The 
reversal Reynolds numbers Re^g, which were deduced from heat 

transfer tests, resulted from the intersection of laminar 

compensating curves with a straight line, which is adjusted to the 
matrix of steeply increasing Stanton numbers in the reversal area 
(see Illustrations 5.15 - 5.17). Reversal locations and matrices, 
which are comparable to information available in literature, are 
thus determined. These transform the three-dimensional process into 
a good approximation of a two-dimensional description which can 

identify the reversal initiation point. 
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The reversal determination in previous studies (see C45]) was 
obtained with the aid of probes. The reversal initiation in my 
analysis of the heat transfer is in context with a noteable 
deformation of the boundary layer profiles. As demonstrated by Owen 
[253 in his own boundary layer tests and with the aid of those 
conducted by Schubauer and Klebartoff [17], a distinct boundary layer 
turbulence formation is required (see Chapter 2.2). Because of 
this, the visible initiation of the reversal in the heat transfer is 
not equal to the formation of boundary layer turbulence. 



1 - reversal completion 

2 - reversal initiation 


Illustration 5.33: Running Length Reynolds Numbers Re at 

M 

Reversal Initiation and Completion 
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Illustration 5.33 shows the running length Reynolds numbers at 
reversal initiation and at the end of the reversal (see Table 5.2), 
which were determined with heat transfer tests, as a factor of free 
stream turbulence. It also makes a comparison with results achieved 
by Abu Ghannam and Shaw 145], Schubauer and Skramstad 1143, Hi si op 
11183 and Blair 1363. The Reynolds numbers refer to the viscosity 
of the free flow. In accordance with a proposal by Abu Ghannam and 
Shaw, average turbulence intensities are used. These are formed 
with the local free stream turbulence Tu at the reversal location 
and with the respective inlet turbulence Tu q at the beginning of the 
plate. This takes into consideration the turbulence intensity 
decrease along the plate and, therefore, the history of the boundary 
layer up to the point of reversal. 

My own determinations of the reversal initiation are in good 
agreement with the data of other authors. Even the fluctuations 
with varying wall cooling can be integrated in other test results, 
which were obtained exclusively under isothermal conditions. This 
shows that, contrary to information obtained from the stability, 
theory, no influence of wall cooling is noted on the reversal 
initiation. The only differences to previous tests are noted in the 
description of the reversal completion. With low turbulence 
intensities, my tests show a tendency toward higher Reynolds numbers 
and, consequently, toward greater reversal lengths while, with high 
turbulence intensities in excess of 5%, the reversal is completed 
after shorter running lengths. 

In order to present a more general picture not dependent on the 
running length, as is required for the development of general 
reversal criteria, Illustration 5.54 provides a reversal analysis 
with the aid of local boundary layer characteristics. Instead of 
running length Reynolds numbers, impulse loss thickness Reynolds 
numbers are studied at the initiation and completion of the 
reversal. The required impulse loss thickness Reynolds numbers were 
determined with the running length Reynolds numbers at reversal 
initiation described above. 


(100 


CiOTiA!. &&S 

OF POOR <jUAUTY 



Illustration 5.34: Impulse Loss Thickness Reynolds Numbers 

at Reversal Initiation and Completion 


As shown by Illustration 5.55 with test values for lowest free 
stream turbulence, the impulse loss thickness Reynolds number at 
reversal initiation can, for all temperature conditions, be 
calculated with the simple Blasius relationship 

Re 2 = 0,664 /Re~ (5.8) 

with the aid of the applicable running length Reynolds number. 
Because, for higher turbulence intensities, laminar boundary layers 
could not be accessed with probes, the Blasius relationship 
(equation 5.8) was also used for these cases to determine the 
impulse loss thickness Reynolds number at reversal initiation. This 
procedure, which is also used by other authors, is supported by the 
experience (see Illustrations 5.16 and 5.17) that laminar boundary 
layers in unaccelerated flows are not subject to turbulence 
intensity influences. The impulse loss thickness Reynolds number at 
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the end of the reversal was determined only for flow conditions 
obtained with boundary layer measurements. To do this, the measured 
Re 2 matrices were either interpolated or extrapolated. 



1 - Re _ from WLi tests; 2 - grid 

A u 

Illustration 5.35; Determining the Reynolds Number Re 2 

at Reversal Initiation 

The application of the reversal Reynolds numbers Re_„, Re„_ in 
Illustration 5.35 again compares my tests to the previous results 
achieved by Hi si op Cl 173, Brown and Burton Cl 183, Martin Cl 193, Well 
C1203, and Abu Ghannam and Shaw C453. The distinct spread of 
available measurement data in this diagram indicates the difficulty 
of an exact determination of the reversal. As in Illustration 5.33, 
my analysis of the impulse loss thickness Reynolds number is also in 
good agreement with the data of other authors determined exclusively 
under isothermal conditions. Because the test results show no 
tendencies which could be attributed to the temperature ratio T /T , 
they confirm that no boundary layer stabilization exists, which 
could be attributed to cooling influences. 
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According to both stability theory and previous tests, there should - 

•for turbulence intensities in excess of 5*/. - be a minimum 

asymptopic value for the reversal Reynolds number (with ~ *60) 

which should not be exceeded, even with maximum free stream 

turbulence intensities (see Illustration 5.33). My tests do not 

confirm this reversal behavior. The heat transfer distributions 

along the measuring plate (see Chapter 5.4) have shown that the 

reversal with high turbulence intensities (Tu = 6.5 and 8.7*/.) 

o 

occurs very early; it is almost complete at the first measuring 
location <x = 47 mm). In accordance with calculations from equation 
5.8 it can be expected that, upstream in this initial measuring 
surface, there are impulse loss thickness Reynolds numbers in the 
reversal initiation which are considerably lower that the asymptopic 
boundary value (with Re 2 g 15 160). Tests with free stream 

acceleration come to the same conclusion. With high turbulence 
intensities, the reversal always initiates directly at the leading 
edge of the plate (see Chapter 5.6). 

The reasons for the discrepancy between the new and the few previous 
test results, which confirm the presence of a stability limit, can 
be found in the varying geometric parameters of boundary layer 
generation (design of the leading edge of the plate). In distinct 
contrast to previous tests, our measurements were taken with a very 
thin, streamlined leading edge of the plate (with a thickness of 
only 2 mm). With thick leading edges of a plate and a respective 
stream-lining, acceleration effects occur which provide additional 
stability to boundary layer development, thus displacing the 

reversal to higher Reynolds numbers. The boundary layers generated 
on thin leading edges are not subject to additional 

stabilization. They can therefore react much more sensitively to a 
fluctuation in the flow caused by free stream turbulence. They may, 
as with our boundary layer flows, experience a premature reversal , 
especially when high turbulence intensities are present. It cannot 
be discounted that there may be a situation where the increased 
sensitivity to free stream turbulence completely prevents the 
formation of an initial laminar boundary layer. As a consequence, a 
turbulent boundary layer may be generated at boundary layer 

initiation. 
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These geometric aspects of reversal behavior and the resulting 
consequences -for heat transfer are especially relevant -for modern 
engine design: Aeronautic engines have, in their initial turbine 
stages, very thin blade profiles which are constantly exposed to 
highly turbulent gas circulation. 

According to recent determinations, there is, when free stream 
turbulence is present, no stabilizing influence of wall cooling on 
the boundary layer which might result in higher reversal Reynolds 
numbers. With the turbulence and cooling intensities, which have 
technical relevance in our case, the influence of free stream 
turbulence dominates over reversal stabilizations which, are 
normally expected when wall cooling is present and which are 
predicted by the linear stability theory for flows without external 
turbulence. This, however, is not a fundamental contradiction to 
stability theory. The existence of cooling effects, which stabilize 
the boundary layer, is also demonstrated in our tests. It is 
clearly described in Illustration 5.26 by the fuller laminar profile 
forms. However, its significance - in contrast to the importance of 
free stream turbulence - is negligable. 

In addition, the analysis of the reversal events measured with wall 
coaling and the comparison with results for adiabatic boundary 
layers (Illustrations 5.33 and 5.34) point out the Reynolds number, 
which is based on the material data of the free flow, as an 
appropriate description measure for the reversal - at least for low 
and medium turbulence intensities. Despite drastic viscosity 
differences between free flow and wall, the introduction of reversal 
Reynolds numbers, which were formed with material values at an 
average reference temperature, has not proven to be practical. This 
is clearly emphasized by the tests with accelerated flow in 
Illustration 5.36 and also applies to highest turbulence 
intensities. By using both pressure gradient contours, slightly 
accelerated flows with varying wall cooling were generated in 
Illustration 5.36a. Highly accelerated flows with varying wall 
cooling were generated in Illustration 5.36b. The plate flow 
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velocities were adjusted to the selected free stream temperature in 
such a manner that, in the test cases of Illustrations 5. 36a and 
5.36b, there are identical plate flow Reynolds numbers Re^ . By 
maintaining constant plate running length Reynolds numbers, which 
were formed with the free stream data at the influx, and with 
varying wall cooling, constant matrices of the acceleration 
parameter k(x) were achieved (see equation 3.1) along the plate, 
i.e. events which are comparable with respect to flow (see 
Illustrations 5.45, 6.15, 6.19). At the same time and in spite of 

varying free stream temperature, identical local running length 
Reynolds numbers occur at each location x, if they are also formed 
with free stream viscosity. Even with maximum turbulence 
intensities, the Stanton numbers (Illustration 5.36) measured under 
these conditions show no influence of the wall cooling on the 
boundary layer development or on the reversal. There are no 
stabilizing effects on the reversal due to wall cooling. Reversal 
events show a distinct correlation to the Reynolds numbers which 
refer to free stream material values. As shown in Illustration 
5.36, the influence of the cooling is limited to the fact that the 
local heat transfer rates can be increased with a de creased 
temperature ratio. Increase rates are comparable to test results 
achieved without pressure gradients (see Chapter 5.4). 















5. 5.2.2 Turbulent Boundary Lavers 
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Chapter 5.3.2 showed the effects of increased free stream turbulence 
on adiabatic turbulent boundary layers. Based on this, the 

following chapter will analyze events with wall cooling as a factor 
of turbulence intensity. Turbulence intensity influences, which are 
almost completely known for other situations, have not been 
researched yet for this combination. The complexity of events 
induced by free stream turbulence, which was shown in previous 
studies, requires that accompanying longitudinal measurement and 
Reynolds number influences must be included in this study. In 
addition to heat transfer and resistance coefficients, the behavior 
of the boundary layer profiles must also be considered. 

Velocity and Temperature Profiles: 

Turbulent boundary layers with identical Reynolds numbers Re^ and a 

different free stream turbulence are shown in Illustration 5.37 for 

the case of moderate wall coolinq (T /T_. = 0.8) and in Illustration 

w o» 

5.58 for the case of maximum wall cooling (T /T = 0.53). The 
logarithmic wall lines were determined by adjusting the temperature 
correcting wall theory, developed in the previous chapter, to 
equation 5.6. As shown by the profile matrices, which were measured 
for wall cooling for the first time, the influence of increased free 
stream turbulence is in complete agreement with the events in 
adiabatic flows. It is, in turn, limited to the decrease of the 
wake portion with increasing free stream turbulence intensity. In 
agreement with the measurements by Blair C36, 393 in a slightly 
heated flow with an almost constant density, this also applies to 
the temperature profile with high wall cooling (Illustrations 5.37b 
and 5.38b) . 

In order to make a quantitative evaluation of the illustrated 
profile deformations and a possible discovery of cooling influences 
not shown by Illustrations 5.37 and 5.38, a comprehensive similarity 
analysis for the kinematic form parameter is conducted in 

Illustration 5.39 . With its respective values for turbulence-f ree 
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The -form paramaters -for boundary layers of negligible free stream 
turbulence, which were not available in my tests, were taken from 
the data by Murlis C333. Illustration 5.39 confirms a distinct 
decrease of the form parameter with increasing free stream 
turbulence, which - for the most part - is in agreement with the 
matrix determined for isothermal flows by Green C51]. The decrease 
of the form parameter can be attributed to fuller velocity 
profiles. Starting with approximately 4%, the form parameters in 
Illustration 5.39 - similiar to the behavior described by Robertson 
and Holt C503 - strive toward a boundary value. When this 
phenomenon occurs, the wake portion has achieved its final condition 
which is marked by complete disappearance (see Illustrations 5.37a 
and 5.38a). Because the tests show no systematic differences for 
cases with and without wall cooling, wall cooling generally has no 
additional influence on the transformation of the boundary layer 
profiles caused by turbulence intensity in the free flow. 



Illustration 5.39; Form Paramater H~as a Function of Local 

Free Stream Turbulence 
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Tangential Stress Coefficients and Reynolds Analogy: 

As already known, free stream turbulence increases the impulse 
exchange within the boundary layer. This results in higher velocity 
gradients at the wall and, therefore, higher wall tangential 
stresses. A quantitative illustration of the measured tangential 
stress increase is contained in Illustration 5.40a and Illustration 
5.40b. For a basic explanation of these circumstances, these 
illustrations compare the resistance coefficients of boundary layers 
with identical impulse loss thickness Reynolds numbers with and 
without external turbulence. The values for turbulence-free 
external flow, which are required for this comparison, were 
calculated on the basis of an analytical relationship (equation 5.1 
or 5.9) for incompressible boundary layers. 

c f = 0,024 • Re 2 ' 0,25 (5.9) 
o 

Equation (5.9) was used because of its excellent agreement with the 
test results for turbulence-free external flows by Murlis, Wieghard 
(see C333) and Coles C 1 153 (max. deviation 27.). 


The influence of the wall cooling, which increases the tangential 
stress, on the resistance coefficient in turbulence-free flows was 
considered with the temperature correction function achieved in 
Chapter 5.5. 1.2 (equation 5.4): 




0,58 


(5.10) 


Illustration 5.40a shows the resulting values c^/c^ for all local 
turbulence intensities and temperature ratios. The illustration 
shows that an increase of the tangential stress coefficient, which 
is a function of the turbulence intensity, exists. This increase 
can satisfactorily be illustrated with a relationship by Simonich 
and Bradshaw C543, which can be applied to our case of flows with 
high wall cooling. The relationship states that, when the free 
stream turbulence is increased by 1%, the wall tangential stress 
increases by approximately 27.. This statement plus the fact that no 
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distinct additional influences caused by wall cooling can be 
observed (see Illustration 5.40) is of great significance to the 
development of boundary layer calculation procedures and to the 
practical application for gas turbine blades. 

Because our tests were based on turbine-relevant boundary layer 
flows with low impulse loss thickness Reynolds numbers of R g = 330 - 
2,000 (see Illustration 5.27), there are Reynolds number influences, 
which superimpose themselves to the effect of external turbulence 
with an opposing sign. In order to demonstrate this, an application 
procedure was selected in Illustration 5.40b, as proposed by Blair 
C36D, which considered the Reynolds number as an additional 
parameter. In this application, the measurements show - as a 
consequence of the corrected Reynolds number influences and, 
therefore, a confirmation of the existence of Reynolds number 
influences - lesser fluctuations and a clear dependency on the 
turbulence intensity. The balancing line in the illustration 
describes the combined influence of Reynolds number and free stream 
turbulence as it was obtained with the tests for the given Reynolds 
number and turbulence intensity areas. The reference relationship 
for the balancing curve (Illustration 5.40b) points out that 
Reynolds number effects counter — act the influence of free stream 
turbulence below Re^ = 1,000. 
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Illustration 5. 40; Dependency of Turbulent Tangential Stress 

Coefficients on Local Free Stream Turbulence 
Tu when Considering the Impulse Loss Thickness 
Reynolds Number Re^ <Re^ = 330 - 2,000) 
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In accordance with recent findings, increases of the wall tangential 
stress as a function of turbulence intensity also depend on the eddy 
structure of the free stream turbulence. According to statements by 
Bradshaw and Hancock [55, 593, turbulence influences are highest, 

when the longitudinal measurements 

— - 3/2 

L“ = - — — <5.11) 


which are formed by the longitudinal component of the free stream 
turbulence have the same size as the reference boundary layer 
thickness. As shown in Illustration 5.41 . this applies for a large 
portion of our tests. In addition, however, there are distinctly 
greater longitudinal measurements which, due to the smaller boundary 
layer thicknesses, can always be found in those measuring surfaces 
which are closest to where the boundary layer is formed. 

In order to understand the influence of free stream turbulence 
longitudinal measurements, Hancock [563 (see [55, 591) - in 
comprehensive tests and comparisons with literature - developed a 
correlation which describes the tangential stress increases as a 
function of turbulence intensity and the ratio L e %. It considered 
the fact that, with largelongitudinal measurement ratios 4 ^/i", the 
tangential stress increase is supressed. It further considered that, 
with very low turbulence intensities below 1 X, the tangential stress 
- according to tests by Meier and Kreplin [583 - the tangential 
stress no longer increases in a linear manner but squarely with the 
free stream turbulence. As an expansion of this study, Blair [393 
developed an empirically deduced supplement, with which Hancock's 
correlation can also be applied to flows with additional Reynolds 
number influences. Illustration 5.42 compares this modified 
correlation to my own measurements. The correlation and the 
surrounding hatched range, which encompasses all test data of other 
authors (see [393), represent the paramaters for this comparison. In 
Illustration 5.42, my own test data are in good agreement with the 
total picture of available results; they also confirm Hancock’s 
correlation. 
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Illustration 5.41; Longitudinal Measurements of the Free Stream 

Turbulence 



1 - Hancock’s Correlation Curve 

Illustration 5.42; Tangential Stress Increase with Free Stream 

Turbulence when Considering Reynolds Number 
and Longitudinal Measurement In-fluences 
(see Illustrations 5.40 and 5.41 for Symbols) 


132 




However, the good agreement in this study was achieved by using the 

Reynolds number correction. Without consideration of the correction (112 

term s (see Illustration 5.42), the measuring points had distinctly 

greater abscissa values which made for a bad agreement. This showed 

that, when very low Reynolds numbers (Re£ = 330-2,000) were used, 

the boundary layer related Reynolds number influences outweigh the 

longitudinal measurement influences (L^/S * 1.1 - 3.0) of the free 

e 

stream influence. 

A determination had to be made of how turbulence intensities, 
identified from comprehensive studies of the friction coefficient, 
can be applied to heat transfer. To do this, the connection between 
heat and impulse exchanges was examined with the aid of measured 
analogy factors. The analogy factors 2 x St/c^, which were obtained 
from the measured Stanton numbers and the tangential stress 
coefficients, are shown in Illustration 5,42 . The results evenly 
fluctuate around the mean value 2 x St/c^. = 1.1. The analogy factor 
is independent of external turbulence and is not influenced by wall 
cooling. Consequently, wall tangential stress and heat transfer are 

increased by free stream turbulence. The Reynolds number and 
longitudinal measurement influences observed for tangential stress 
coefficients are also present in the heat transfer. 



Illustration 5.43: Analogy Factors with Varying Wall Cooling and 

Free Stream Turbulence 
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The answer to the original question whether wall cooling alters 
turbulence intensity influences can be summarized as •fallows: No 
additional effects on the wall cooling were observed in analyses on 
the influence of free stream turbulence on the reversal, wall 
tangential; stress and heat transfer. Turbulence intensity effects 
and material value influences caused by wall cooling are thus not 
connected. They can therefore be studied independently or they can 
be combined by a simple super imposition for a total effect. This is 
confirmed by Eriksen et al C79, 803 and Cheng et al C843 by direct 
turbine tests in boundary layers with very low or negligible free 
stream turbulence. These tests showed that, in turbulent boundary 
layers with highly varying material values, the turbulence structure 
is the same as in the isothermal case. Therefore, the influence of 
wall cooling can be regarded separately. 

5.5.3 Influences of Negative Pressure Gradients 

In gas turbines, the influences of external turbulence and wall 
cooling are superimposed by high pressure gradient effects, which 
play an important role in the development of heat transfer and the 
boundary layer along the blade surface. Only very few experimental 
studies on the influence of pressure gradients on transitional 
boundary layers with external turbulence have been conducted. These 
studies were only made for low acceleration intensities (k < 1.0 x 
10 (see reference list, Chapter 2.2). The goal of my studies was 
to analyze flows with free stream acceleration as close to realistic 
turbine conditions as possible. The extraordinarily high increases 
of laminar heat transfer (see e.g. Krisnamoorthy C993) observed in 
blade tests led to the conclusion that there are effects, which are 
yet unclear and which occur, when very highly accelerated boundary 
layers are admitted. 

My analysis is made exclusively as a result of measured heat 
transfer distributions. Due to the very high accelerations, the 
boundary layer estimates made in design calculations resulted in 
very thin boundary layer thicknesses. Therefore, using probes for 
boundary layer measurements was considered impractical. 
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5.5.3. 1 Basic Influences of Free Stream Acceleration 
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In order to be able to demonstrate the basic influences of highly 
negative pressure gradients in transitional boundary layers, I 
conducted comparative tests with and without free stream 
acceleration. The comparison was made with identical inlet 
conditions at the leading edge of the plate and acceleration events 
with varying intensity along the test plate. 

Illustration 5.44 shows situations measured with a constant inlet 

Reynolds number (Re Lo = 6.4 >: 10 5 ) and intensive wall cooling (T w /T M 

= 0.64). The Stanton numbers applied as a factor of the plate’s 

running length were formed with local free stream velocity and the 

local free stream material values. The lower acceleration occurred 

with contour 1 and an acceleration parameter along the plate which 

was almost constant (k = 1.2 x 10 > . The intensive acceleration 

with k _ = 3.2 x 10 was realized by usina contour 2 (Illustration 

max 

5.45) . In contour 2, flow conditions on the vacuum side of gas 
turbine blades are simulated. The contour generates an alternating 
velocity matrix; the acceleration in the forward portion of the 
plate greatly increases and then experiences a slight retardation 
(Illustration 5.45). 

The heat transfer coefficients applied in Illustrations 5.44a and 
5.44b represent the test results with very low free stream 
turbulence (Tu = 1.6% and 2.3'/.). They describe flow events with a 
highly transitional character which, under the influence of negative 
pressure gradients, show varying reversal behavior. As clearly 
shown in Illustration 5.44b, the low acceleration (k = 1.2 x 10 ^) - 
in contrast to the unaccelerated situation (k = 0) - results in an 
extension of the reversal area with an almost unchanged reversal 
initiation point. The higher acceleration, with k = 3.2 x 10 
shows a later reversal initiation with the minimum Stanton number 
further down-stream. The reversal delay and the expansion of the 
reversal area can be attributed to the known effect of negative 
pressure gradients, which stabilizes the boundary layer. This 
results in fuller boundary layer profiles and steeper wall gradients 
(see Chapter 2.2). In accordance with classical reversal theory 
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(stability theory), the higher velocity gradients at the wall are 
mainly responsible -for the greater stability of the laminar boundary 
layer (see Schlichting C22]). 
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1 - grid; 2 - laminar boundary layer 

Illustration 5.44; Pressure Gradient Influences with Identical 

Plate Flow Reynolds Number (Re^ = 6.4 x 
10 J ) and Varying Free Stream Turbulence 
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In addition to the changed reversal behavior, there are - in the 
•forward portion of the plate - distinct increases of the laminar 
heat transfer coefficients (Illustration 5.44a and b), which can 
also be attributed to the influence of free stream acceleration. 
Its cause can be explained by steeper profile gradients at the wall 
which occur in accelerated flows. In comparison to test values 
without pressure gradients, these increases are approximately 20-30% 
near the leading edge of the plate in the laminar heat transfer. 
However, they decrease as the running length increases. The decay 
behavior of the Stanton number with pressure gradients is slower and 
more viscous than in unaccelerated situations. The reason for this 
effect is also due to increased stability of accelerated laminar 
boundary layers. 

Illustration 5.44c rounds out the basic study on the influence of 
negative pressure gradients. In contrast to the above examples, 
flows with very high external turbulence are analyzed. The high 
turbulence intensities generated with grid 4 show no stabilizing 
effects on the reversal caused by the pressure gradients. In all 
cases, the reversal occurred immediately behind the leading edge of 
the plate. Due to intense free stream turbulence, there is no re- 
laminization tendency, even with high acceleration of k = 3.2 x 

^ IDdX 

10 . Therefore, Illustration 5.44c describes only the turbulent 
heat transfer. These tests show that, with increasing free stream 
acceleration, there is a decrease of the local turbulent Stanton 
number. Similar pressure gradient effects were identified by Kays 
and Moffat C863 in turbulent boundary layers without external 
turbulence. Therefore, the results confirm that, even with maximum 
external turbulence, acceleration causes an attenuation of the 
turbulent exchange events in the boundary layer (especially near the 
wall), which results in a decrease of local heat transfer 
coefficents. If the description of turbulent exchange events in the 
boundary layer is made with mixing method equations (see intengral 
procedures, Chapter 6), these changes are theoretically recorded in 
the boundary layer turbulence with increased van Driest wall 
attenuation constants. 
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With a detailed comparison of measured heat transfer distributions 
both with and without pressure gradients it can be demonstrated that 
pressure gradient influences generate significant changes in the 
reversal behavior and in laminar and turbulent heat transfer. They 
are therefore very important for the computation of accelerated 
blade boundary layers. The causes can be attributed to dampening 
effects of negative pressure gradients on boundary stability and 
boundary layer turbulence, as they are described in classic boundary 
layer theory and as demonstrated in previous experimental studies 
with low free stream acceleration. Based on this, a study is made 
in the following chapters on how these effects are dependent on the 
respective flow conditions (inlet Reynolds number, free stream 
turbulence). I will also show additional effects if very high 
acceleration intensities, as they occur in gas turbines, are applied 
to a boundary layer with high free stream turbulence. 
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5. 5. 3. 2 Reynolds Number and Turbulence Intensity Influences 
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In gas turbines, Reynolds number effects always occur when operating 
conditions are changed. They are therefore of great significance to 
the turbine’s operation. 

The acceleration parameter k is inversely proportional to the 
Reynolds number 

k s 1/Re 

An increased Reynolds number will, therefore, result in increased 
acceleration values. This is shown by Illustration 5.45a for the 
tests with contour 1, and by Illustration 5.45b for the tests with 

contour 2. When the flow Reynolds number decrease from Re. = 6.4 x 

5 5 Lo 

10 to Re. = 3.3 or 3.5 x 10 , the k-values of contour 1 increase 

L ° -6 -b 

from 1.2 x 10 to 2.3 x 10 , and the k -values for contour 2 

, • max , 

from 3.2 x 10 to 5.7 x 10 . The flow with k =5.7 shows 

max 

those acceleration intensities, which are representative for 
realistic blade flows. Events in the heat transfer and boundary 

layer behavior caused by Reynolds number changes are shown for the 
flow with contour 1 in Illustration 5.46 and for the flow with 
contour 2 in Illustration 5.47 . 

Tests conducted with contour 1 with varying plate flow Reynolds 

numbers and varying free stream turbulence are based on low 

(Illustration 5.46a) and medium (Illustration 5.46b) acceleration 

-6 -6 

intensities with k = 1.2 x 10 and k = 2.3 x 10 . A special 

characteristic of these flows is that there are acceleration 

parameters (see Illustration 5.45, sink-flow) which are almost 

constant in the test section. As expected, the increase of the 
pressure gradient parameter k, which is caused by the decreased 
Reynolds number, results in an increasing boundary layer stability. 
With unchanged flow turbulence, this boundary layer stability 
results in an additional down-stream re-location of the reversal. 
It may, as shown in situations with low free stream velocity (Tu q = 
1.6% and 2.3%), even supress the reversal completely. The absolute 


level of laminar and turbulent Stanton numbers will also be 

changed. With lower Reynolds numbers (Illustration 5.46b), the heat 

transfer coefficients - in the laminar case (see Tu q = 1.6 - /.) - may 

be up to 30% higher. In turbulent boundary layer sections, however, 

there is a comparatively low increase of approximately 5%. We 

already know about lower Reynolds number dependency of turbulent 

heat transfer in unaccelerated flow situations: In unaccelerated 

laminar boundary layers the local Stantonnumber is proportional to 
-0 5 

Re ' . In turbulent boundary layers without pressure gradient, 

X - 0.2 

there is only a dependency on Re (see my test results in 

M 

Illustrations 5.15 - 5.17, Chapter 5.4). 




Illustration 5.46: Reynolds Number Influences with Low 

Acceleration and Varying Free Stream 
Turbulence (Contour 1) 
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In flows with pressure gradients (Illustration 5.46), this 
difference is emphasized by the following: In accordance with 
Illustration 5.44, the increase of the pressure gradient parameter 
k(x), which is a function of the decreasing Reynolds number, results 
- in laminar boundary layers - in an additional increase and - in 
turbulent boundary layer sections - in an additional decrease of the 
Stanton numbers. 

For both Reynolds numbers, the heat transfer distributions shown in 
Illustration 5.46 give no indication of noticeable dependencies of 
the turbulent heat transfer on free stream turbulence. 
Additionally, it was determined that, due to the thin leading edge 
of the plate and despite acceleration with high flow turbulence, the 
reversal occurs in the area of the uncooled leading edge (see 
Chapter 5.5). These observation are in complete agreement with the 
results obtained in unaccelerated boundary layers with a local 
application of the Stanton number (see Illustrations 5.15 - 5.17). 
However, events in the laminar boundary layer sections during 
acceleration in contrast to the tests without pressure gradients. 
Illustration 5.46 (x < 0.12 m) shows increases of the laminar heat 
transfer as a factor of turbulence intensity. This was observed in 
previous studies with lower pressure gradient intensities (see C35, 
383 and Chapter 2.2.1). Increase rates are up to 157. if, as shown 
in Illustration 5.46b, the flow turbulence is increased from 1.67. to 
4%. These turbulence intensity influences are greatly dependent on 
the relative strength of free flow acceleration and, therefore, on 
the respective flow Reynolds number. As shown by the measurements 
with Tu q = 2.37. in the forward section of the plate, turbulence 
intensity dependent increases of the heat transfer (compared to Tu q 
= 1.67.) are almost completely suppressed by higher acceleration 
(Illustration 5.46b). 
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The increases of the laminar heat transfer, which are induced by 
external turbulence, exist only when negative pressure gradients are 
present. This effect was not observed in tests with constant 
external velocities. Without free stream acceleration, the required 
stabilization of the laminar boundary layer does not exist . 

Therefore, the disturbance caused by external turbulence results 
directly in a laminar-turbulent boundary layer reversal. 

In order to provide a better classification of this new flow 
situation with free stream acceleration, I would like to refer to 
the studies by Junkhan and Serovy C353. Junkhan and Serovy C351 
measured similar pressure gradient effects with comparable increases 
of the laminar heat transfer in less accelerated boundary layers. 
The boundary layer tests conducted under these conditions show that 
fuller velocity profiles occur which, despite the influence of outer 
turbulence, maintain the typical boundary profile form described by 
Hartree 11081. Due to the comparable heat transfer increase << 
1551) , this result can be applied to our tests. This means that, 
under the conditions described here, the laminar properties of the 
boundary layer continue to prevail despite the influence of free 
stream turbulence. Therefore, the flow can be assigned to laminar 
boundary layers. 

I have presented a detailed discussion of Reynolds number and 

turbulence intensity influences for simple flow forms with constant 

acceleration parameters and moderate pressure gradient intensities. 

The question now arises how the measured effects can be applied to 

the more complex flow conditions on the vacuum side of gas turbine 

blades. In contrast to the cases shown in Illustration 5.26, this 

flow is characterized by pressure gradient matrices with down-stream 

variations and by maximum acceleration intensities. In order to 

demonstrate the resulting additional effects, free stream velocity 

distributions were generated with contour 2 on the vacuum side 

(Illustration 5.45). In addition, an analysis was made of heat 

transfer distributions measured in a comparison with varying plate 

5 5 

flow Reynolds numbers (Re. = 6.5 x 10 and Re. = 3.5 x 10 ) . The 
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results are summarized and discussed in Illustration 5.47. 


Illustration 5.47b describes boundary layer flows which, due to 


their maximum acceleration intensities 
particularly relevant to turbines. 


<k 


max 


= 5.7 x 10 i , are 


The comparison in Illustration 5.47 shows that, with a high Reynolds 

number (Illustration 5.47a), an increase of the flow turbulence from 

Tu = 2.3'/. to Tu = 47. will result in a noticeable relocation of 
o o 

reversal initiation toward the leading edge of the plate. With a 

lower Reynolds number (Illustration 5.47b), this is completely 

suppressed by the higher pressure gradient. In the entire 

acceleration distance (x < 0.25 m) , all measurement points with less 

than 4% flow turbulence are in accordance with the matrix described 

in Illustration 5.46 with low turbulence intensity dependent 

increases and a laminar heat transfer. Only the high turbulence 

intensities of grids 3 and 4 (Tu = 77. and 107.) - with the turbine 

o 

pressure gradient conditions in Illustration 5.47b - result in a 

higher influence of free stream turbulence. Starting with the 

beginning of the plate, they effect a drastic increase of the local 

Stanton numbers (up to 1007.). However, in their path along the 

plate, they continue to be similar to laminar heat transfer. The 

high turbulence intensities induce distinct disturbances in the 

boundary layer which cannot be attenuated despite maximum 

acceleration (k = 5.7 x 10 . ) . In our tests, a balance was 

observed as result of the destabilizing effect of external 

turbulence and the stabilizing influence of accceleration. This 

balance causes large portions of the boundary layers to be in a 

stable, laminar-turbulent intermediate condition; it therefore has 

the typical properties for boundary layer reversal. The 

distribution of the Stanton numbers measured in Illustration 5.46a 

for 4% flow turbulence can also be regarded as confirmation of the 

potential occurrence of stable, laminar-turbulent intermediate 

conditions. In this flow situation, the measurement points are 

evenly distributed throughout the entire forward section of the 

plate (x < 0.12 m) in the center between laminar (Tu = 1.67.) and 

o 

turbulent (Tu = 107.) heat transfer. 
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Illustration 5.47: Reynolds Number Influences with High 

Acceleration and Varying Free Stream 
Turbulence (Contour 2) 


According to our experiences, this special behavior of the boundary 
layer, which was not observed in tests with constant acceleration 
(see Illustration 5.45), occurs when high turbulence intensities 
exist simultaneously with intensive acceleration. However, another 





■factor for these stable, reversal-1 ike boundary layer conditions is 
that, as with flows modelled on the vacuum side, the acceleration 
intensity increases with increasing running length (see Illustration 
4.45b). This is the main reason why further excitation of the 
boundary layer turbulence is suppressed. This excitation, in an 
unaccelerated or constantly accelerated situation (Illustration 
5.46), will result in a continuation of the reversal process. 

A systematic comparison of tests with varying free stream 
temperature conditions has thus shown that, for the initiation of 
heat transfer on gas turbine blades, absolute turbulence and 
pressure gradient intensities are not the only factors to be 
considered. A significant role is also played by the acceleration 
matrix along the blade surface. Certain combinations of these 
influencing parameters may result in flow situations which will have 
to be re-classified. In these situations, the heat transfer can 
neither be identified as laminar nor as turbulent. In large 
sections, it will obtain a laminar-turbulent intermediate 
condition. This knowledge can be used to explain the 
extraordinarily high increases of local heat transfer, which were 
observed in blade tests with high free stream turbulence. Due to 
their similar matrix they have, so far, been interpreted as 
"increases of the laminar heat transfer" (see E993). They are long, 
highly turbulent boundary layer events with a reversal character. 
They are a result of an interaction between free stream turbulence, 
highly negative pressure gradients and, especially, an acceleration 
matrix with a down-stream variation. Their reliable prediction is 
presently the main problem in the design of cooled gas turbine 
blades. The complexity of this flow event is not recorded with 
sufficient accuracy in the numeric computation procedures used for 
turbulence models. 
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5. 5. 5. 3 Correlation Equations -for Calculating Turbulent Heat 
Transfer in Accelerated Flow 

Because correlations are the simplest and quickest method to 

calculate heat transfer, there is a desire to use them in as many 

flow situations as possible. The numerous tests in this study 

provide the opportunity to examine and further develop correlation 

equations for turbulent heat transfer with respect to their use in 

highly accelerated flows with high external turbulence. My studies 

are based on the previously described heat transfer distributions 

with varying acceleration and high free stream turbulence, which 

were measured with a temperature ratio of T /T = 0.64. As stated 

w « 

in my analysis of reversal events, intensive wall cooling has no 
additional influence on boundary layer behavior (see Chapter 5.5.2, 
Illustration 5.36). Its presence only increases the absolute heat 
transfer, so that this additional effect can be considered 
separately. 

In Illustration 5.48 . the measured heat transfer coefficients are 

again collated as a function of the local running length Reynolds 

number. Illustrations 5.4Ba and 5.48b provide comparison data for 

the tests conducted with contour 1 and acceleration values of up to 

k = 2.3 x 10 . In the forward Reynolds number range, the Stanton 

numbers for Tu = 47. in Illustration 5.48a and for Tu = 77. in 

o o 

Illustration 5.48b represent portions of the reversal process. 

Illustrations 5.48c and 5.48d contain available comparison 

situations for medium and maximum acceleration intensities. They 

can be assigned to the tests with contour 2 (k =3.2 x 10 6 and 

max 

5.7 x 10 ). With these flows, acceleration of the boundary layer 

5 

is replaced by a subsequent delay, starting with Re « 3-5 x 10 . 

To correlate the turbulent heat transfer coefficients, I used a 
relationship by Kays (see C663, p. 248). It was originally deduced 
from Moretti’s Cl 123 data for almost isothermal boundary layers 
without external turbulence. It describes the local Stanton number 


St = 


as a -function o-f the Prandtl number Pr, the running length Reynolds 
number Re, and the local pressure gradient parameter k(x). This 
correlation, expanded by the additional influence of wall cooling, 
is as follows: 



1 - correlation equations in acc. with eq. 5.12 

Illustration 5.4B: The Use of Correlations for the 

Calculation of Heat Transfer in Accelerated 



The value F., which occurs in equation (5.12), represents the effect 

of the pressure gradients. Fj represents the temperature decrease 

at the wall. is assigned the mean value 1.08 for test 

comparisons because, according to the heat transfer distributions 

without acceleration (Illustrations 5.15 - 5.17), the increase of 

local Stanton numbers with T /T = 0.64 must be approximately 6-107.. 

w ™ 

Kays C663 specifies the Reynolds number in equation 5.12 by an 
integral reference value 
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and assigns the value 165 to the pressure gradient factor F^. 
Illustration 5.48 shows the resulting matrix of the Stanton numbers 
for all test situations (hatched line). As demonstrated by the 

comparisons with tests, the heat transfer is calculated in a 
satisfactory manner for situations with constant acceleration 
(Illustrations 5.48a and b). In the higher accelerations with 
contour 2 of Illustration 5.48c and, especially in Illustration 

5.48d, greater deviations from the test values can be observed. In 
accordance with the statements made in the previous chapter (see 
Illustration 5.47b), measurements in Illustration 5.47d do not 
represent flow events in the forward Reynolds number range, but flow 
events which are similar to reversal in the boundary layer. With 
this comparison, Kays’ relationship for turbulent boundary 
layers predicts heat transfer coefficients which are too low. It 
therefore becomes obvious that the correlation considers the 

atten uating effects on boundary layer turbulence with very high 
acceleration as too large. The reason for this is the influence of 
free stream turbulence, which was not considered. External 

turbulence counteracts this atten uation process which results in a 
delayed re-laminization initiation. 


The observed inaccuracies and the fact that an integration to 
determine the Reynolds number according to Kays [663 requires great 
effort were the reasons I sought a development which eliminates the 
calculation of the reference Reynolds number, and which accurately 
describes the turbulent heat transfer as a direct function of the 
local running length Reynolds number. This was achieved by an 
adjustment of the pressure gradient factor F^, with unchanged use of 
the formula context of equation (5.12). As shown by Illustrations 
5. 48a-c, the results obtained are in very good agreement with the 
test data by introducing the constant value = 90. With this (127 
procedure and the use of the local running length Reynolds numbers 



equation (5.12), in Illustrations 5.48a-c, provides an accurate 
matrix and a very realistic absolute value for the turbulent heat 
transfer (max. deviation 6-10'/.). The remaining increases in the 

forward Reynolds number range can be explained with post- 

transitional effects, which are a result of the high up-stream 
increases of the laminar heat transfer caused by acceleration (see 
Illustration 5.44a and b). 

In contrast to the procedure used by Kays [663, which requires a 
greater effort, this expansion of the correlation also improved 
calculation results. Illustration 5.48c clearly shows that, in high 
acceleration, the influence of free stream turbulence considered by 
adjusting F. and the use of local Reynolds numbers results in more 

K 

accurate reproductions of the measured heat transfer coefficients. 
It is further pointed out that these measures resulted in a much 
better agreement with the test data for the transitional flow 
situation with highest (turbine-relevant) pressure gradients 

(Illustration 5.48d) in the forward acceleration range (Re < 3 x 

5 x 

10 ). 
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In order to re-confirm the practicability of explicitlyconsidering 
pressure gradient influence in correlations, Illustrations 5.48a-d 
shows the heat transfer coefficients (hatched lines), which were 
obtained with equation (5.12) by neglecting the pressure gradient 
terms (F^ = 0). All Stanton numbers calculated with this 

simplification are too high and their matrix is not accurately 
described. A satisfactory agreement with test data exists only with 
the very weakly accelerated boundary layer with k = 1.2 x 10 ^ 
(Illustration 5.48a). This no longer applies for flows with higher 
velocity gradients. In this case, the correlation which considers 
the pressure gradient as an additional factor is clearly more 
favorable (Illustrations 5.48b-d). 

With this systematic study of pressure gradient influences plus the (128 
examination and further development of correlations for computing 
turbulent heat transfer in accelerated flows with high external 
turbulence, my analysis of experimental results (see C1233, C1243) 
is concluded. Essential new information achieved on the influence 
of free stream turbulence, wall cooling and pressure gradients was 
based on thorough tests. In order to provide a good overview of 
this subject, all test data and the results obtained were thoroughly 
documented. This thoroughness fulfilled an important prerequisite 
for the examination and further development of complicated boundary 
layer calculation methods (integral and differential procedures). 

My goal in the following chapters of this dissertation was to 
develop an integral procedure based on these experimental results, 
and to use it as an alternative to differential methods for an 
accurate, flexible computation of heat transfer in cooled gas 
turbine blades. 
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6*. DEVELOPMENT OF AN INTEGRAL PROCEDURE 10 COMPUTE TWO- 

DIMENSIONAL BOUNDARY LAYERS WITH REVERSAL. HEAT TRANSFER. 

PRESSURE GRADIENTS AND FREE STREAM TURBULENCE 

6. 1 Procedure Selection 

The complexity of experimental -flow events under the influence of 
free stream turbulence, wall cooling and pressure gradients requires 
highly flexible calculation procedures. These properties are 
contained in field methods, which provide a numerical solution for 
differential conservation equations in connection with a turbulence 
model (see C33). This evidence has not yet been established for 

integral procedures. So far, these procedures failed due to the 
description of the long reversal events and due to the separate 
calculation of velocity and temperature boundary layers. According 
to previous studies (see Chapter 2.2) and my own determinations, 
these two properties are an important factor in a calculation 
program, which is to predict the blade heat transfer in extremely 
accelerated boundary layers. 

In this dissertation I will show for the first time that the basic 
prerequisites for calculating transitional boundary layers in gas 
turbine blades can also be applied to integral procedures. I will 
also demonstrate that these integral procedures, just like 

differential procedures, are suitable for an accurate description of 
flows in blade boundary layers with high free stream turbulence, 

intense wall cooling and maximum pressure gradients. 

The starting point for this calculation model areknown conservation 
equations for impulse and kinetic energy used to determine the flow 
boundary layer. Two additional integral conservation equations were 
developed for the separate calculation of the temperature boundary 
layer and heat transfer. In analogy to the impulse and kinetic 

energy equations, they describe the enthalpy exchange in the 
boundary layer. In order to complete the equation system, I used a 
deary arranged 1-equation turbulence model, which requires little 
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effort and which was tested on numerous occasions in connection with 
differential methods (see C4, 35, 473). The model can be 

incorporated in the calculation algorithm, because it performs the 
integration of the balance of turbulence exchange events in the 
boundary layer. Turbulent shear stress distributions were computed 
with the aid of a mixing method hypothesis. 

With integral procedures, formulations for the velocity and 

temperature profiles in the boundary layer must be specified for a (130 
solution of the conservation equations. They form implicit 
definitive equations for wall tangential stress and the heat 
transfer at the wall. For this calculation, universal profile 
formulations were made for laminar, reversing and turbulent boundary 
layers. They can be used even in the most extreme acceleration 
situations found in turbines. When compared to field methods, the 
availability of these boundary layer matrices is the only 

flexibility restriction with these newly developed integral 
procedures. 

The boundary layer in the reversal is modelled as a function of the 
boundary layer intermittance matrix. Profile formulations as well 
as the intensity of the boundary layer turbulence are determined by 
the respective value of the intermittence factor. 

The description of the calculation procedure begins with the 
derivation of integral conservation equations for velocity and 
temperature. I will then introduce the profile formulations and the 
turbulence model. On this basis, I will model the boundary layer 
reversal . 

6.2 Basic Equations 

The theoretical description of boundary layers is based on the 
differential conservation equations for mass, impulse and energy. 

Using the prerequisite of a two-dimensional, stationary flow - which 
applies to blade circulation - we obtained the following parabolic 
differential equation system: 
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continuity equation: 


3(P'U) + 3(p‘V) c 
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impulse equation: 
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energy equation: 
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( 6 . 1 ) 


( 6 . 2 ) 


(6.3) 


The variables in the conservation equations, like the velocities u 
and v, the fluid density p and the total enthalpy h Q , represent 
chronological mean values. The diffusion terms on the right side of 
the balance equation are mainly determined by the tangential stress 
r, the pressure gradient dp/dx (which is constant throughout the 
boundary layer) and the local heat flow distribution q. 

The required integral conditions for flow and temperature boundary 
layers and the prerequisite parameters can be totally derived from 
these basic equations. The impulse and energy equations (equations 
6.2 and 6.3) are multiplied by "weighted functions". They are then 
integrated through the wall distance and replaced by two common 
differential equations of the first degree. A large amount of 
literature (e.g. C20, 673) is available for obtaining integral 
equations in this manner for the computation of the flow boundary 
layer. Referred to as impulse and kinetic energy equations, they 
are assigned to the "dissipation" integral procedure and are 
expressed as follows: 


impulse equation: 
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kinetic energy equations 


d63 

dx 
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The respective definitions for integral boundary layer values and 
form parameters are as follows: 
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(6.6) 


The integral relationships required for describing the temperature 
boundary layer can be obtained with the basic equations in a similar 
manner. In accordance with a proposal made by Patankar and Spalding 
[1231, the energy equation (6.3) is expanded with the weight 
functions (h Q -h oM ) n (n = 0; 1). After eliminating the diagonal 

velocity v, it is partially integrated with equation (6.1) over the 
wall distance. This results in two enthalpy equations, which are in (132 
accordance with the flow equations (6.4) and (6.5): 


enthalpy equation 1: 
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enthalpy equation 2: 
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with 



The turbulent exchange events are only considered in the dissipation 
terms c D and c^ on the right side of the balance equations (6.5 and 
6.8). In addition to the molecular transport, the tangential stress 
r and the heat -flow q contain the transport components caused by the 
turbulent -fluctuation movements. They are as -follows: 


( 6 . 10 ) 


du , i 

• (3y ' ’ U V 

/ \ du 

= + u t> ' ay 
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ay P 
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( 6 . 11 ) 


The turbulent components are expressed by the cross-correlation 
terms -pu’v' and -pc^’v’. They occur as an additional unknown 
quantity of the equation system. In analogy to laminar transport, 
they are replaced by the introduction of a turbulent viscosity 
and a turbulent heat conductivity . A connection between the two 
values is made by determining the turbulent Prandt number 


Pr 


t 



( 6 . 12 ) 


which describes the relationship between the turbulent impulse and 
heat exchanges. When the turbulent Prandtlnumber is known - which 
can be replaced with Pr^ s 0.9 = const, for all boundary layers 
flows - the determination of turbulent exchange events is reduced to 
the description of u’v 5 or the apparent viscosity In our 
procedure, this is accomplished with the aid of a 1-equation 
turbulence model (see Chapter 6.4). 
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In addition, -formulations for the velocity and temperature profiles 
are required for a solution of the differential equation system 
based on integral calculation. Velocity and temperature matrices 
are required for the calculation of the integral boundary layer 
values in equations (6.6) and (6.9) and for the determination of the 
dissipation integral and c^. At the same time, the profile 
formulations provide the definitive equations for the friction 
coefficient c^. , which is required in equation (6.4), and the Stanton 
number St, which enters the integration equation (6.7). 

6.3 Profile Formulations 

6.3. 1 Boundary Laver Profiles ill Flows with Constant Material 
Values 

Integral procedures are characterized by the fact that suitable 
profile formulations are required for the determination of velocity 
and temperature matrices. Similar to differential procedures with 
’wall functions’ (see C783), known parameters enter the calculations 
with the profile formulations. In accordance with conservation 
equations (6.2) and (6.3), they must be completed in the following 
manner for the flow near the wall (y = 0, u = v = 0): 
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, 3’U, dU 6 
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(6. 14) 


The observance of 

these ’wall conditions’ 

is achieved by a square 

formulation 

for 

the velocity profile 

at the 

wall; for the 

temperature 

profile 

it is obtained by a linear matrix 

as fallows: 


u(y) = k x 

y + k 2 x y 2 

(6.15) 



T - 

T = k 
w 3 

; x y 

(6.16) 
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Regardless whether the boundary layer i laminar, transitional or 
turbulent in character, this determines the viscous ■flow near the 

4 * ^ 4 * 

wall. In the similarity coordinates u , T and y , which are use-ful 
for a general illustration of boundary layer profiles, the velocity 
and temperature matrices, which are described by equations (6.13) 
through (6. 16) and the paramaters 


. - t 3U> 

T w * (■* ' 


ay' 

J y=o 


% - • w> 


y=o 


are in accordance with the followings 


(6.17) 

(6.18) 


(6.19) 

( 6 . 20 ) 

With the applicable definitions 


u + = y + + k 2 *-y +2 
T + = Pr • y + 



(T-T.j t )-p-U T .c 

( 6 . 21 ) 

q w 

( 6 . 22 ) 

(6.23) 


the analysis of the boundary layer profiles is initially 
concentrated on the critical case of flows with constant material 
values. 


Based on equations (6.19) and (6.20) for flows near the wall, 
analytical profile matrices were determined separately for laminar 
and turbulent flows which describe the total flow up to the edge of 
the boundary layer. By superimposing the laminar and turbulent 
profile formulations we are now able to model the boundary layer 
profiles for the laminar-turbulent reversal range. 
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6.3. 1. 1 Laminar Profiles 


A modelling of the reversal and an inclusion of the acceleration 
intensities, which are representative for blade circulation, as 
proposed by Polhausen and Hartree (see C1093), were not possible 
with known profile formulations for laminar boundary layers. In the 
early development stages of the calculation program I determined 
that these matrices, which are based on simple polynomial 
formulations, are too inflexible. Therefore, they result in 
inappropriate overshoots of the boundary layer profiles, even in 
flows with relatively low acceleration. This overshoot behavior 
also occured when these formulations were used for reversal 
modelling (see following chapters). 

In my newly developed version, the laminar boundary layer profiles, 
which can also be used as starting profiles, are in accordance with 
the schematic shown in II lustration 6.1. Temperature and velocity 
profiles were determined with identical formulas. 



Illustration 6. 1; The Laminar Profile Formulation 
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In the section near the wall, the laminar velocity pro-files are in 
agreement with the asymptotic matrix of equation <6.19). This was 

confirmed by my boundary layer tests (see Illustration 5.7). With a 

given wall distance y*, this matrix is replaced by a polynomial of 
the 3rd degree with logarithmic variables. The coefficients of the 
polynomial are determined by the locations and gradients near the 

validity limit for y*" and at the edge of the boundary layer. The 

conditions for equation (6.19) are achieved at the y = y+ 
location. At the edge of the boundary layer, the condition 


U = U. 



must be considered. 


(6.24) 
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After solving the polynomial, the formulation for the entire laminar 
velocity profile is as follows: 


+ + . , * 
u = y + k 2 


(0 ' y + < y+) 
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+ 


+ 

u 6 


K-n + (3-2-K).„ 2 + ( K- 2 ) - n 3 


(y T > y^) 


(6.25) 


with 


(1+2- kj-y+) . ln(6 + /y!) 


K - <) 


and the dimensionless wall distance 


In (y + /y*) 
In (« + /y+) 


(6.26) 


(6.27) 


With minor changes to the velocity profiles, which are mainly due to 
the differences in equations (6.19) and (6.20), the laminar 
temperature profile can be illustrated in a similar manner: 


r = Pr • y 


(0 < y < y! 


(6.28) 
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with 


K t = 


Pr-y* • In («W ) 
s t t 

a: -t!) 


(6.29) 


and 


In (y + />s ) 
^ = 1" 


(6.30) 


The range limit y+ ( y +^) of the dual layer boundary layer pro-file 
determined in the above manner was, in comprehensive test 
calculations, determined for both unaccelerated and accelerated 
flows. To avoid an overlap of the profile in the outer range with 
the asymptote given in equation (6.19) or (6.20), the following 
distinction is made: 

y* = max { 11 ; y + (u + = 0,25-u*)} 

(6.31) 

yt = max { 11 ; y + (T + = 0,25 T* )} 

s t 6 t 


6.3. 1,2 Turbulent Boundary Laver Profiles 

The formulations for turbulent boundary layer profiles are based on 
the 3-layer model proposed by Walz and Neubert C1261. This model 
determines all three boundary layer area - laminar bottom layer, the 
area of the logarithmic wall theory, and the wake portion - 
cohesively and by including the wall relationship given by equation 
(6. 19) : 

(6.32) 

with C = 5.1; K= 0.41; a = 0.3. 


u + - {(1- I - c-a).y + + K* *y +2 - C}-e 
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C and form the constants of the logarithmic wall theory. In this 
profile formulation, the buffer layer constant "a" results in an 
even transition from the laminar bottom layer to the logarithmic 
range. Its value a =0.3 was determined by Walz and Neubert with 
the aid of a test comparison. The follow-up function w(y/6) is 
modelled according to Hinze 

w(y/«) = j (l - cos (n |)) (6.33) 

In equation (6.32) it is weighted with the applicable wake parameter 
P/t\ . P/ v* represents the increase of the velocity profile to the 
logarithmic wall theory at the edge of the boundary layer: 

p - = u + 5 - (i • ln6 + + C) (6.34) 

With the pressure gradient term K^xy"**, the profile formulation 
(equation 6.32) meets the complete wall connection theory (equation 
6.19) and can therefore be used for accelerated flow. Based on a 
great variety of experimental studies (see Chapter 2.2), no change 
to the logarithmic wall theory is introduced if the calculation is 
made for flows with pressure gradients. The same applies to the 
temperature profiles, which is discussed in following chapters. 

Because turbulent temperature profiles - like the velocity profile - 
are characterized by laminar bottom layers, the logarithmic wall 
theory and the wake area, equation 6.32 can be used for their 
description. By considering the flow conditions near the wall <T + = 
Pr x y + , equation 6.20), which are different from those in equation 
6.19, the following formulation is established for the turbulent 
temperature profile: 

T+ * { ( pr - ^ - c t -a)-y + - C t }.e 3 y + i_ . in (y + + 1) + 

+ C T + (1 - cos(n f-))-[T+ - j|- • lns + - c t ] 

t t t 

with Cj = 3.B; = /Pr^j Pr^ = 0.9. 


(6.35) 
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The constants of the logarithmic wall theory C T and v\j., which are 
contained in equation 6.35, are determined independent of the 
respective cooling condition in flows with almost constant material 
values (cooling is introduced to the profile formation by 
transformation formulations). Experimental profile studies (see 
Chapter 5.5. 1.2) have shown that, in this case, C^ has the value 
3.8, as proposed by Sill C783 and Simonich and Bradshaw C543. 

6. 3. 2 Boundary Laver Profiles in Flows with Changing Density (139 


The formulations for temperature and velocity profiles were 
initially determined only for flows with almost constant material 
values. However, profile studies conducted by me and other authors 
showed that, with the aid of suitable profile transformations, these 
formulations can, in an unchanged form, also be applied to highly 
cooled boundary layer flows. This procedure will be described in 
the following paragraphs. 


In order to perform profile transformations, the profile 
formulations introduced with the similarity coordinates u + = f(y + ) 

4 * + 

and T = f(y ) must be expressed in a more practical manner: 


u 





(6.36) 


This can be achieved, as proposed by Walz [126], by forming the 
easily reconstructable interrelations for u r , y + and u + 


u , = v /c f /? 



JJ 

v /c 7 2 
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Re 2 -/y2 


(6.37a) 


(6.37b) 


They are then inserted in the original profile formulations u + = 
f(y + ) (equations 6.25 and 6.32). The laminar and turbulent 
formulations for the velocity profiles thus exhibit a structure 
which is implicitly dependent on and Rej 

U7 * * Re z • c f) (6.38) 
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and -for which available transformation procedures can be directly 
used. 

The transf ormation for the calculation of boundary layers in flows 
with changing density relates to the method described by Walz 
[126]. Within the profile context given by equation (6.38), c^ and 
Re ? are replaced by material value correcting expressions. The 
Reynolds number Re 2 is corrected (with respect to temperature) with 
the free stream density and the wall viscosity in the following 
manner s 


Re 2 


'W 4 ? 


(6.39) 


The actual tangential stress coefficient c^. is replaced by an 
incompressible reference value c^ . 


In accordance with van Driest (see Kiister C1283), the c^ 
transformation for high free stream velocities is performed as 
follows: 


r - 

c f 



(6.40) 


with 


A = 






= 1.4 


2A 2 - B 


/4A 2 + B 2 




B 


/4A 2 + B 2 


i) ; 


r = Recovery-Faktor = 0.88 


In flows with non-constant material values and Ma ■* 0 (in our case: 
Ma < 0.3) Kiister’s C1283 simplified transf ormation is used for c^.: 



(6.41) 
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The transformation equations for c^ can also be used for laminar 
boundary layer flows, although these interrelations were only 
developed for turbulent boundary layers. The simplification is 
based on available literature (see Kays C643), which states that 
cooling influences in laminar and turbulent boundary layers are 
generated in an approximately identical manner. 


Similar to procedures used for the velocity profile, the inclusion 

of variable material values and the transformation of the 

temperature profile also occurs within the context described by 

equation (6.36). The relationship provided by equation (6.37a) for 

u p in included in the development of this general temperature 

profile formulation in an unchanged manner. The wall distance 
+ 

coordinate y is adjusted to the reference values of the enthalpy 
equation (equation 6. 7). The impulse loss thickness 6 0 , which 

characterizes the flow boundary layer, is replaced by the enthalpy 
loss thickness S,_ (equation 6.9), The boundary layer thickness 6 is 

J 

replaced with the thickness of the temperature boundary layer 6. . 
✓ t 

This results in the -following: 


y 


+ 




Re 5 • 


(6.42) 


The transf ormation of T + , which is required for the determination of 
the temperature profile in accordance with equation (6.36), will - 
by including the Stanton number St defined by equation (6.7), result 
in: 


* 



(6.43) 
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This transforms the profile formulations T* = f (y + ) (equations 6.28 
and 6.35) into 


T - T w 

n -t 
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■ F (f 


Re* 


st x ) 


(6.44) 
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When considering variable material values, the transformation 
principle used for the velocity profile can be used. In the new 
profile illustration according to equation (6.44), Re,, is formed 

J 

with an average viscosity 



The Stanton numbers and c^ values in reference to free stream 
density are replaced by the values St^ and c^., which were 
transformed with equations 6.40 or 6.41. 

6, 4 The Turbulence Model 

A great number of turbulence models have been introduced for the 
determination of turbulent exchange events in boundary layers. 
Depending on the number of transport equations used for 
characteristic turbulence values, they are categorized into 0- 
equation, 1-equation, 2-equation and Reynolds stress models. The 
complexity and the general validity increase with the number of 
model equations. To find a good compromise between general 
validity, complexity, availability and user friendliness, I would 
like to use a 1-equation model developed by McDonald and Fish C51 
and McDonald and Kreskovsky C493. This model presents numerous 
expansion possibilities and its numeric solution is very simple. 
This type of model uses F'randtl’s mixing method for the 
determination of the turbulent tangential stress distribution. The 
mixing path at the outer edge of the boundary layer is deduced from 
the integrated balance equation for turbulent, kinetic energy. This 
method has already been successfully used in the calculation of heat 
transfer in gas turbine blades (see C53). In the following 
calculation procedure I will use the unchanged turbulence model for 
turbulent boundary layers. I will, however, introduce an expansion 
to the model for the calculation of the boundary layer reversal . 
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The starting point for the turbulence model is the known (142 
differential transport equation for turbulent kinetic energy q. 

After a partial integration over the wall distance y, it is 
transformed into an integral condition for turbulence energy. 
According to McDonald et al C5, 49], it is as follows: 


1 d 

2 <Ix 


p*u*q dy 
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Konv. 
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1 -pu ' v ' dy - | p-f -dy 
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( p( 72 - 72 ) . |H d y + E ( 6 . 46 ) 

o 

Norm. -Sp.- Prod. Entr- 



The individual terms of equation (6.46) are, from left to right: 
convection, production, dissipation, normal stress production and 
’entrainment’ of turbulent kinetic energy. The entrainment portion 
E consideres the turbulence energy (q. = 3 x Tu 2 x u^ - ) , which 

w w 

penetrates the boundary layer from the outside. It explicitly 

registered influences on free stream turbulence. 

The connection between the kinetic turbulence energy q, which is 

contained in equation (6.46), and the additional turbulence 

components u’, v’ , u’v’, and € is defined with additional model 

equations: 

(d - f (^) • q { ) • «i = - FF 1 

F 2 * a 2 • q ; v T2 = a 3 • q ; f(£) = (1 - cos (n • |))/2 j (6.48) 

e = (-^V) 3/2 / L 
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With these -formulations, the integral energy equation (equation 
6.46) is transformed into a correlation which is only dependent on 
the shear stress u’vMy/6). In turbulent boundary layers with high 
Reynolds numbers, the structural parameters a^, and a^ have the 
values 0.15, 0.5 and 0.2 which were determined in experiments. In 

non-equilibrium boundary layers with low Reynolds numbers, a^ 
experiences an additional Reynolds number correction, while and 
a^ remain the same. This correction will be explained later. 

The solution of the shear stress distribution u’v’ (y/S) is continued (143 
with a F'randtl mixing method equation. 

- u 7 * 7 1 (* • |jj> 2 (6.49) 

The mixing path £and the dissipation length L, which is required to 
determine the dissipation f (equation 6.48), together with the wall 
distance are variable values. In fully turbulent flows, JL and L are 
modelled with the relationships 

j = 0.1 • tgh (q^) • 0 (6.50) 

I ■ j 1 {t 9 h + (1 ’ fc 9 h <ij-) * f ‘ D (6.51) 

is the von Karman constant (K = 0.41) and «£. establishes the 
mixing path at the edge of the boundary layer. 


The damping factor D considers the fact that turbulent exchange 
events disappear in the laminar bottom layer. The relationship, in 
this case, depends on the pressure gradient in the free flow. In 
this calculation, it is performed with the van Driest equation 
expanded by von Cebeci C 1291s 


D = l-e 


26. Ll-11,8 p T ] 



(6.52) 
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With Z. , the distribution of the mixing path in accordance with 
equation (6.51) indicates a single-parameter context. 

To calculate the -free stream mixing path with the aid of turbulent 
kinetic energy equations is the immediate goal of this turbulence 
model. When is known, the turbulent shear stress distribution 

required for the solution of the basic equation can be determined. 
The turbulent heat flow in the boundary layer can be determined with 
the relationship for the turbulent F'randtl number (equation 6.12). 
In order to determine the free stream mixing path, the turbulent 
kinetic energy (equation 6.46) is, with the additional relationships 
in equations (6.48) through (6.52), transformed to a common 
differential equation of the 1st degree for the unknown variable 
It can then be numerically integrated with very simple 
algorithms. 

Reynolds Number Influences 

In the above form, the turbulence model is suitable only for 
computing turbulent equilibrium boundary layers with high Reynolds 
numbers. Bondary layers with low Reynolds numbers, which are the 
main focus of this dissertation, are accompanied by additional 
viscosity influences which must be considered separately. McDonald 
and his colleagues C5, 493 specified this flow situation with 

Reynolds number dependent modifications of the structural paramater 
a^ and the dissipation length L. Comparative calculations with the 
flat, unaccelerated turbulent plate boundary layer and a simplified 
analytical solution of the turbulent kinetic energy equation 
determined the following dependencies on the Reynolds number Re^ 
(see C473 ) : 


a l = V (Re 2 /Re 2,0 ),C1 + 6 » 666 ‘V (Re 2 /Re 2,0 ‘ 1)1 

L(y/S) = L^ly/fi) • [1 - exp (- 1,63 • In Re 2 + 9,7)] 


(6.53) 

(6.54) 
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The applicable values 0.0115 and 100 are applied to the unknown 
values a Q and Re^ q. L^iy/i) is in accordance with the original 
context given by equation (6.50). 


It can be seen in the above that, despite the relationships of 
equations (6.53) and (6.54), the turbulence model continues to be 
suitable only for conditions in the flat plate boundary layer 
without pressure gradient and heat transfer. Additional 
general izations are required to enable a calculation of arbitrary, 
turbulent boundary layer flows. By making substitutions for Re 0 in 
the relationships (equations 6.53 and 6.54), McDonald et al C5, 493 
introduces a more basic dependency on local boundary layer 
turbulence. The following turbulence Reynolds number was selected 
to describe boundary layer turbulence: 


R = j 

T 6 


v t dy 


lli 

ay 


(6.55) 
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To determine in what manner R^ must replace Re,, in equations (6.53) 
and (6.54), McDonald and Fish [53 used experimental studies of 
turbulent boundary layers and their own comparative calculations 
with the turbulence model. The following definitve equation was 
determined as a result: 


(6.56) 


The range between R r = 1 and R p = 40 is bridged by a polynomial of 
the 3rd degree whose coefficients are obtained from the values given 
by equation (6.56) and derivations from applicable validity limits. 

(NOTE: As a simplification to these boundary layer calculations, I 
have computed R' r with an analytical relationship (equation 6.59) 
developed by Forest C93) 
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6.5 Modelling the Boundary Laver Reversal 


6.5.1 Determination o f the Boundary Laver Intermittencv 

The transition -from a laminar into a turbulent boundary layer 
matrix can be computed with a better agreement to actual conditions, 
if this transition is not treated as an abrupt event but as a 
location with finite length. To prove this in an integral 
procedure, we must define a transient function which describes the 
physical distribution of the intermi ttency factor in the reversal. 
The intermittency factor T characterised the chronological fraction 
of turbulent flows in the transition range which alternates between 
laminar and turbulent flow. It has a zero value in laminar flow and 
the value one in fully turbulent boundary layer conditions. 
Intermediate values model the conditions in the reversal range. The 
intermittency factor T is used to determine boundary layer 
formulations in the reversal and to specify the applicable boundary 
layer turbulence. 

The matrix of the intermittency factor T is determined by the model 
equations by Forest C93. A simple differential equation applies for 
the reversal : 


dy Y e ~ v 

dx * (6.57) 


With r , this retardation equation (6.57) is based on a presumed 
equilibrium intermittency. In a quasi -stationary reversal 
situation, the boundary layer flows toward this equilibrium 
intermittency. The actual delay is determined by the factor L^. 
According to Forest C93, L depends on the boundary layer values 
Re^, and on the average effective viscosity Y of the existing 
boundary layer flow. 


L 

y 


a • Re 2 * 6 2 
(v/v) 


a = 2,0 ; 


(6.58) 


171 


(146 



v is recorded with a semi-empirical correlation which considers the 
influences of pressure gradient and free stream turbulence. 


± = X + Re, • (0,02 • y + 0,82 • c T • Tu(x)) 

V 


(6.59) 


with 


= 0,75 - B . e = /o ,0625 m 2 + y 2 

' B + 0,01 

(m = effective Polhausen parameter) (6.60) 

In this procedure, the equilibrium i ntermittency T e is identified 
with the intermittency behavior in unaccelerated flows, which is 
known from reversal tests. Referred to the standardised reversal 
coordinate z = Re^/Re^g, it has the matrix shown in Illustration 6.2 
which is described in equation (6.61) 


Rp 2s = 163 + exp(F(m) - • Tu) 



Illustration 6.2; Intermittency Matrix in the Reversal Area 
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(6.61) 


(147 


Y e = 1 + exp(l-4,4*z) ; z = Re 2 /Re 2 S 

Re,, represents the impulse loss thickness Reynolds number at 
reversal initiation. With z = 1, the reversal begins with an 

initiation intermittency of approx imatelty 1%. When the impulse 
loss thickness Reynolds number has tripled, it is regarded as 
completed with z = 3 (Illustration 6.2). This will guarantee a 
realistic, long reversal process. 

To determine the reversal initiation (more accurately: the 

intermittency initiation), Forest [93 proposes the use of suitable 
correlation functions, which were determined on an experimental 
basis. In our case, I selected a relationship developed by Abu 

Ghannam and Shaw C453 for boundary layers with pressure gradient and 
external turbulence. This relationship demonstrates good agreement 
with tests performed for this dissertation (see Chapter 5.2.1, 
Illustration 5.2.1): 


with 


F(m) = 6,91 + 12,75 m + 63,64 nt 2 for m < 0 
= 6,91 + 2,48 m - 12,27 m 2 for ra > 0 


with "m" in accordance with equation (6.60) and 

Tu = T H * I * + 2,7 /7 m 

2 

The turbulence intensity used for the determination of the reversal 
Reynolds numbers in equation (6.62) meets the criteria established 
by Abu Ghannam and Shaw C433 • and the model hypothesis by Forest 
[93. In order to record the attentuation behavior of the free 
stream turbulence (see Chapter 5.2), Abu Ghannam correlated the 
reversal initiation with average turbulence intensities between the 
running length x and the location of the boundary initiation. 
Forest also considered the influence of existing boundary layer 
intermittency. To consider the destablilizing effect in the 
reversal and after estimating the turbulence conditions in 
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the external boundary layer area, he modelled a supplement which is 
a function of J“T. According to this concept, the turbulence 
intensity Polhausen parameter (equation 6.60) is also included in 
the reversal condition (equation 6.62). The boundary layer 
intermittency and the Reynolds number for reversal initiation, which 
was determined with the correlation (equation 6.62), occur at 
locations with varying running length. The intermittency 
initiation, which is the basis for the reversal model is clearly 
ahead of the reversal Reynolds number obtained with equation (6.62) 
(see Chapter 5.5.2. 1). To balance this fundamental difference in 
new calculation programs, comprehensive development computations 
were made for the adiabatic, flat and unaccelerated boundary layer 
by using various free stream data (free stream velocity, free stream 
turbulence). The required up-stream displacement of the 
intermittency initiation was determined in such a manner that the 
applicable increase of wall tangential stress occurred at a reversal 
location for Re,,., which is described in equation (6.62). The 
result is summarized in a corrective function which is dependent on 
external turbulence. It provides the following relationship for 
the impulse loss thickness Reynolds number in the intermittency 
initiation (Re^^) required for the calculation program: 

Re 2s = Re ?s ( G1-6 . 62 ) ' (°’° 5 ‘ Xu + °’ 1 )) (6.63) 

Equations (6.47) through (6.63) provide all information required for 
determining the intermittency matrix during the laminar-turbulent 
boundary layer reversal. Intermittency has an influence on the 
turbulence model as well as on the determination of the reversal 
profiles for the velocity and temperature boundary layers. 


nb 
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If the calculation of transitional boundary layers is made for 
highly cooled surfaces, the impulse loss thickness Reynolds number 
Re£ in equations (6.58) and (6.52) is determined with the free 
stream viscosity. This incorporates the important experimental 

determination that reversal events, even with maximum wall cooling, 
must be correlated with the Reynolds numbers which refer to free 
stream material values. 

6.5.2 The Turbulence Model Purina Reversal 

It is the function of the turbulence model to define turbulent 
transport events even in the reversal area. It is, in this case, 
also assumed that, during reversal , the turbulent mixing path# 
continuously increases from zero to the value of a fully turbulent 
boundary layer. The dissipation length L increases in the same 
manner during the transition from the laminar to the turbulent 
flow. The increase of both values (j£,L) is an immediate function 
of the intermittency factor matrix. The required connection with T 
is defined for the turbulent shear stress during reversal with an 
expanded mixing path equation: 


- uV = y • (i • fy) 2 (6.64) 

This context represents an averaging of the turbulence events which, 
during reversal , occur only intermittently. In order to achieve the 
turbulent shear stress portion described in equation (6.64), the 
mixing path JLky/b) described in Chapter 6.4 is replaced with Tr 
x Z (y/S). Because mixing path and dissipation length are 
subjected to an almost identical behavior in the boundary layer, 
L(y/6) is also supplemented with J"T. The remaining relationships of 
the turbulence model, including the Reynolds number corrections, are 
used unchanged for the calculation of the turbulence values during 
reversal . 
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6.5 


Boundary Laver Profiles during Reversal 


The laminar and turbulent flow event, which chronologically 
alternate by the value of the intermittency factor T, are modelled 
by a superimposition of the laminar and turbulent profile 
formulations. This procedure of forming chronologically average 
reversal profiles was used in previous publications by Dhawan and 
Narashima C263 and Owen C253. 

The superimposition of the laminar (equations 6.25 and 6.28) and 
turbulent (equations 6.31 and 6.34) profile formulations is 
determined by the intermittency factor T according to the following 
formula: 


u uinschl ~ u lam + y u turb 
^ umschl " ^ lam’^ ~ y * + y ‘ ^ turb 


(6.65) 


Illustration 6.3 is a sketch of the profile formation described in 
equation (6.35) in the reversal area. The laminar and turbulent 
partial profiles u turb* y+ * 11 w ^^ c * 1 are shown in the 
illustrations, always fulfill identical parameters near the wall and 
at the edge of the boundary layer. The dotted line represents the 
reversal profile formed for an arbitrary intermittency factor in 
accordance with equation (6.65). 


When the reversal occurs (0 < T < 1), the velocity profile 
approaches, with increasing intermittency, the known f ul ly-turbulent 
form. An additional expansion of the range limit y£ (or y^ ) , which 
also depends on the intermittency, was required to cause 
participation of the laminar profile component u i^ n ( V’ + * * n this 
process. 
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1 _ logarithmic wall theory 

Illustration 6.3; Boundary Layer Profile during Reversal 

Equation (6.31) was supplemented in such a manner that, at the end 

of the reversal process (T = 1), that portion of the laminar partial 

profile which is near the wall and which is limited by y"' r (or y"'" ). 

s 'st 

is in accordance with the size of the laminar bottom layer of a 
turbulent boundary profile (limited by y + =6). 


+ 



+ 

lam 


u -y) + y 


6 


+ 



t lam 


(1 -y) + y 


6 


y+ s lam, yt , (6.66) 

t lain 

n.Gl . (6.31) —| 


1 “ in accordance with equation (6.31) 
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In initial calculations it was observed that, in special cases with 
reversal events -far down-stream (Tu <<) and very high u^ - and y+ 
values, inaccuracies in the model led reversal profile formation in 
the outer area of the boundary layer caused unrealistically high 
temperature and velocity gradients. In order to avoid this, the 

following limitations were introduced for the area y/6 > 0.5: 

• • 9U 3U| 3T j>T| 

y/6 > 0,5 : t < t w ; q .< % ; y y < 3y | w S y ' Sy‘| w 

6.6 Solution Procedures (151 

a’ 

The solution of differential equations, which describe the problem, 
contains the numeric integration of the basic equations for flow and 
temperature boundary layers (equations 6 . 4 , 6 . 5 , 6 . 7 , 6 . 8 ), the 

intermittency function (equation 6.67), and the determination of the 
integral balance (equation 6.46) for the turbulent kinetic energy. 

In connection with this, an evaluation of the profile equations for 
velocity and temperature matrices in the boundary layer must be 
made. This is done with the definitive equations for the Stanton 
number St and the friction coefficient c^.. 

6.6.1 Evaluation of the Profile Equations 

An evaluation of the profile equations is based on procedures used 
in integral equations. The thickness parameters 6^ and 6^ from the 

jL. %!> 

conservation equations of the flow boundary layer (equations 6.4 and 
6.5) are used for the velocity profile. As profile form parameters 
Re,, (in accordance with equation 6.39) and H,„ (in accordance with 
equation 6.6) they have a direct influence on the solution process 
and represent conditions which must be met by the profile. By 
including these form parameters, an expanded equation (6.38) is 
achieved for the velocity matrix in the boundary layer: 

= c f W * Re 2< x >> h 32( x > } <6.67) 
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Cp (for cooled boundary layers replaced by the transformation value 
c^., see Chapter 6.3.2) and the boundary layer thickness 6 (x) are 
established as implicit values due to the profile evaluation. This 
process, which occurs in stages, is as follows: 6 is determined 

with a starting value for c^., which is obtained from the solution of 
the preceeding calculation stage. The boundary layer thickness 6 is 
adjusted in such a manner that, after adding the profile with an 
integration, impulse loss thickness is in agreement with the 

setpoint value from the impulse equation. In the second stage of 
the profile adjustment an examination is made as to whether the form 
parameter of the velocity profile, which has just been 

determined, is in agreement with the setpoint value in the boundary 
layer equations. When this is not the case, c^. is adjusted, which 
is followed by another determination of the boundary layer thickness 
6 and another examination of This procedure is repeated until 

c^ and 6(x) converge in a solution. 

The temperature profile is determined in accordance with the above (152 
procedure. Therefore, I will only show the most important 

differences to the velocity profile. The temperature profile has 
the solution structure: 

T^T^ = F{ 6^fky’ R( M x )> st ’ h 6s(*)> <6.68) 

It is the definitive equation for the temperature boundary layer 
thickness S^(x) and the Stanton number St. Re^ (equation 6.45) and 
(equation 6.9) represent the form paramaters of the temperature 
profile from the enthalpy equations. 
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6.6.2 Integrating the Dif ferential Equation 

The calculation program performs a successive integration of the 
available differential equation system which consists of impulse and 
energy equations (equations 6.4 and 6.5) for the computation of the 
flow boundary layer, the two enthalpy equations (equations 6.7 and 
6.8) to determine the heat transfer and the temperature boundary 
layer, the intermittency functions (equation 6.57), and the 
turbulent kinetic energy equation (equation 6.47) to determine 
boundary layer turbulence. All basic equations are common 
differential equations of the 1st degree which can be integrated in 
a simple manner with the standard differential method. The forward 
direction is determined by the respective size of the impulse loss 
thickness. For a solution of rapid local parameter changes 
(pressure gradient, free stream turbulence, etc.) 

fix = 0.2 x 6 7 (x) (6.69) 

was selected for all equations. The integration of impulse and 
energy equations occurs in connection with the evaluation of the 
velocity profile. The same applies to the enthalpy equations and 
the determination of the temperature profile. The differential 
equations are solved in an implicit manner to guarantee numeric 
stability. 

6.7 Testina the Calculation Program with the Aid of 
Experimental Results 


This chapter describes testing of the integral procedure with 
experimental data and other tests obtained from literature. The 
procedure test is performed in stages, starting with an elementary 
analysis of the flat, unaccelerated plate boundary layer and ending 
with the flow with maximum pressure gradients and intense wall 
cooling. 
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6.7.1 Boundary Laver without Pressure Gradient 


6.7. 1.1 The Laminar-Turbulent Boundary Laver Reversal with Low 
Free Stream Turbulence Intensity 

The comparison between calculated and actual test results begins 

with the elementary case of a flat, unaccelerated, adiabatic 

boundary layer with low free stream turbulence (Tu = 1.67., u = 40 

6 o 

m/s, Re^ = 1.0 x 10 ). Similar comparative data were provided by 

Arnal C 153 and Aby Ghannam C451. The flow situation in our tests, 
where the reversal occurs throughout large plate areas, serves for 
an evaluation of the selected reversal model in the calculation 
program. Comprehensive boundary layer and traversion tests on this 
subject (see Chapter 5.3.1) indicate a three-dimensional reversal 

process which starts in the center of the tunnel and expands down- 

stream (including in diagonal direction). The Reynolds numbers at 
reversal initiation measured in the center of the tunnel is in 

agreement with data previously determined by other authors (see 
C453 ) . 

Illustrations 6.4a and 6.4b show a comparison between test and 

calculation for the boundary layer values and In this case, 

the calculation was started with an impulse loss thickness Reynolds 
number of Re,, = 140. The respective distance of the starting point 
from the leading edge of the plate was determined with the known 
analytical Blasius relationship (equation 5.8). At 20 mm it is thus 
very close to the boundary layer origination point. The tangential 
stress coefficient c^. and the respective boundary layer thickness 6 
(see profile equations, Chapter 6.3.1), which are required to 
determine the laminar starting profile, were obtained with Thwaites’ 
[1303 procedure and the given impulse loss thickness Reynolds 
number . 
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The behavior o-F the form paramater is shown in Illustration 6.4b 
for the same flow situation. According to Blasius’ theory (see e.g. 

C1081) there is an asymptopic value of 2.6 for unaccelerated laminar 
boundary layers. This value is reproduced both in the test and in 
the calculation. The calculation shows noteable deviations from the 
test only in surface 2 (x = 0.21 m). Experiment and model 

calculation are in good agreement again at reversal termination (x > 

0. 3 m) . 
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1 - test; 2 - calculation; 

3 - wall distance 

Illustration 6.5: Velocity Profiles in the Reversal Area - 

Isothermal Flow 
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Illustration 6.5 is another depiction of this test situation. It 
provides more details with the aid of velocity profiles in the 4 
measurement surfaces with x = 0.11; 0.21; 0.31 and 0.41 m. The 
laminar profile in test surface 1 and the reversal profile in test 
surface 2 are in excellent agreement with the calculation program. 
There are slight deviations with x - 0.31 m. The physical reversal 
process starts earlier in the calculation and its profile is not as 
full. The agreement is better in the subsequent test surface <x = 
0.41 m) . It is obvious that there are turbulent boundary layer 
matrices in the test as well as in the calculation. In addition, 
the calculated boundary layer thickness increase is in good 

agreement with observed data. 

6.7. 1.2 Boundary Lavers with Varying Free Stream Turbulence 

The most significant influence of external turbulence is that the 
laminar-turbulent boundary layer reversal occurs at an earlier 
point. The reasons for this is heat transfer rates which are 
considerably higher. In addition, heat transfer and tangential 

stress coefficients are slightly increased in turbulent boundary 
layers (see Chapter 5.5.2). 

Illustration 6.6 is a comparison of calculation and test where the 

turbulence intensity at the beginning of the plate experiences a 

step-by-step increase from 2.37. to 9 . 7 ’/.. The external velocity in 

this example is 55 m/s. In contrast to the above described 

experiment, a free stream temperature of 470 K was realised and the 

test plate was cooled in such a manner that there was a constant 

temperature ratio of T /T = 0.64. This value represents the 

w ® 

conditions in highly convectively cooled gas turbine blades. 
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The calculation was started with x = 20 mm, immediately behind the 
uncooled leading edge of the plate which had a length of 15 mm. As 
in all test situations, the applicable starting Reynolds number Re^ 
was determined with the Blasius relationship (equation 5.8). The 
initial velocity profile was established in accordance with the 
previous example. A Stanton number St wQ and a temperature boundary 
layer thickness 6^. was given for the initial temperature profile. 
The context for this procedure is based on analytical equations for 
flows with varying velocity and temperature boundary layer 
initiations by Kays [643 and Jischa C683. It is as follows: 


St - 0,332 • Pr' 2/3 Re '°’ 5 [1 - (i-) 0,75 ]' 1/3 


r 13 1 n ,c ,3/4, -,-1/3 

= [ T5 * Tr t 1 ‘ ) )J 


(6.70) 


( ^ = uncooled starting length) 

( x o = starting point for the calculation) 

The model calculation provides a good description of the advance of 
the reversal as a function of turbulence intensity. The reversal 
length is also in good agreement with the test. With maximum 
turbulence intensity (Tu q = 8.77.), the Stanton numbers in the 

reversal increase starting with an impulse loss thickness Reynolds 
number of Re^ = 135. The early reversal thus starts below the 

limiting Reynolds number (indifference Reynolds number) for reversal 
initiation, which is known from linear stability theory and which - 
according to available literature (see C22, 453) - has a value of 
Re ? ~ 160. This early reversal is caused by the Reynolds number 

correction for the intermittency initiation, which was introduced 
with equation (6.63). It is an important contributor to a good 
reproduction of the test data. Test and calculation are in very 
good agreement in all cases in the turbulent range. 
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Tests which, with the aid of heat transfer distributions, show the 
influence of free stream turbulence on the reversal (even with 
isothermal conditions) were recently conducted by Blair and Werle 
[373 on a slightly heated plate with constant heat flow (q w <x) = 
const.). Illustration 6.7 compares the data with model calculations 
and shows that a good agreement was achieved. With high turbulence 
intensities of 6.2’/. and 7.5'/., the laminar -turbulent reversal already 
occurs within the unheated leading edge of the plate, which was 
attached for design reasons. It was not immediately available to 
the calculation, because the procedure for both uncooled and 
unheated starting lengths is designed for the flow immediately 
behind this plate. Even in the above cases with high free stream 
turbulence, the calculation procedure was started in a laminar 
manner. This is the reason why, in the illustrations of the 
calculation results, reversal events are still shown. However, the 
good agreement with the data for low free stream turbulence confirms 
that the integral procedure provides an accurate reversal 
description. 
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1 - tests, Blair and Werle (1980) 

2 - calculation 

Illustration 6.7; Stanton Number Matrices on the Flat, 

Unaccelerated Plate with Varying Free Stream 
Turbulence and a Slightly Heated Wall 
(Tests by Blair and Werle C373) 

Illustration 6.8 shows the significance of external turbulence for 
the boundary layer parameters tangential stress, impulse loss 
thickness and form parameter for various turbulence intensities. In 
these tests, wall and free stream temperatures had the same level 
(no heat transfer, flow condition 2 in Table 5.1). The comparison 



between test and calculation refers only to turbulent boundary layer 
flows. Because reversal behavior was still present in test surface 1 
<x = 0.12 m) with a minimum free stream velocity of (Tu = 2.37.), 
the re-calculation was begun at the subsequent test location (with x 
= 0.22 m) . In all cases, the program was initiated in a turbulent 
area; the values for Re^ and H„ were derived from the tests. An 
adjustment of the profile equation at this point determined the 
additional start values for c^. and H^. A distribution proposed by 
Voges <[1313, see l 43), which also records the influence of external 
turbulence, was selected as a starting profile for the boundary 
layer turbulence. This given profile was used to solve the integral 
turbulence equation (equation 6.46) and to determine the mixing path 
distribution at the starting location. 

Illustration 6.8a shows the measured and calculated c r matrices. 
The effect of tangential stress increase induced by free stream 
turbulence is also present in the calculation. The deviation from 
the test data is only minor (maximum deviation approx. 67.). The 
applicable impulse loss thicknesses are re-cal cul ated in 
Illustration 6.8b. A special situation occurs in this test: an 
increasing free stream turbulence results in a deformation of the 
boundary layer profiles without a noticeable change of the local 
impulse loss thickness (Illustration 6.8, see Illustration 5.12, 
Chapter 5.3.2). The calculation describes a general situation where 
increasing turbulence results in an increase of the impulse loss 
thickness ( II lustration 6.8b). However, as shown by Illustration 
6.8b, the resulting differences are only minor. 

A mirror-image trend of the tangential stress coefficient occurs in 
the form parameters (Illustration 6.8c). Increasing external 
turbulence generates fuller boundary layer profiles which result in 
a decrease of the form parameters. This behavior is accurately 
reproduced by the calculation procedure. When calculating with 
minimum free stream turbulence (Tu q = 2.3'/.) , a build-up process is 
observed during which the form parameter increases slightly. 


1 
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The so-called wake portion, i.e. the deviation -from the logarithmic 
wall theory in the external boundary layer sane, decreases with 
increased -free stream turbulence, while the validity of the wall 
theory remains unchanged. Both effects are recorded by the integral 
procedure in good agreement with the test data. 



1 - test; 2 - calculation 

Illustration 6,9; Dimensionless Turbulent Velocity Profiles 

during Varying Free Stream Turbulence 

6.7. 1.3 Transitional Boundary Layers with Wall Cool i no 

According to experimental findings, wall cooling results in an 
increase of heat transfer and wall tangential stress. At the same 
time, there are changes in the velocity and temperature matrices 
within the boundary layer which can be attributed to fuller 
profiles. Contrary to the rules of linear stability theory, I did 
not, in my tests, discover a stabilization and delay of the reversal 
due to wall cooling, when free stream turbulence was continously 
present (see Chapter 5.52). 
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Illustration 6.10 is an example of three test series conducted with 
varying temperature conditions T^/T^ and the applicable calculations 
■for the relatively low turbulence intensity of 2.3'/.. 
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Illustration 6.10; Heat Transfer and Boundary Layer Development 

with Varying Wall Cooling and Low Free Stream 
Turbulence (dp/dx = 0) 

The flow velocity was variable in the tests, so that a variation of 
the Reynolds number is an additional factor. In the case of maximum 
wall cooling (T /T = 0.53) and minimum flow velocity (li = 38 m/s), 
a laminar flow over wide areas of the plate occurs in the tests. 
Starting with x = 0.2 m, this flow, with slowly increasing Stanton 
numbers, is replaced by a reversal . This type of reversal behavior 
is reflected in the calculations; however, the transition to 
turbulent flow is more rapid than in the experiment. There are also 
differences between test data and calculations in the laminar 
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boundary layer section. The discrepancy is caused by calculated 

heat transfer coefficients, which are too low. This is caused by 

the applied material value transf ormation for the determination of 

velocity and temperature profiles during wall cooling; it is 

insufficient in this case. This was to be expected, because the 

tests showed a higher increase of the laminar heat transfer with 

wall cooling than is known from literature, and because it is 

predicted by the material value transf ormations implemented in the 

calculation program (see Chapter 5.5.1). The differences between 

test and calculation data in the laminar boundary layer area are a 

function of the intensity of the wall cooling. This is confirmed by 

the lesser discrepancies during reduced wall cooling, which is 

identical to an increased temperature ratio T /T . 

w « 

In Illustration 6.10, the increase of external velocity from 38 to 

55 m/s and the reduction of wall cooling from the value T /T = 0.53 

w « 

(which is hardly passible under actual conditions) to T T = 0.64, 

w « 

which is a realistic condition in turbines, shows a much earlier 
boundary layer reversal . This is confirmed by the model with 
respect to reversal initiation as well as to the length of the 
reversal. There is very good agreement in the turbulent zone. 

In comparison to the above test, a further reduction of wall coaling 

to T /T = 0.8 and external velocity (from 55 to 47 m/s) - as shown 
w w 

in the lower half of the illustration - causes location and length 
of the laminar-turbulent reversal to remain almost unchanged. This 
is verified by the calculation. In this case, the accuracy of the 
calculation is especially good in all boundary layer areas. 

Illustration 6.11 compares tests with a turbulence intensity of 
8.7%, which is a typical condition in turbines, to their respective 
calculations. Due to the high turbulence, the reversal - in all 
temperature conditions - occurs very close to the leading edge of 
the plate. This is also shown in calculations. Identical results 
for the distribution of heat transfer in turbulent boundary layer 
zones are achieved by both the test and the calculation. The 
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influence of wall coaling in turbulent boundary layers is thus 
verified by applied profile transf ormations in a very accurate 
manner. 

The tangential stress coefficients for the heat transfer 
coefficients during high free stream turbulence and variable wall 
cooling, as described in Illustration 6.11, and for an additional 
adiabatic situation (flow conditions 1.1, 3, 4 from Table 5.1) are 
shown in Illustration 6.12. The decreasing Reynolds number is a 
contributing factor to the increase of Cj.-values with increasing 
wall cooling. The calculated matrices are in agreement with the 
test results. In all cases, variations are less than 107.. 
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0 0.1 0.2 0.3 x [ml 0.4 


1 - test; 2 - 

Illustration 6.11; Heat Transfer and 

with Varying Wall 
Turbulence (dp/d>: 


calculation 

Boundary Layer Development 
Cooling and High Free Stream 
= 0 ) 
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The velocity and temperature pro-files shown in Illustrations 6.15 


and 6.14 must be assigned to the c ^-values in test surface 4 (;•: = 

0.42 m) . Illustration 6.13 shows the influence of wall cooling on 
the logarithmic wall theory. Due to varying density and viscosity 
as a function of the wall distance, the area of the wall theory 
obtains greater y + - and u + - values if, as shown in Illustration 
6.13, these coordinates are derived from material data at the wall. 

This event is simulated by the calculation program in an almost 
perfect manner, which is proof of the practicality of applied 
profile transformations. 
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1 - test; 2 - calculation 

Illustration 6.12; Resistance Coefficients in Unaccelerated Flows 

with Varying Wall Cooling 
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1 - test; 2 - calculation 

Illustration 6.13: Dimensionless Turbulent 'velocity Profiles 

with and without Wall Cooling 

The only differences (max. 10 %) exist in the wake area. This flow 

area is mainly determined by the material data at the edge of the 

boundary layer. This is not considered in the material value 

determination of the entire boundary layer, which was introduced by 

the transf ormation in the profile equations. Because of this, 

greater deviations are shown on the edges of the boundary layer. 

Another factor, which contributes to an over-emphasis of the wake 

area, are the friction coefficients c^ , which are calculated too low 

by up to 10% (see Illustration 6.12, x = 0.42 m) . Because of this, 

excessive edge values u+ are used for calculating the profile matrix 

(see profile equations, Chapter 6.3). This applies even in 

situations without wall cooling (T /T = 0.42 m). 

3 w « 
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The two velocity matrices -for cooled situations are shown again in 
Illustration 6.14a, where the wall distance was not standardized. 
Despite the deviations in the boundary layer’s external area 

observed in Illustration 6.13, which is a very sensitive semi- 

logarithmic illustration, turbulent pro-files are in good agreement 
with the tests. 

The respective temperature matrices are shown in Illustration 6.14b 
in a dimensionless -form and as a function of velocity. The 

calculated matrices are in agreement with a linear distribution 
which, with an incompressible flow (T /T^ •* 1.0, Ma -» 0) and a 

F'randtl number of 1, forms an asymptote. The measured deviation 
from this ideal temperature theory is not substantiated by the 
integral procedure with given profile equations. However, as shown 
by the applicable Stanton numbers in Illustration 6.11, this 

imperfection does not have an appreciable negative affect on the 
quality of calculated heat transfer. 









6.7.2 


Flows with Pressure Gradients 
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In my experiment , I studied boundary layers under the influence of 
two different pressure gradient matrices, which simulate the 
circulation on the pressure and vacuum sides of a turbine blade. 
The various pressure gradient matrices were superimposed with 
variations of the Reynolds number, free stream turbulence and wall 
cooling. In the following paragraphs, integral procedures will be 
used to re-calculate the characteristic flow situations. 

6. 7.2.1 Comparative Calculations for Constant Free Stream 
Acceleration 

Illustration 6.15 provides a summary of the test situations with 
almost constant free stream acceleration <k(x) = const.). It shows 
distributions obtained with contour 1 of the free stream velocity 
u^tx) and the acceleration parameter k(x) for flow conditions with 
varying wall cooling <T /T = 0.64 and 0.53), the acceleration 

intensity <k = 1.2 x 10 ; 2.3 x 10 and 2,5 x 10 ) and the 

Reynolds number of the plate's running length which, in Illustration 
6.15, was determined at the flowoff conditions at the plate end. 
Additional information on these flows can be found in Table 5.1 
(flow conditions 6, 7, 8). As in all other re-calculations, the 

turbulence intensity matrix in the free flow is given by the test 
data (see Chapter 5.2) 

The measured Stanton numbers for Re = 1.7 x 10 6 are, for various 

turbulence intensities, compared to model calculations in 

Illustration 6.16. In all cases, the calculation starts with x = 

0.02 with an initial Reynolds number of Re,, = 112. Using the 

Blasius relationship for laminar boundary layers without pressure 

gradients (equation 5.8), initial Reynolds numbers were also 

approximated in accelerated flow conditions. The best agreement is 

shown in Illustration 6.16 for the test situation with maximum free 

stream turbulence, Tu = 117.. Despite acceleration (k = 1.2 x 10 ^) 

o 
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the turbulence causes the boundary layer reversal to occur very 
early, which is simulated by the calculation. There is also a very 
good agreement between test and calculation data in the turbulent 
boundary layer section. The calculation predicts a short laminar 
boundary layer matrix with an inlet turbulence of 47., while the test 
already shows a reversal flow at that location (see Illustration 
5.46a). The differences between test and calculation are, as usual, 
minor in the fallowing boundary layer area. In tests with a 
turbulence intensity of only 2.37., the reversal to turbulent flow 
occurs very slowly and seems to be completed at the plate end. The 
reason for this is the superimposed acceleration. 
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of Re-Calculated Test Cases with Contour 1 
(Contour on the ’Pressure Side’)' 
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In this case, the calculation predicts a reversal initiation -further 
down-stream, which is -followed by a relatively -fast reversal 
process. This is the reason- -for deviations between test and 
calculation results in the area between x = 0.15 m and 0.26 m. 

Illustration 6.17 shows analogous tests with higher acceleration (or 
smaller Reynolds numbers, see Illustration 5.46b) with k=2.3 x 10 
The more intense acceleration causes a stronger atten nation of the 
turbulence movement in the boundary layer. Because of this, there 
is no reversal with minimum external turbulence intensity. The 
calculation is in excellent agreement with this behavior. With a 
free stream turbulence of 47., a balance is shown in the test between 
the destabilizing effect of external turbulence and the stabilizing 
influence of free stream acceleration, which results in constant 
Stanton numbers throughout large areas of the plate. In contrast to 
this, the calculation model indicates a discrete reversal process. 
This causes deviations from the measured Stanton numbers of up to 
157.. 


Turbulent boundary layer flows are again present throughout the 

entire test area with maximum free stream turbulence, Tu = 117.. 

o 

This boundary layer behavior is described by the integral procedure 
with good accuracy. 
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i - test; 2 - calculation 

Illustration 6.16; Stanton Number Matrix in Accelerated Flows 

(Contour on the ’Pressure Side’, high 
Reynolds Number) 

A situation similar to Illustration 6.17, with identical -flow 

turbulence intensities is shown in Illustration 6.18 . The Reynolds 

number - and therefore the acceleration intensity k(x) matrix - is, 

compared to the above described situation, almost identical (see 

Illustration 6.15). The only difference is the increased coaling 

(T /T = 0.53 versus 0.64). With this flow parameter constellation 

there is no boundary layer or reversal behavior change compared to 

the test with lower wall cooling (T /T = 0.64; Illustration 6.17). 

w « 

This phemonenon was discussed in detail in Chapter 5.5.2 in the 
analysis of test results. The reason for it is that the wall 

cooling does not result in any noticeable boundary layer 

stabilization. The calculation is in accordance with this reversal 
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theory. All other comparative calculations show the same reversal 

behavior as with low cooling intensity <T w /T w = 0.64; Illustration 

6.17). The best agreement between test and calculation is again 

with minimum and maximum -free stream turbulence (Tu = 2.3 and 11%), 

o 

while with medium -flow turbulence of Tu = 4.0%, the above 

o 

weaknesses of the reversal calculation can again be observed. 



1 - test; 2 - calculation 

Illustration 6.17; Stanton Number Matrix in Accelerated Flows 

(Contour on the ’Pressure Side’, Low 
Reynolds Number) 
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1 - test; 2 - calculation 

illustration 6.18; Stanton Number Matrix in Accelerated Flows 

with High Wall Cooling (Contour 1, Low 
Reynolds Number) 

6.7. 2.2 Comparative Calculations -for Flows Accelerated on 
the Vacuum Side 

The distribution o-f the external velocity and the acceleration 
parameters k(x) o-f this example are summarised in Illustration 
6. 19. Again, we must distinguish between three situations: a test 
with a high Reynolds number (Re^ =2.1 x 10 6 ) and a temperature 
ratio of T /T = 0.8; another test with T /T = 0.8 but reduced 

w ot W » 

Reynolds number (Re = 1.0 x 10 ); and an experiment with the lowest 

L & 

Reynolds number Re. = 0.93 x 10 and intense wall cooling of T T 

L w w 

= 0,53 (see Table 5.1, test numbers 9, 12, 13). 

1 C’ 
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Freistromgeschwindigk eit 



J'*Platlenlauflange x[m] 


1 - free stream velocity 

2 - running length of the plate 

3 - acceleration parameter 

Illustration 6.19: Velocity and Acceleration Parameter Matrices 

of Re-Calculated Test Cases with Contour 2 
(’Vacuum Side’) 
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In the last two cases, the relaminization criteria of k z 3 x 10 

is clearly exceeded by the acceleration intensities k =5.7 x 10 ^ 

max 

and k = 6.3 x 10 . Such intensive acceleration conditions are 

max 

very typical for actual gas turbine blade circulation. They place 
highest demands on the calculation program. 


Illustration 6.20 compares model calculations with high Reynolds 
numbers to tests where the flow turbulence intensity of 2.37. was 
increased to 107.. The alternating effects of pressure gradients 
with transitional events are described in an accurate manner. The 
reversal initiation is accurately shown with minimum turbulence 
intensity, Tli q = 2.37., while the following reversal process occurs 

more rapidly in the test. With 4.07. flow turbulence, differences 
exist only in the forward plate section (x < 0.2 m) . Again, they 
can be attributed to the fact that the test shows reversal behavior 
of the boundary layer starting with the beginning of the plate, 
while the laminar calculation shows the reversal at a later point, 
after a distinct laminar initiation. 



1 - test; 2 - calculation 

Illustration 6.20: Stanton Number Matrix in Highly Accelerated 

Flows (Contour on the ’Vacuum Side’ High 
Reynolds Number) 
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1 - test; 2 - calculation 

Illustration 6.21: Stanton Number Matrix with Highly Accelerated 

Flow (Contour on the ’Vacuum Side’, Low 
Reynolds Number) 

In the cases with maximum -flow turbulence (Tu = 77. and 107.), there 

o 

is satisfactory agreement with the test. However, compared to 
previous cases with and without constant acceleration, there are 
certain inaccuracies in the turbulent flow areas. There are maximum 
deviations of 10-157.. 



In cases of maximum acceleration with k = 5.7 x 10 ^ , 

max 

Illustration 6.21, even more distinct effects of the pressure 
gradient can be observed. With the lowest turbulence intensity of 
2.37., the turbulence in the accelerated flow area is attenuated so 
much that, even in the delayed area (x > 0.26 m) , the reversal is 
only gradual. The reversal initiation, which is located near the 
end of the plate, is accurately reproduced by the calculation 
model. However, the reversal to the turbulent condition occurs 
somewhat earlier in the calculation. A similar behavior is observed 
with a higher turbulence intensity of 4.07.. The test with the 
highest free stream turbulence, Tu q = 107., represents the very 

complicated flow situation where, starting with the beginning of the 
plate, there is a quasi-stable laminar-turbulent intermediate 
condition during the entire acceleration phase. In contrast to the 
laminar situation, the heat transfer is, in this case, higher by 100 
percent or more. This boundary layer behavior was discussed and 
analyzed in detail in Chapter 5.6 and Illustration 5.47b. Due to a 
lack of accurate initiation profiles (this starting location was not 
permissible for profile tests), the applicable calculation was begun 
in a laminar manner, i.e. contrary to the experimental result. The 
same applies to all other cases with low free stream turbulence in 
Illustration 6.21. The laminar boundary layer initiation of the 
calculation is followed by a discrete reversal area, which is 
completed after approximately 100 mm of plate running length and 
which, down-stream of this location, results in heat transfer 
coefficients which are exclusively turbulent. The discrepancy 
between test and calculation is up to 507. in the accelerated flow 
section. However, a careful adjustment of the calculation (6.57), 
which describes the i ntermi ttency matrix, to the given flow 
behavior, will result in better agreement. Because no general 
criteria could be found for the reproduction of these complex 
boundary layer events (this requires comprehensive boundary layer 
tests), I did not perform such an expansion of the calculation 
program. It should be noted, however, that the connection between 
turbulence model and i ntermi ttency function, which was used in the 
calculation procedure, is the best method for carrying out such a 
procedure. 
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The spirit of the comparison described in Illustration 6.21 can be 
applied to Illustration 6.22 . In Illustration 6.22, the calculation 
program is examined -for an identical Reynolds number, and thus the 
same acceleration conditions, but with -far more intensive wall 
cooling (T /T = 0.53 versus T /T = 0.8). As expected, this final 
example is in agreement with the test. In addition, the calculation 
shows the same boundary layer behavior as in the previous situation 
with low wall cooling (Illustration 6.21). The boundary layer 
behavior which, as shown by the calculation, confirms again that a 
modelling of the reversal with wall coaling and external turbulence 
is only performed correctly when, as in our program, the reversal 
Reynolds number, which causes the laminar-turbulent boundary layer 
transition, is applied to the free stream viscosity and not to a 
viscosity computed from a reference temperature. 
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1 - test; 2 - calculation 

Illustration 6.22; Heat Transfer Coefficients in Highly Accelerated 

and Highly Cooled Boundary Layer Flow 
(Contour 2) 
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This chapter described the calculation of transitional boundary 
layers, both with and without pressure gradients with varying free 
stream turbulence and varying wall cooling. The examples introduced 
above showed that the influences of these parameters are accurately 
recorded by the newly developed integral, procedure. The best 
agreement between test and calculation is always in turbulent 
boundary layer sections. There are differences in the prediction of 
the reversal length, if the flow has pressure gradients 
(acceleration) and low turbulence intensities. In this case, a more 
careful adjustment of the reversal model should result in better 
agreement. In some cases with pressure gradient and increased flow 
turbulence, differences occurred with heat transfer coefficients at 
the beginning of the plate. They can, on the one hand, be 
attributed to the more rapid reversal behavior in the calculation 
model, and, on the other hand, to different initiation conditions. 
An initial laminar boundary layer was assumed in the calculation, 
while reversal conditions were already present in the tests. A more 
exact starting condition for the calculation was not possible, 
because the flow in this area was not accessible. Pertinent 
information on this subject is not available in literature. With 
similar flow situations, the calculation of cooled blade boundary 
layers presents identical starting value problems. Therefore, it 
should be the goal of further studies to analyze the velocity and 
turbulence structure in initial boundary layers, which are 
influenced by these external turbulences and pressure gradients. 

The remaining inaccuracies of the calculation program signify the 
present level of numeric boundary layer calculation. They should 
not be attributed to the selection of the calculation method - 
integral or differential procedures - but to the unreal i ability of 
present turbulence models. Even highly developed k-< turbulence 
models, which were incorporated in differential procedures, are 
unable to give a better description of long reversal processes. 
This was clearly demonstrated by parallel calculations for the test 
situations in this chapter, i.e. the Lam-Bremhorst-k-f turbulence 
model, which requires a great deal more effort (see C4, 1103). With 
an almost identical quality of the results in laminar and turbulent 
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boundary layer sections, the reversal length was predicted even more 
accurately with the simple integral procedure and its turbulence 
model . 


7. SUHMftRY 
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This dissertation discussed the influences of -free stream 
turbulence, wall cooling and pressure gradients on the heat transfer 
and the boundary layer development in cooled surfaces. In order to 
determine the significance which is to be attributed to these 
parameters, individually or when superimposed, comprehensive 
experimental and theoretical studies were conducted with 
transitional boundary layers. The reasons for this were a lack of 
basic knowledge and insufficient data on heat transfer and boundary 
layer behavior in cooled gas turbine blades. 

The experiments, which consisted of an analysis of local heat 
transfer and a study of applicable boundary layer profiles, were 
performed in an even boundary layer tunnel. Free stream turbulence 
intensities varied between 1.6 and 11 In order to determine the 
influence of pressure gradients, accelerated flows were generated. 
Expressed by the acceleration paramater k, they exist in the blade 
circulation of gas turbines with: 

„ du 

- 1.2 - 6.3 • ID' 6 


To simulate realistic turbine conditions, a step-by-step increase of 

the free stream temperature was performed, which resulted in the 

temperature ratios T /T = 1.0; 0.8; 0.64; and 0.53. 

w w 

The test results confirm and supplement effects as they are known 
from previous studies. These studies were conducted with mostly low 
free stream turbulence, low -free stream acceleration and without 
consideration to wall cooling. Based on this and by including 
turbine-like flow conditions, additional problems were solved which 
are important for heat transfer and boundary layer 
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development in cooled gas turbine blades. This includes the 
•following: the influence of wall cooling on the laminar-turbulent 
reversal process; the significance of wall cooling in turbulent 
boundary layers with free stream turbulence; and boundary layer 
behavior in flows with very intensive acceleration and high external 
turbulence. 

The tests show that, in boundary layers with free stream turbulence, 
even the most intensive wall cooling will not result in a change of 
reversal behavior. Consequently, the influence of external 
turbulence dominates over the effect of wall cooling, which 
stabilizes the boundary layer. 

Heat transfer and boundary layer tests further revealed that, in 
turbulent boundary layers, the effects of free stream turbulence are 
not dependent on the wall cooling. The deformation of boundary 
layer profiles caused by external turbulence and the resulting 
increase of wall tangential stress and heat transfer are identical 
with events at the uncooled wall. 

Without changing the effects of free stream turbulence, the wall 
cooling results in a deformation of the boundary layer profiles and 
an increase of the heat transfer. The test with maximum cooling in 
laminar boundary layers resulted in a heat tranfer increase of up to 
20’/.. Under the same conditions, local heat transfer coefficients 
increase by only about 10-12’/. in the area of the turbulent boundary 
layer. 

A study of the influences of negative pressure gradients revealed 
that, in a combination of high free stream turbulence and high free 
stream acceleration - which is typical for gas turbine blades - 
there are flow situations where, over large areas, the boundary 
layer obtains an almost stable laminar-turbulent intermediate 
condition. In contrast to laminar flow, these conditions result in 
drastic increases in the heat transfer (up to 1007.). A closer 
analysis of this event showed that the occurrenceof stable laminar- 
turbulent intermediate conditions is a function of high free stream 
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intensities, intensive pressure gradients and the acceleration path 
along the wall. Because of its complexity and the present level of 
available information, this boundary layer flow cannot be accurately 
described and characterized. Further boundary layer studies are 
therefore recommended. 

In order to reconstruct experimental results in a theoretical 
manner, a numerical calculation procedure was developed which is 
based on an expanded integral method. As with differential 
procedures, this new procedure conducts a separate calculation of 
flow and temperature boundary layers. This eliminates the use of 
analogy equations for the determination of heat transfer, which is 
common practice with integral procedures. The velocity and 
temperature distributions in the boundary layer are established with 
suitable profile equations. These boundary layer matrices are the 
only flexibility limitation as compared to field methods. A single 
equation turbulence model, which was further developed in this 
dissertation, was used - in connection with a mixing path equation - 
to describe the turbulent exchange events in the boundary layer. 
The results of the numeric calculation are in good agreement with 
experimental data. Re-calculations show that the influences of free 
stream turbulence, wall cooling and pressure gradients are 
accurately represented. The goal of assigning the efficieny of 
field methods to integral procedures, which are required for an 
accurate prediction of the heat transfer in cooled gas turbine 
blades, was thereby achieved. 

The best agreement between calculation and test is in flows with 
high external turbulence. Reversal initiation and reversal length 
are very accurately described in most cases. This is mainly due to 
the expanded reversal model. There are deviations in the 
calculation of the reversal length with accelerated flow and a low 
free stream turbulence. Under these flow conditions, the reversal 
process is faster than in the test. An improvement of the agreement 
between calculation and test requires additional fine-tuning of 
reversal and turbulence modelling. 
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In summary, it can be stated that the above studies resulted in 
discoveries on the influence of free stream turbulence, pressure 
gradients and wall coaling, which were unknown before and which 
expand our understanding of heat transfer and boundary layer 
behavior on cooled surfaces, especially on cooled gas turbine 
blades. Newly developed calculation procedures provide theoretical 
support for experimental findings. The results can be used in the 
design of cooled gas turbine blades. 
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